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Section  1 


INTRODUCTION  AND  SUMMARY 

Q- switched  laser  pulses  have  found  wide  applications  in  diverse  areas  of  pure 
and  applied  research.  Some  fields  which  have  been  opened  to  investigation  by  the 
availability  of  such  pulses  include  ga”  breakdown  at  optical  frequencies,  plasma 
production  for  thermonuclear  research,  optical  harmonic  generation  and  parametric 
amplification,  stimulated  scattering  effects  and  coherent  propagation  effects.  Q- 
switeched  pulses  have  also  been  used  for  ranging  and  guidance  systems,  high  3peed 
photography,  medical  research,  precision  machining  and  other  applications.  The  min¬ 
imum  pulse  durections  obtainable  with  the  various  existing  Q-switching  techniques  are 
limited  to  approximately  10"  seconds  because  of  the  time  required  for  the  buildup 
of  the  pulse  in  the  laser  cavity. 

The  advent  of  the  mode-locked  laser  has  brought  abouc  a  four  to  five  order  of 
magnitude  decrease  in  pulse  duration,  from  10"®  sec  to  less  than  lO"1^  sec  and  a  three 
order  of  magnitude  increase  in  peak  power,  from  109  watts  to  10^  watts  or  greater. 
Detection  systems  capable  of  resolving  events  on  a  time  scale  as  short  as  iO"‘L®  sec 
have  been  developed.  It  is  expected  that  these  developments  should  open  more  areas 
of  investigation  and  should  lead  to  additional  scientific,  military  and  commercial 
applications.  The  short  time  duration  of  the  pulses  allows  the  investigation  of 
atomic  and  molecular  processes  and  coherent  interaction  effects  on  a  time  scale  that 
was  previously  inaccessible  to  direct  observation  and  high  power  available  with 
these  pulses  allows  the  investigation  of  nonlinear  effects  that  were  previously 
unobservable. 

The  United  Aircraft  Research  Laboratories  have  been  conducting  under  the  present 
contract  a  continuing  investigation  of  short  time  duration  laser  pulses  and  their 
interactions  with  matter.  This  investigation,  while  primarily  concerned  with  the 
picosecond  duration  pulses  produced  by  mode-locking,  has  also  been  extended  to  tran¬ 
sient  phenomena  involving  longer  pulses.  During  the  period  covered  by  this  report, 
work  has  been  conducted  in  the  following  areas:  (a)  the  analysis  of  the  propagation 
of  ultrashort  pulses  in  a  resonant  medium,  (b)  the  analysis  of  the  generation  of 
Cerenkov-like  radiation  by  nonlinear  optical  effects,  (c)  the  development  and 
demonstration  of  a  novel  technique  for  the  measurement  of  a  phase  structure  of 
picosecond  pulses,  (d)  the  study  of  stimulated  transient  Raman  Scattering  in  liquids 
and  gases, (e)  the  development  of  extremely  fast  pumping  sources  for  organic  dye 
lasers,  (f)  the  mode  locking  of  a  flashlamp  pumped  organic  dye  laser,  (g)  the  demon¬ 
stration  of  traveling  wave  laser  action  in  an  organic  dye,  (h)  the  investigation  of 
the  possibility  of  obtaining  stimulated  emission  from  laser  produced  plasmas,  (i) 
the  investigation  of  nonlinear  polarization  effects  in  anisotropic  molecular  liquids. 

In  the  area  of  the  analysis  of  pulse  propagation,  the  major  portion  of  time  has 
been  directed  toward  the  formation  of  a  unified  synthesis  of  advances  thus  far 
achieved  by  various  workers  in  obtaining  analytical  results  in  the  field  of  ultra- 
short  optical  pulse  propagation.  In  the  course  of  reformulating  the  results  of  these 
workers,  a  number  of  minor  extensions  and  simplifications  of  their  work  have  been 
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obtained.  Although  the  results  of  some  investigators  have  not  been  incorporated  at 
this  time,  it  is  expected  that  the  material  reported  herein  will  constitute  a  major 
step  toward  the  compilation  of  a  complete  summary  of  the  current  status  of  this  field. 

During  this  period,  consideration  was  given  to  the  "dc"  component  of  the  nonlinear 
polarization  induced  in  a  dielectric  by  the  passage  of  a  high  power  optical  pulse. 

This  component  constitutes  a  polarization  source  that  moves  through  the  dielectric  with 
the  group  velocity  of  the  optical  pulse.  This  polarization  can  lead  to  Cerenkov  radia¬ 
tion  at  frequencies  whose  velocity  of  propagation  is  less  than  that  of  the  pulse.  This 
effect  is  of  fundamental  interest  in  itself  as  it  represents  a  completely  new  technique 
for  the  generation  of  Cerenkov  radiation.  It  could  have  application  to  millimeter  and 
submill imeter  wave  generation  and  noise  sources.  Of  even  more  general  importance  is 
the  realization  that  many  experiments  in  beam-wave  interactions  that  can  be  performed 
with  difficulty  using  megavolt  electron  beams  can  be  performed  easily  using  the  non¬ 
linear  polarization  associated  with  a  propagating  picosecond  pulse. 

Since  the  first  measurement  of  the  time  duration  of  the  pulses  produced  by  mode¬ 
locking  of  neodymium-glass  lasers,  it  has  been  apparent  that  the  duration  of  the  pulses 
are  considerably  in  excess  of  the  minimum  possible  duration  as  determined  from  the 
reciprocal  of  the  bandwidth.  A  time-bandwidth  product  in  excess  of  unity  implies 
a  phase  or  amplitude  structure  in  the  pulses.  The  first  measurement  of  this  structure 
was  obtained  at  this  laboratory  in  1968.  It  was  found  that  the  pulses  could  be  com¬ 
pressed  in  time  by  passing  them  through  a  dispersive  optical  system.  The  dispersion 
of  the  system  that  was  capable  of  compressing  the  pulses  was  such  that  the  transit 
time  increased  with  increasing  wavelength.  The  compression  of  the  pulses  indicated 
the  presence  of  a  linear  component  of  the  frequency  vs  time  of  the  pulses,  with 
longer  wavelengths  or  lower  frequencies  coming  earlier.  The  compression  also 
indicated  that  a  significant  fraction  of  the  observed  bandwidth  was  due  to  the  fre¬ 
quency  sweep.  While  the  compression  experiments  serve  to  demonstrate  the  presence 
of  a  frequency  sweep,  it  is  difficult  to  extract  more  detailed  information  about 
the  phase  structure  from  them.  During  this  reporting  period  a  technique  has  been 
devised  to  extract  the  maximum  available  information  on  the  phase  structure  of  the 
pulses  subject  to  the  constraint  of  present  detectors.  The  information  that  can  be 
obtained  from  this  technique  consists  of  essentially  /U  resolvable  points  on  a  fre¬ 
quency  vs  time  curve.  Here  N  represents  the  time-bandwidth  product  of  the  pulse. 

Initial  experimental  results  have  been  obtained  and  indicate  that  the  technique 
works  as  designed.  These  results  also  are  consistent  with  the  compression  experiments 
and  reconfirm  the  presence  of  a  linear  component  of  the  frequency  sweep. 

Experiments  in  stimulated  scattering  with  picosecond  pulses  have  continued. 

During  this  period  emphasis  has  been  placed  on  the  investigation  of  stimulated  Raman 
scattering  in  liquids  and  gases,  and  a  number  of  new  results  have  been  obtained.  The 
essential  difference  between  stimulated  Raman  scattering  with  nanosecond  duration 
pulses  from  a  Q- switched  laser  and  that  produced  with  picosecond  duration  pulse  from 
a  mode-locked  laser  lies  in  the  finite  build-up  time  of  the  associated  phonons.  In 
the  case  of  picosecond  pulse  excitation,  the  phonon  population  does  not  reach  its 
steady  state  value  and  the  transient  nature  of  the  scattering  must  be  considered.  An 
important  consequence  of  the  transient  nature  is  the  fact  that  in  the  extreme  tran¬ 
sient  limit,  the  Raman  gain  is  determined  not  by  the  peak  value  of  the  spontaneous 
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Raman  cross  section  but  by  the  total  area  under  the  curve  of  cross  section  vs  fre¬ 
quency.  It  is  possible  therefore  to  excite  Raman  lines  in  the  transient  limit 
that  cannot  be  excited  in  the  steady  state  due  to  competition  with  other  narrower 
lines  that  have  a  higher  peak  cross  section  but  lower  integrated  cross  section. 

Using  a  mode-locked  ruby  laser  that  produced  pulses  of  5  -  10  picoseconds 
duration  and  approximately  56  watts  of  power,  transient  stimulated  Raman  scattering 
was  observed  in  all  liquids  that  were  tested  including  liquids  such  as  water,  carbon 
tetrachloride  and  methanol  in  which  stimulated  scatter’ ng  is  normally  difficult  to 
obtain.  Energy- conversions  as  high  as  20$  were  observed  together  with  multiple  order 
Stokes  generation  ar  beam  trapping.  Stimulated  scattering  was  also  observed  in  a 
variety  of  gases  including  N  ,  0 2,  CO^,  N^O,  SFg,  CH,  ,  C^H,  ,  C  Hg,  HC1  and  HBr.  Prior 
to  these  experiments,  stimulated  Raman  scattering  had  been  observed  only  in  H  , 
and  CH^.  Self  trapping  in  a  collimated  beam  was  observed  for  the  first  time  In 
gases,  and  evidence  of  a  strong  optical  Stark  shift  was  observed.  It  is  felt  that 
stimulated  Raman  scattering  could  have  important  application  in  the  determination 
of  relaxation  rates  of  materials  of  interest  for  chemical  lasers.  The  scattering 
provides  a  means  of  selectively  exciting  a  non-equilibrium  population  on  a  very 
short  time  scale.  The  subsequent  decay  can  then  be  monitored  to  determine  the 
relaxation  rates. 

The  work  on  stimulated  Raman  scattering  was  carried  out  in  collaboration  with 
N.  Bloembergen,  R.  L.  Carman,  F.  Shimizu  and  J.  Reintjes  of  Harvard  University.  All 
of  the  experimental  results  on  stimulated  scattering  reported  herein  were  obtained 
at  United  Aircraft  Research  Laboratories. 

The  continuing  effort  in  organic  dye  laser  technology  has  led  to  improved  pumping 
and  modulation  capabilities.  Dye  lasers  are  capable  of  producing  energetic  short 
pulses  at  a  wide  variety  of  wavelengths  and  should  have  significant  military  and 
industrial  as  well  as  basic  research  applications. 

Initial  investigations  of  the  feasibility  of  obtaining  stimulated  emission 
from  laser  produced  plasmas  have  been  carried  out.  Picosecond  pulse  excitation 
provides  an  extremely  fast  pumping  rate  and  might  be  used  to  populate  levels  having 
a  very  rapid  decay  time.  With  such  a  pumping  technique  it  might  be  possible  to 
obtain  stimulated  emission  at  ultraviolet  or  shorter  wavelengths.  Further  experi¬ 
ments  in  this  area  will  be  carried  out  upon  completion  of  a  corporation  owned  10  joule, 
.1  nanosecond  neodynium  laser  facility. 

The  extension  of  the  ideas  of  magnetic  adiabatic  passage  into  the  optical  region 
and  the  first  experimental  demonstration  of  the  optical  adiabatic  rapid  passage  were 
carried  out  under  the  present  contract  and  are  reported  in  detail  in  the  Third  Annual 
Report,  H920479-13,  March  31,  19^9 •  A  separate  research  program  on  optical  adiabatic 
passage  in  gases  has  been  established  and  is  being  funded  under  contract  with  the 
Army  Research  Office,  Durham,  North  Carolina  (Contract  DAHC04-70=C-00l4) .  Work  in 
this  area  under  the  present  contract  has,  therefore,  been  discontinued. 
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Section  2 

ANALYTICAL  DESCRIPTION  OF  ULTRA SHORT  OPTICAL  PULSE  PROPAGATION 

IN  A  RESONANT  MEDIUM 

2.1  Introduction 

Recent  advances  in  laser  technology  have  led  to  the  production  of  coherent 
optical  pulses  having  durations  in  the  picosecond  (10"^  sec  )  regime^'^.  Such  time 
intervals  are  comparable  to  or  shorter  than  the  relaxation  times  associated  with  the 
energy  levels  of  many  atomic  systems.  The  interaction  of  radiation  with  matter  on 
such  short  time  scales  gives  rise  to  phenomena  whicfy  as  a  result  of  quantum  mechanical 
coherence  effects,  cannot  be  described  by  the  rate  equation  analysis  developed 
previously  (3-5)  for  the  treatment  of  much  longer  pulses. 

The  novelty  of  the  effects  which  may  appear  in  the  short  time  resonant  response 
of  atomic  systems  has  been  brought  out  quite  strikingly  by  the  recent  discovery  of 
self- induced  transparency  (6,7,8).  jn  this  effect,  the  leading  edge  of  an  optical 
pulse  is  used  to  invert  an  atomic  population  while  the  trailing  edge  returns  the 
population  to  its  initial  state  by  means  of  stimulated  emission.  The  process  is 
realizable  if  it  takes  place  in  a  time  that  is  short  compared  to  the  incoherent 
damping  time  of  the  resonant  atomic  systems,  i.e.,  to  the  homogeneous  broadening 
time  of  the  medium.  When  all  conditions  for  this  process  are  met,  it  is  found  that 
a  steady  state  pulse  profile  is  established  and  that  this  pulse  envelope  then  propa¬ 
gates  through  the  medium  at  a  velocity  which  may  be  considerably  less  than  the  light 
velocity  in  the  medium.  What  is  perhaps  most  remarkable  is  that  even  those  atoms 
that  are  off  resonance  due  to  inhomogeneous  broadening  can  partake  of  this  process 
in  such  a  way  that  they  are  returned  to  their  initial  state.  Within  the  theoretical 
framework  that  has  been  used  to  describe  this  effect  it  has  been  shown  (?>®)  that 
such  steady  state  propagation  can  take  place  only  if  the  profile  of  the  electric 
field  is  of  a  special  form,  namely  that  of  a  hyperbolic  secant.  Many  salient  features 
of  this  effect  have  been  considered^)  since  its  discovery.  The  possibility  of 
analogous  effects  in  semiconductors  has  also  been  proposed(^3)  and  somewhat  similar 
effects  have  been  observed  in  the  study  of  neuristor  waveforms. 

In  addition  to  the  anomalous  transmission  property  of  ultra short  optical  pulses, 
the  amplification  of  such  pulses  has  also  drawn  considerable  attention.  A  number  of 
analytical  results  have  been  obtained  here  as  well(9,10,ll) .  One  expects  that  ampli¬ 
fication  processes  will  ultimately  be  limited  by  non-resonant  loss  mechanisms.  If 
these  are  introduced  in  a  phenomenological  way  through  a  conductivity,  then  the  ad 
hoc  assumption  that  there  is  a  steady  state  pulse  propagating  at  the  light  velocity 
may  be  verified  by  direct  computation. 


Whenever  it  becomes  necessary  to  extend  the  range  of  validity  of  a  theory  to 
encompass  new  phenomena,  it  is  useful  to  seek  limiting  cases  of  the  new  formal ism(®> 12, 
1  *  1  '  which  admit  of  exact  solutions  of  the  type  referred  to  above.  While  the 
experimentalist  is  rarely  moved  by  theoretical  descriptions  that  fail  to  provide 


for  all  facets  of  a  phenomenon  as  it  is  known  to  exist  in  the  experimental  situation, 
such  as  relaxation  times,  inhomogeneous  broadening,  etc.,  it  should  be  emphasized 
that  many  of  the  most  interesting  effects  in  ultra-short  pulse  propagation  appear 
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already  In  much  simpler  theoretical  contexts  and  an  understanding  of  even  these 
simpler  equations,  as  opposed  to  the  new  generation  of  numerical  results,  is  far  from 
complete.  The  most  promising  method  of  attack  is  clearly  through  a  simultaneous 
application  of  both  analytical  and  computational  techniques  and  the  field  of  opticnl 
pulse  propagation  provides  an  ideal  opportunity  for  the  application  of  this  "synergeti<V 
approach.  In  fact,  self  induced  transparency  was  discovered  from  an  analysis  of 
numerical  solutions  of  the  appropriate  equations. 


The  present  paper  summarizes  the  success  that  has  thus  far  been  achieved  in 
describing  the  novel  aspects  of  ultrashort  optical  pulse  propagation  by  analytical 
treatment.  In  addition  to  the  above-mentioned  results,  two  relatively  simple  models 
have  been  devised  which  describe  a  number  of  other  effects  that  have  been  observed 
both  experimentally^®**^)  and.  Qs  output  from  machine  computations  based  upon  more 
complete  theoretical  descriptions ' 12-1  '.  The  first  model^'^'  one  in  which  in¬ 
homogeneous  broadening  is  neglected.  The  physical  situation  most  closely  related  to 
such  a  model  is  that  of  propagation  under  conditions  of  extreme  saturation  broaden¬ 
ing'1®*1^.  Although  the  problem  under  consideration  involves  a  coupling  between 
radiation  and  matter  that  is  too  strong  to  be  treated  by  perturbation  theory,  a  fairly 
extensive  analytical  treatment  of  this  model  is  still  possible  since  it  expresses  this 
interaction  in  terms  of  a  single  nonlinear  partial  differential  equation  which  arose 
long  ago  in  differential  geometry.  The  techniques  developed  about  the  turn  of  the 
century  for  obtaining  solutions  to  this  equation  may  be  employed  to  great  advantage. 

(20-23) 

Certain  other  phenomena,  notably  that  of  photon  echo  ,  require  for  their 

explanation  the  relative  dephasing  of  atoms  that  results  when  inhomogeneous  broadening 
is  present.  This  effect  is  also  an  example  of  the  collective  superradiant  state'*  ' 
in  which  energy  is  radiated  coherently  into  the  electromagnetic  field.  Here  again,  it 
is  possible  to  construct  a  soluble  model^2^)  in  which  the  reaction  of  stimulated 
emission  back  on  the  incident  wave  is  taken  into  account.  If  one  is  willing  to  forgo 
consideration  of  the  detailed  structure  of  pulse  shapes,  the  time  dependence  of  the 
pulse  may  be  assumed  to  be  that  of  a  delta  function  and  interest  confined  to  the 
spatially  dependent  amplitude  of  such  delta  function  pulses.  Only  the  time  integral 
of  such  a  pulse  shape  is  meaningful,  of  course,  but  such  a  time  integral  has  been 
shown  to  be  precisely  the  quantity  of  interest  in  the  treatment  of  ultrashort  pulses. 
The  area  theorem^?*®' ,  which  is  so  useful  in  understanding  short  pulse  phenomena,  is 
also  found  to  govern  the  spatial  evolution  of  the  amplitude  functions  introduced  in 
this  model. 


Although  much  of  the  physical  insight  required  for  an  understanding  of  these 
propagation  effects  may  be  obtained  from  a  consideration  of  the  interaction  of  light 
with  a  system  of  two-level  atoms,  it  should  be  emphasized  that  the  results  thus 
obtained  may  require  modification  when  level  degeneracy  is  included^®**-®) . 

2.2  Basic  Equations 

We  begin  by  summarizing  the  standard  semiclassical  description  of  the  interaction 
of  an  electromagnetic  wave  with  an  assembly  of  two-]evel  systems.  The  optical  field 
in  the  form  of  a  plane  polarized  electromagnetic  pulse  may  be  characterized  by  its 
electric  field  vector  E(_r,t)  which  satisfies  the  usual  wave  equation 
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£E  I  d*£ 


4w 

“c5"  at2 


(2.1) 


where  o  is  a  conductivity  that  ic  introduced  to  cimulatc  nonresonant  losses  in  the 
medium,  c  is  the  velocity  of  light  in  the  medium  and  P  in  the  polarization  of  the 
medium  that  ic  induced  by  the  electromagnetic  wave.  For  a  medium  consisting  of  an 
assemblage  cf  noninteracting  two-level  systems  distributed  with  a  uniform  density  n0, 
this  polarization  is  noP  where  p  ic  the  polarization  of  an  individual  two-level 
system. 


The  polarization  of  an  individual  system  may  in  turn  be  obtained  from  its  micro¬ 
scopic  description  by  the  usual  prescription 


p=Tr(*P) 


(2.2) 


where  p  is  the  density  matrix  of  the  two-level  syctcm(^5)#  and  £  is  the  polarization 
operator.  The  time  dependence  of  p  is  given  by  the  quantum  mechanical  Liouvillc 
theorem 


(2.3) 


where  X  is  the  total  Hamiltonian  of  an  individual  two-level  system.  The  time  depen¬ 
dence  of  an  arbitrary  operate-,  0  is  governed  by  the  relation 


i  h 


(2.4) 


For  later  use,  it  proves  convenient  to  recognize  that  operators  not  containing 
explicit  time  dependence  also  satisfy 


(2,5) 

The  Hamiltonian  of  a  two  level  system  interacting  with  a  classical  electromagnetic 
field  may  be  adequately  represented  by 

»  ■  *0  ♦  V  (2.6) 

where  X  is  the  Hamiltonian  of  the  isolated  two- level  system  and 
o 

v  .  -g  p  (2.7) 

is  the  interaction  energy  in  dipole  approximation. 

The  wave  function  for  the  isolated  two-level  system  may  be  written 

*  o0lt)u0(£)  +  Ob(t)ubl£,)  (2.8) 
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where  and  ore  eigenfunctions  of  the  system  and  satisfy 

#QUa  '  couo  ,0  ‘  0,b 

Such  *  rgy  eigenfunctions  provide  the  specific  representations 


(2.9) 


(2.10) 


where 


P«/3  * 


l  Ufi 


(2.11) 

po$  ’  °o#  °0 

and  the  levels  are  labeled  such  that  Er  >  E^.  The  vanishing  of  the  diagonal  elements 
in  P  signifies  the  assumed  absence  of  any  permanent  dipole  moment  in  the  system  under 
consideration. 


In  addition  to  the  polarization  P,  the  difference  in  population  between  upper  and 
lower  states  n,  is  also  of  interest  and  may  be  expressed  in  the  form 

n  *  n0  IPqo  *  ^bbJ  *  n0  Trip*,)  (2.12) 


where  e_  is  the  Pauli  spin  matrix 


V 


(2.13) 


The  tine  dependence  of  p  can  be  conveniently  obtained  by  taking  the  trace  of  the 
operator  equation 


2  •• 

*  p 


lie-' 


(2.14) 


which  follows  from  Eq.  (2.4).  With  the  representations  given  above, 
side  of  this  eouation  vanishes  while 


e.*|  *  [p,»0j  ■ 


the  right-hand 


(2.15a) 


and 


l’\ 


I  fobfbo 


t>o  ®  \ 

-bO'Ob' 


(2.15b) 


(2.15c) 
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where 


(2.16) 


Equation  (2.12)  therefore  reduces  to 

e  +  f  ■ 


fi 


£  ‘  Pob  Pbo 


(2.17) 


Application  of  the  trace  operation  converts  this  operator  equation  into 


(2.18) 


where  it  has  been  assumed  that  =  P^&  =  ft* .  The  factor  of  l/3  in  parenthesis  is 
to  be  included  if  all  possible  spatial  orientations  of  the  two  level  systems  are 
permitted(26,7 5  ). 


From  a  direct  multiplication  of  the  quantities  involved  there  follows 

*wob 


tfiP  =  [p,/ 


m  ■■  - 


Pi » 6 


(2.19) 


On  the  other  hand,  multiplication  of  Eq.  (2.3)  by  n  3  and  application  of  the  trace 
operation  yields 


•rhri^0E-Tr|/>[P.crz]| 


(2.20) 


Calculation  of  the  trace  of  Eq.  (2.19)  then  finally  leads  to  the  equality 

27L  „  ' 

91 s  -r— 2-  E-  P  .  , 

tuuob  (2.21) 

Equations  (2.l8)  and  (2.2l)  provide  a  convenient  starting  point  for  describing  the 
response  of  a  two-level  system  to  an  external  electromagnetic  field. 

Having  formulated  the  response  of  the  two-level  systems  to  the  incident  optical 
field,  we  now  turn  to  a  consideration  of  the  reaction  of  the  medium  back  upon  the 
incident  wave.  It  is  the  self-consistent  evolution  of  these  two  processes  that  is 
the  essence  of  the  problem  under  consideration.  On  the  right-hand  side  of  Eq.  (2.1), 
the  term  d2p/dt2  may  be  approximated  by  -ui|p  where  u>0  is  the  carrier  frequency  of 
the  incident  pulse.  This  follows  immediately  if  is  replaced  by  u)0  in  Eq.  (2.18) 

and  the  term  of  the  right-hand  side  of  that  equation  is  neglected.  Neglect  of  this 
term  is  equivalent  to  a  neglect  of  the  backscattered  wave  that  is  produced  as  the  in¬ 
cident  wave  traverses  the  medium.  From  Eqs.  (2.1)  and  (2.18),  one  sees  that  this 
neglect  is  permissible  provided  that  Tim  «  1  which  will  be  nearly  always 

satisfied. 
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Since  even  the  shortest  pulses  produced  to  date  contain  many  optical  cycles, 
it  is  appropriate  to  write  the  electric  field  in  terms  of  a  carrier  wave  as  well  as 
envelope  and  phase  functions  g,(r,t)  and  $(r,t)  respectively  which  vary  slowly  on 
the  length  and  time  scales  of  the  carrier  wave.  Hence  we  write 


E(r,t)  =  EJ?(r,t)cos 


(2.22) 


and  assume  u>0  8  »  9  E  /dt,  kQF  »  9  t  /9x  plus  similar  inequalities  for  $.  A  character¬ 
istic  field  magnitude  EQ  has  also  been  extracted  from  the  amplitude.  Because  of  the 
assumption  that  £(r,t)  and  $(r,t)  vary  slowly  compared  to  the  carrier  wave,  Eq.  (2.1) 
may  be  reduced  to  a  much  simpler  form.  In  particular,  only  first  derivatives  of 
g  (r,t)  and  $(r,t)  need  be  retained  on  the  left-hand  side  of  this  equation  when  it 
is  expressed  in  terms  of  F  and  $.  The  solution  thus  obtained  is  customarily  referred 

i 

to  as  the  solution  in  the  slowly  varying  envelope  approximation. 


In  general,  it  is  appropriate  to  consider  not  a  single  transition  frequency  u>ab 
but  a  spread  in  transition  frequencies  about  uub,  so  called  inhomogeneous  broadening. 

It  is  convenient  to  analyze  the  situation  in  which  this  distribution  is  symmetric  about 
cuab  and  the  carrier  frequency  of  the  incident  optical  pulse  is  at  this  frequency, 
i*e*>  U)Q  =  (%b. 


Equation(2.l),  specialized  to  a  plane  wave  traveling  in  a  positive  x  direction 
now  becomes 


where 


^  +  2  TT-o-^sm  $(x,t)4^C0S  <5Hx,t)  =  27r^Q  JdAcug(Auj)p(Acj,x,t) 
<*>(x,t)=k0x-w0t  +  <Mx,t)  .  df=“^~+c 


(2.23) 

(2.23a) 


and  Aid  =  oi  -  u)0.  Although  the  idealization  of  an  infinite  plane  wave  front  is 
convenient  for  theoretical  purposes,  it  should  be  emphasized  that  transverse  mode 
structure  may  be  important  in  experimental  situations. 

The  spectrum  g(A'u)  which  characterizes  the  inhomogeneous  broadening  is  assumed 
to  be  normalized  so  that 


/•CO 

J  dAuiglAui)  = 

-00 


(2.24a) 


The  causal  Green's  function  for  Eq.  (2.l8),  i.e.,  the  solution  of 

+  uJ0bZG  =  -S(t-t') 
dt2 

which  satisfies  G  =  0  for  t  <  t'  is 


(2.24b) 


G(t|t‘)=  -  V  (t— t') sin ajQb(t— 1‘) 


(2.25) 
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in  which  u(t)  is  the  unit  step  function.  Multiplication  of  Eqs.  (2.l8)  and  (2.25)  by 
G(tjt')  and  p  respectively,  subtracting  and  integrating  over  all  time  yields 


p(Acu(x,t)  =-2p£1  fdt'sin cu(t— t')A  (Aui,x,t')£(x,t')cos  <£(x,t') 

J-  oo 


(2.26) 


where  A  =  n/nQ,  and  Q,  the  Rabi  frequency,  is  given  by  Q  =  EQ  m-  The  factor  of 
1/3  due  to  orientational  averaging  is  neglected.  This  expression  for  the  polarization 
may  now  be  decomposed  into  parts  which  are  in  phase  and  tt/2  out  of  phase  with  the 
electric  field.  One  finds  that 


where 


p  =  P  [P  (Aw.x.t)  sm($(xt))+  2  (Aaj,x,t)cos($(x,t))] 


(?  =  (x,t')  A(Acu,x,t')cos  jAuj(t-t')  +<£(x,t)-<£(x,t')j 

2  -  -Zlj  dt'^lx.t'jA  ( Acu.x.t’)  sin  [a  u>  (t— t')  +  <£  (x,t)-  <Mx,t')J 


(2.27) 


(2.28a) 


(2.28b) 


In  obtaining  this  result,  terms  near  the  second  harmonic  of  have  been  discarded. 
However,  it  should  be  emphasized  that  in  obtaining  these  results  there  has  been  no 
assumption  that  A  ,  2  ,  and  A  vary  slowly  compared  to  the  carrier  wave. 

Equation  (2.23)  may  now  be  decomposed  into  the  pair  of  relations 


where 


|jy+27rcr£  =  (ca7^)J*dAujg(Acu)/:?(Aa;,x,t) 

(2.29a) 

Z  =  -(ca'/Sl) >fdAug(Aj2(Aco,x,1) 

(2.29b) 

2 

a1  =  2vn0uJ0P  /he 

(2.30) 

The  functions  P  and  2  are  readily  shown  to  satisfy  the  differential  equations 

dP  „  /  \ 

Wrn£A  +  ra+7TF  (2.31a) 


ff  +  (*«♦#)* 
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When  time  dependence  near  the  second  harmonic  of  u)Qis  also  neglected  in  Eq.  (2.21) 
it  is  equivalent  to 


M.  --  azi  0-  — 

at  -  U^V'  at  /  (2.32) 

When  the  slowly  varying  envelope  approximation  is  used,  the  second  term  on  the  right- 
hand  side  may  be  discarded.  It  is  convenient  to  introduce  the  transformation 


t  =  ft  (t-x/c),£  =  (&« /c)x  (2.33) 

The  constant  e  is  essentially  the  ratio  of  energy  stored  in  the  medium  to  the  energy 
of  the  wave,  i.e. 


*  5"o1W<Eo2''2’r> 


(2.34) 


Neglecting  the  conductivity  at  this  point,  Eqs.  (2.29)  are  transformed  to 

dZ  _  rco 


Ar  =  f  dAwg(Aw)/9(Aw,C,r) 
•'-(JO 

=  j[®  dAwg(M2(Autf,r) 


(2.35) 


(2.36) 

while  Eqs.  (2.31)  and  (2.32),  in  the  slowly  varying  envelope  approximation,  become 


0-  -cp 

(2.37a) 

(2. 37h) 

a? 

(2.37c) 

where  f  =  Au)/ft.  Equations  (2.37)  describe  how  the  field  amplitude  f  and  phase  $ 
determine  P  ,  2  and  A  for  a  two-level  system  that  is  off  resonance  by  an  amount 
Ao).  Equations  (2.35)  and  (2.36)  show  how  the  polarization  due  to  a  distribution  of 
such  systems  reacts  back  on  the  amplitude  and  phase. 

Equations  (2.37)  also  arise  in  nuclear  magnetic  resonance  studies  in  which  an 
oscillating  magnetic  field  interacts  with  an  assemblage  of  two  level  systems  which 
possess  a  magnetic  moment.  Such  studies  have  been  confined  to  samples  that  are  of  a 
sufficiently  small  size  that  the  reaction  of  the  induced  field  back  on  the  exciting 
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field  could  be  ignored.  Equations  (2.37)  may  then  be  solved  for  a  specified  external 
field.  This  is  not  the  case  in  the  situation  envisioned  here.  A  satisfactory 
description  of  optical  pulse  propagation  is  only  obtained  when  Eqs.  (2.36)  and 
(2.37)  are  solved  self  consistently. 

At  this  point,  it  may  be  noted  that  Eqs.  (2.35)  and  (2.37)  contain  two  conserva¬ 
tion  laws.  Conservation  of  energy  follows  upon  multiplication  of  Eq.  (2.35)  by  8 
and  introduction  of  Eq.  (2.37a)  which  yields 

2 

T  r  00  dAwg(Aw)A  =  0  (2.38 

2  oT  J-cd 

Multiplication  of  Eqs.  (2.37)  by  A  , P  and  2  respectively  along  with  summation  of 
the  resulting  equations  yields  an  exact  differential  which  is  equivalent  to 

PZ+22=  I  (2-39) 

In  this  result,  a  constant  of  integration  has  been  set  equal  to  unity  since  in  the 
usual  applications  of  the  theory  one  has  P  (£,  -a>  )  =  2  (e,-~)  =  0,  A  (c,  -»)  =  ±1. 

The  form  of  Eq.  (2.39)  enables  one  to  interpret  the  response  of  a  two  level  system 
in  terms  of  the  motion  of  a  vector  on  the  surface  of  a  sphere  in  a  P  ,2  ,  A  space. 


2.3  Self-Induced  Transparency  and  the  Area  Theorem 

Up  to  the  present  time,  the  full  set  of  equations  given  by  Eq.  (2.35  -  2.37) 
has  received  little  attention.  However,  if  one  adopts  the  consistent  set  of  assumptions 
that  the  phase  term  cp  is  initially  zero,  that  the  carrier  frequency  is  at  the  center 
of  a  symmetrically  broadened  line  (i.e.,  g(Aiu)  =  g(-Atu),  and  that  2  is  an  odd  function 
of  Au5.  Then  from  Eq.  (2.36)  one  sees  that  the  source  term  governing  variations  in  cp 
is  zero  so  that  will  remain  zero.  This  form  of  the  theory,  particularly  with  the 
aid  of  numerical  computations,^)  has  provided  considerable  insight  into  the 
subject  of  ultrashort  pulse  propagation. 

Even  this  specialized  form  of  the  basic  equations  has  yielded  only  steady  state 
solutions.  These  include  both  the  solitary  wave  solution  of  self-induced  trans¬ 
parency^^)  and  infinite  wave  train  solutions > 28,29)  contain  the  solitary 

wave  as  a  limiting  case.  Only  the  former  will  be  discussed  here;  infinite  wave  train 
solutions  will  be  discussed  later  in  connection  with  a  somewhat  more  specialized 
theoretical  model. 

For  a  steady  state  solution,  one  may  assume  that  £  ,  P  ,  2  ,  and  A  are  functions 
of  a  single  dimensionless  variable  w  =  (t-x/V)/xp  where  V  is  the  velocity  of  the 
pulse  and  xp  is  a  parameter  having  the  dimensions  of  time.  It  will  be  shown  that 
x  may  be  directly  related  to  the  pulse  width. 

Equations  (2.35)  and  (2.37a)  may  now  be  combined  and  integrated  to  yield 
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—  £  (w)  +  €  il  J  ^  dAojfl  (Acu)j2  (Acu  ,w)  /Aoi=  0 


(2.40) 


An  integration  constant  has  been  set  equal  to  zero  in  this  result  since  £  and  2 
are  zero  before  arrival  of  the  pulse.  If  this  equation  is  divided  by£(w)  and 
differentiated  with  respect  to  w,  one  obtains 


_d_ 

dw 


=  0 


(2.41) 


Now  the  function  g(Aiu)  contains  a  parameter  such  as  T*  that  determines  the  width 
of  the  inhomogeneous  broadening,  e.g. 


g(A“)  =  ^rexp[-(AwT2*/2)2]  (2.42) 

Although  Tp  and  V,  and  hence  id  depend  in  an  implicit  way  upon  Tp,  there  is  no  explicit 
dependence  of  £  upon  T*.  Also,  since  2  is  the  response  of  an  individual  two-level 
system,  it  is  also  independent  of  T2  except  for  the  implicit  dependence  contained  in 
w.  Consequently,  the  function  in  parentheses,  the  integral  of  Eq.  (2.4l),  does  not 
contain  explicit  dependence  upon  TJ.  Now  if  it  is  assumed  that  the  theory  must  be 
valid  for  arbitrarily  large  values  of  Tj),  then  one  may  invoke  Lerch's  theorem^0) 
to  justify  the  conclusion  that  the  term  multiplying  g(Au>)  in  the  integrand  must 
itself  be  equal  to  zero.  It  then  follows  that 


.2(A cu,w)  =  X  (Aw)£(w) 


(2.43) 


where  X  (a<d)  is  an  as  yet  unknown  function  of  the  detuning.  V/hen  this  result  is 
introduced  into  Eq.  (2.40),  one  finds  that  the  velocity  of  the  envelope  function 
is  given  by 


1 

V 


^1  +  €  si  j  dAtug  (Aw  )x(Aw)/Ao;j 


(2.44) 


To  obtain  the  form  of  the  envelope  function  £,  it  is  first  noted  that  Eqs.  (2.37^, 
c),  with  the  phase  term  cp  set  equal  to  zero,  are  equivalent  to  the  linear  equation 


||-  +  f2 

where  \  =  P  +  i2 •  Setting 

<j>  (w)  =  ^rp  f  dw'^(w') 

•'-GO 


(2.45) 


(2.46) 
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and.  introducing  Eq.  (2.43),  one  finds  that  the  solution  of  Eq.  (2.45)  which  reduces 
to  P{  -a>  )  =  0,  A  (-=>)  =  -1  is 


P  =  -(i  -f  x)  sin  4> 


(2.47) 


ft  -  “fx-^“^x)cos^ 


Substitution  of  Eq.  (2.43),  (2.47)  and  (2.48)  into  Eq.  (2.39)  yields 


/a  \  flrpAqj 

*  ^  l+ (Acdtd)2 


Also,  from  Eq.  (2.38),  the  energy  conservation  lav, 


(2.48) 


(2.49) 


=  2  y2  ( I  -  cos  <£)  =  4  yl  sin2  (4>/2.) 

,dw  / 


where 


(2.50) 


yz=er(nTpf/ (■%•-!) 


(2.51) 


and 


r  =  1  -jj dAu;g(Acij)  Aoix^Aui) 


The  solution  of  Eq.  (2.50)  that  vanishes  as  w-»±®  is 

-i  y* 

=  4  tan  e 


(2.52) 


(2.53) 


Hence  the  electric  field  envelope  is 


—  i  C  -  —  =  — sech  y  w 

h  rp  d  w  r p 


(2.54) 


and  one  sees  from  the  definition  of  w  that  Tp  will  play  the  role  of  pulse  half  width 
provided  y  is  set  equal  to  unity. 

If  the  maximum  pulse  height  is  equated  to  2EQ,  then  Eq.  (2.54)  yields 


^Tp  = 


(2.55) 


which  determines  the  pulse  width  in  terms  of  its  amplitude  and  ^  which  characterizes 
the  medium. 
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From  Eqs.  (2.43),  (2.47), 
given  two-level  system  is' 3) 


(2.48)  and  (2.49)  one  finds  that  the  response  of  a 


where 


yy  =  -i+2Dsin2(<£/2) 

P  -  -  Dsin<£ 

2=  Z  AuiTp Dsin(<^2) 

(2.56) 

D=  [l+(Au»rp)2] 

(2.57a) 

Also,  it  is  now  possible  to  return  to  Eq.  (2.44)  for  the  envelope  velocity  and 
consider  its  dependence  upon  T2*  Combining  Eqs.  (2.44)  and  (2.49),  one  finds  that 
V  is  given  by 


J_ 

C 


a'  rD2  /  dAu»g(Aai)D 

y  J  —  00 


(2.57b) 


V=(l~r\v0"c 


(2.58) 


Here  Vo  is  the  velocity  in  the  limit  of  no  inhomogeneous  broadening.  For  g(Auu) 
given  by  Eq.  (2.42),  one  obtains 


j jAwg(Aai)0  =  2kexp(K2)E*  fc(«) 


(2.59) 


where 

K=T2/2TP  (2.60) 

and  Erfc(k)  is  the  compliment  of  the  error  functional).  For  T2  »  t^,  F  -*  1  and 
■X*  T  P 

V- V0.  As  T2  becomes  much  less  than  Tp,  a  much  smaller  percentage  of  the  atoms  are 

on  resonance.  One  then  finds  I  -*•  0  and  hence  V  -»  c. 

Finally,  it  has  been  noted(^2) 

that  if  the  carrier  frequency  is  not  confined  to 
the  center  of  a  symmetric  line  but  is  in  fact  far-off  resonance,  the  expression  for 
the  velocity  given  in  Eq.  (2.44)  goes  over  to  the  usual  result  for  the  velocity  of 
a  wave  in  a  dispersive  medium. 


The  above  results,  along  with  the  infinite  wave  train  solutions  which  correspond 
to  libration  and  oscillation  solutions  of  Eq.  (2.50)  are  the  only  analytical  solutions 
of  the  inhomogeneously  broadened  version  of  Eqs.  (2.35  -  2.37)  that  have  been  reported 
to  date.  However,  further  analytical  progress  is  still  possible  if  one  confines 
attention  to  the  area  under  the  envelope  curved®).  Defining 


0(x)=  ft 


(2.61) 
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the  equation  governing  the  variation  in  9  is  readily  obtained  by  integrating  Eq.  (2.29a) 
over  all  time.  Using  Eqs.  (2.28b)  and  (2.31b)  with  <p  again  set  equal  to  zero  one  finds 


66 

dx 


(2.62) 


where  k  =  2tt  o/c.  Since  the  indicated  limiting  process  is  somewhat  delicate,  a 
method  of  carrying  it  out  is  now  given  in  some  detail.  First,  one  may  introduce  the 
well  known  forms 


Lim  sinAcut 
t— oo  Acu 


=  7t8(Acl») 


Lim  cosAcut  _/  j_  _  _p_ ' 
f— Acu  lAcu  Acu, 


(2.63) 


where  P  denotes  a  principal  value.  Using  another  standard  representation  for  the 
delta  function  and  the  principal  value, 

Lim  COSAcut  _  Lirn  r  J _ Acu  ~| 

t-*'00  Acu  €—0  L  Acu  c2+(Acl)2J 


=  7T 


5  (Acu)  Lim  e 
cu  e— 0 


0 


When  these  results  are  used  in  Eq.  (2.62),  one  obtains 


66 

dx 


+  K  6 


2 


J  dt'£(x,t')/V(o,x,t‘) 


(2.64) 


(2.65) 


where 


a  =  2  7rg(o)a' 


(2.66) 


and  o1'  is  as  defined  in  Eq.  (2.30).  As  will  be  shown  later,  the  population  of  on- 
resonance  atoms,  i.e.,  those  represented  by  A(0,x,t),  may  be  expressed  as 


A(o,x,t)=  ±cos  Jdt'^x.Oj  (2.67) 

where  the  upper  sign  is  to  be  used  if  the  population  is  initially  inverted  while 
the  lower  sign  is  used  if  the  population  is  initially  in  the  lower  level.  Hence 
Eq.  (2.65)  finally  takes  the  form 


:p  +  K0  =  +  f-sin0 

OX  d 


(2.68) 


Equation  (2.68)  is  the  area  theorem^®),  it  contains  the  key  to  an  understanding 
of  many  of  the  effects  which  occur  in  the  propagation  of  ultrashort  optical  pulses. 
Again, if  orientational  averaging  is  included,  the  factor  a  should  be  replaced  by  0/3. 
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For  o  =  0,  the  solution  of  Eq.  (2.68)  which  satisfies  9  =  0O  at  x  =  xQ  is 

ton  -f-  Mon  -^exp[±-2-(x-x0)]  (2.69) 

and  is  depicted  schematically  in  Fig.  1.  Modifications  of  this  result  due  to 
nonvanishing  conductivity  have  been  inferred'^®'  from  numerical  solutions  of  Eq.  (2.68). 

Since  Eq.  (2.68)  contains  a  choice  of  signs,  it  is  actually  two  distinct  differen¬ 
tial  equations.  The  two  solutions  are  obtained  from  Fig.  1  by  reading  the  diagram 
from  right  to  left  for  the  plus  sign  (amplifier)  and  from  left  to  right  for  the  minus 
sign  (attenuator) .  Hence  one  sees  that  an  infinitesimal  area  will  grow  to  tt  in  an 
amplifier  while  any  area  less  than  tt  will  evolve  to  zero  in  an  attenuator.  This 
second  result  allows  for  not  only  the  well-known  decay  of  a  pulse  as  it  propagates 
in  an  attenuator,  but  also  for  evolution  into  a  nonvanishing  zero  tt  pulse,  i.e., 
one  in  which  the  total  area  under  the  pulse  envelope  is  zero  but  the  area  under  the 
pulse  energy  ( — ^ )  is  not  zero.  This  is  of  course  possible  if  the  positive  portions 
of  a  pulse  envelope  are  equal  in  area  to  the  negative  portions.  Physically,  the 
regions  of  positive  and  negative  envelope  are  merely  regions  in  which  there  is  a 
relative  difference  of  180  degrees  in  the  phase  of  the  carrier  wave.  In  an 
attenuator,  initial  pulse  areas  between  tt  and  3tt  will  evolve  into  the  steady  state 
2tt  pulse  of  self- induced  transparency.  One  also  sees  that  the  2tt  pulse  is  unstable 
in  an  amplifier  and  will  evolve  into  either  a  tt  or  3tt  pulse. 

Figure  1  refers  only  to  the  total  area  of  a  pulse  and  gives  no  information  at 
all  about  the  possible  breakup  of  a  pulse  into  two  or  r  ir e  pulses  with  the  same 
total  area  or  of  whether  a  continually  amplifying  pulse  will  retain  an  area  of  n  by 
virtue  of  pulse  narrowing  or  by  developing  negative  regions  in  the  pulse  envelope. 

2.4  Steady  State  Pulse  in  an  Amplifier 

In  addition  to  the  self-induced  transparency  solution  in  an  attenuator,  a  some¬ 
what  similar  steady  state  result  may  be  obtained  in  an  amplifer  if  the  loss  term  0 
is  retained  in  Eq.  (2.29a).  This  was  first  recognized  by  observation  of  machine 
computations (13)  and  subsequently  obtained  analytically^).  Both  results  have  been 
obtained  in  the  limit  of  no  inhomogeneous  broadening  although  certain  cases  in  which 
homogeneous  broadening  is  retained  have  also  been  treated'^'.  From  Eqs.  (2.29),  (2. 31) 
and  (2.32)  the  relevant  equations  are 


(2.71) 


as  well  as  Eqs.  (2.37)*  The  ad  hoc  assumption  which  renders  the  analysis  tractable 
is  that  both  £  and  <p  travel  at  the  velocity  c.  The  differential  operators  in 
Eqs.  (2.70)  and  (2.71)  then  vanish  identically  and  the  problem  is  greatly  simplified. 
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It  should  be  emphasized  that  there  is  no  rigorous  Justification  for  this  assumption. 
However,  it  has  recently  been  notedU^)  that  if  a  steady  state  pul3e  at  any  velocity 
vjc  is  assumed,  then  the  resulting  numerical  solution  is  unstable  and  evolves 
asymptotically  into  a  pulse  propagating  at  v  e  c. 

From  Eq.  (2.71)  one  then  sees  that  2=0  and  hence  from  Eq.  (2.37b)  that 


(2.72) 


Consequently, 


(2.73) 

and  from  Eq.  (2.22)  one  sees  that  the  frequency  of  the  steady  state  pulse  is  always 
equal  to  uj  the  transition  frequency  of  the  two-level  svstem(9).  Since  2  =  0, 

Eq.  (2.37a)  follows  from  Eq.  (2.32)  without  the  use  of  the  slowly  varying  envelope 
approximation.  The  remaining  equations  are  now 


U'P 

(2.74) 

p--SiZH 

(2.75) 

a  --step 

(2.76) 

where  the  dot  signifies  differentiation  with  respect  to  t  -  x/c  and 


C=o-  Eo^Vow0!f> 


(2.77) 


Since  Eq.  (2.75)  and  (2.7 6)  have  the  parametric  representation 


with 


P  -  sin  0 

(2.78a) 

A  =  cos 

(2.78b) 

I* 

(2.79) 


Eq.  (2.74)  is  equivalent  to  the  differential  equation 
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-  sin  4> 

The  solution  is  immediate  and  leads  to 


(2.80) 


Z  =  sech[  ft(t-x/s)/c]  (2.81) 

Also,  the  population  is  seen  to  be  inverted  by  the  pulse  since  Eq.  (2.78b)  is  now 
equivalent  to 


A  =  —  tanh  [51(t— x/s)/C] 

Since 


(2.82) 


(2.83) 


and  also  because  the  vector  whose  components  are  P  and  2  is  rotated  through  an 
angle  n  during  the  passage  of  the  pulse,  the  result  given  in  Eq.  (2.8l)  is  customarily 
referred  to  as  a  n  pulse. 


Steady  3tato  pulse  propagation  in  an  amplifier  has  also  been  analyzed  without 
the  assumption  of  a  slowly  varying  envelope  and  phase^1^).  The  assumption  of  propa¬ 
gation  at  the  light  velocity  of  the  medium  is  retained,  however,  and  provides  the 
simplification  that  is  sufficient  to  permit  an  exact  solution.  For  pulses  that  arc 
many  optical  cycles  in  duration,  there  is  very  little  difference  between  the  pulse  shape 
obtained  with  this  more  exact  treatment  and  the  method  described  above,  as  is  to  be 
expected.  What  is  of  great  interest,  however,  is  the  prediction  of  phase  variation 
in  the  carrier  wave.  The  "chirp"  predicted  by  the  theory  is  proportional  to  the 
square  of  the  ratio  of  optical  period  to  pulse  width.  Such  a  result  could  not  be 
obtained  in  the  slowly  varying  phase  and  envelope  approximation  since  it  is  equivalent 
to  an  expansion  to  only  first  order  in  this  ratio. 


,(H) 


The  method  has  subsequently  been  extended' ■L'L/  to  include  the  effect  of  dispersion 
in  the  host  medium.  In  the  limit  of  large  dispersion  it  was  found  that  a  monotonic 
frequency  sweep  is  predicted.  Such  chirping  of  ultrashort  pulses  has  been  observed 
experimentally  and  offers  new  opportunities  for  pulse  compress  ion'  ->->>->  '  and  population 
inversion'  .  The  chirp  that  devclopos  in  the  presence  of  large  dispersion  is  found 
to  be  proportional  to  the  first  power  of  the  ratio  of  optical  period  to  pulse  width. 
Tliis  suggests  that  such  a  result  can  be  obtained  with  the  framework  of  the  slowly 
varying  phase  and  envelope  approximation  and  in  the  following  development,  the 
problem  will  be  approached  from  this  point  of  view. 


When  no  dispersion  is  present  and  the  pulse  is  assumed  to  propagate  at  the  velocity 
of  light  in  the  hoot  medium,  Eq.  (2.21)  reduces  to 
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r  -  <*P 


(2.84) 


where  w  =  t  -  x/c.  When  E  and  P  are  written  in  the  form  of  Eq.  (2.22)  and  (2.27) 
respectively  but  without  any  assumption  that  £  and  cp  are  slowly  varying,  Eq.  (2.80) 
yields  the  pair  of  relations 


■1  (pd±  +  d2) 

\  a  w  +  a  w ) 

(2.85) 

(w  o“aw)= 3 

(2.86) 

Eqs.  (2.37),  with  Ao)  again  set  equal  to  zero,  are  also  applicable, 
of  Eqs.  (2.85)  with  (2.37c)  and  Eq.  (2.86)  with  (2.37b)  leads  to 

U--P 


u02 z  -  Si£  ft 

Eqs.  (2.37b)  and  (2.37c)  may  be  combined  to  yield 


Combination 

(2.87a) 

(2.87b) 


(2.88) 


in  which  Eq.  (2.39)  has  been  employed.  From  Eqs.  (2.87),  (2.37)  and  ( 2.39 )>  one 
obtains 


2Z-  QZft%~ftZ)/  [\+QZft  2  ]  |^-  = 


-to- 


a  [1  -A2] 
i  +a2ft  2 


(2.89) 


Finally,  Eq.  (2.88)  yields 


d4>  _  2 

aw  "woa 


a  =  U/C  tL>c 


ft‘ 


i+(a/V)2 


\-ftZ  ) 

[i+(aA)2]j 


(2.90) 


(2.91) 


Unfortunately,  one  cannot  solve  Eq.  (2.89)  explicitly  for  the  explicit  time 
dependence  of  ft  .  Solution  of  Eq.  (2.89)  merely  leads  to  the  inexplicit  relation 
equation  for  ft  which  has  the  solution 


aw0w  =  ft  -  (a2+i )  tan  h"'/V 


(2.92) 
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For  a  »  1  this  reduces  to  Eq.  (2.82)  as  is  to  be  expected. 
The  "instantaneous"  frequency  is 


d  4> 

“inst 8  “o  - 

For  A  =  ±1,  i.e.,  at  both  ends  of  the  pulse 


(2.93) 


“w="°[|-aZ/(l+oZ)] 


(2.94) 


while  at  the  center  of  the  pulse  A  =  0  and 


winsts%(|+<>2) 


(2.95) 


The  fractional  shift  in  frequency  is  2a2/(l  +  a2).  For  a  picosecond  pulse  at 
lp,,  a  ~  10_3  and  the  fractional  frequency  shift  is  ~10  .  The  absolute  frequency  shift 
is  300  mHz. 


The  calculation  outlined  above  was  subsequently^  ^extended  to  include  dispersion 
in  the  host  medium.  The  dispersion  was  treated  by  standard  methods  of  linear  wave 
propagation.  It  was  found  that  in  the  limit  of  large  dispersion,  the  chirp  becomes 
proportional  to  the  first  power  of  the  ratio  of  optical  period  to  pulse  duration 
rather  than  to  the  second  power  as  was  found  in  the  previous  calculation.  As  is  to 
be  expected,  then,  this  limiting  case  can  be  treated  in  the  slowly  varying  envelope 
and  phase  approximation  and  this  formulation  is  developed  below. 


The  wave  equation  given  in  Eq.  (2.1)  is  readily  modified  to  include  effects 
arising  from  the  presence  of  a  host  medium.  Since  the  effect  of  the  host  is  merely 
to  provide  an  additional  contribution  to  the  polarization,  one  need  merely  introduce 
an  additional  polarization  term  Pnr  to  describe  this  nonresonant  contribution.  The 
total  polarization  in  Eq.  (2.1)  is  then 


p=Vpm  (2.96) 

where  the  first  term,  P  ,  is  the  resonant  polarization  resulting  from  the  interaction 
of  the  wave  with  the  two-level  system  suspended  in  the  host  medium.  The  frequency 
dependence  of  the  nonresonant  polarization  is  conveniently  described  in  terms  of  a 
susceptibility  x  ((u)  by  writing 
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where  E(id)  is  the  Fourier  transform  of  the  electric  field  vector 

E (u>)  =/  dteM  E(t)  (2.98) 

In  the  neighborhood  of  the  carrier  frequency  tuQ,  the  susceptibility  may  be 
approximated  by 


2 

47rv(uj)  =  a0+a2  (cu0 /uj) 

(2.99) 

where  aQ  and  a2  are  real.  Hence  X  is  real  and  x(uj)  =  x(-tu)  which  assures  the  reality 
of  Pnr.  Note  that  any  absorption  associated  with  an  imaginary  part  of  X  is  ignored. 
The  nonresonant  polarization  then  satisfies 


An 


(2.100) 


For  an  electric  field  polarized  in  the  Y  direction  and  travelling  in  the  positive  X 
direction,  it  follows  from  the  Maxwell  equations  that  the  associated  magnetic  field 
vector  is 


where  k  is  a  unit  vector  in  the  Z  direction. 


Energy  conservation  is  expressed  by 


(2.101) 


V-(EXH)  +  ^  4-  (E2+H2)+  4Z.  E 


.  _  4tt  f  E2  ,  tJut  i£L) 

at  '  T\fE  2  at  1 


(2.102) 


where  Eq.  (2.21)  has  been  employed.  If  one  assumes  that  a  steady  state  pulse  is 
maintained  by  a  balance  between  ohmic  losses  and  resonant  gain,  then  the  right-hand 
side  of  Eq.  (2.102)  will  vanish.  The  left-hand  side  may  be  simplified  and  one  finally 
obtains 


E  • 


_Lii 

c  at 


Vxjhj  + 


Att 

C 


*5*-)  =0 

at  / 


(2.103) 


If  it  is  assumed  that  E  is  of  the  form  of  a  steady  state  pulse  with  envelope 
velocity  ve  and  phase  velocity  Vp  then  one  may  write 


r 

)£(t-x/Vp)C0S  w0(t-x/vp)+  <f>(t-x/ve) 


(2.104) 
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From  Eq.  (2.101) 


where 


vXtt'  f  +  ^ 


(2.105) 


Zc-  Zoos  <p 
ZS  =  Z  sin  <fi 


(2.106) 


4/  =  OJ0(t-X/Vp) 


(2.107) 


When  second  harmonic  terms  are  neglected. 


171- V 


1  d 


ve2  dt  ~  ~ 


and  Eq.  (2.103)  takes  the  form 


(2.108) 


2c(l  +  °o  +  &  (v'e  Vp)  C2 


(2.109) 


Since  the  two  terms  in  this  equation  have  a  completely  different  time  dependence, 
one  must  require  that  they  vanish  separately.  This  yields 


t!,S« 

Vp  ’  >/*  +a0“  y/ a2 

Turning  now  to  the  wave  equation  for  the  medium  under  consideration,  one  has 


(2.110) 


i+ao  c  1 

a*2  c2  at2  c2  ~  c2  dt  at  , 


(2.111) 


With  the  velocities  as  determined  above,  the  left-hand  side  of  this  equation  reduces 
to  a  perfect  time  derivative  and  one  finds 


oj0AE0(£scos \p  -£csin<//)  =  crE  +  yp 


(2.112) 
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where 


2?rA=  -  ^2(l+a0) 


(2.113) 

The  resonant  polarization  is  of  the  form 

V  n0^[-Psinty+<p+2cos(^+<£)]  (2.n4) 

where  \|j  is  as  defined  in  Eq.  (2.107).  Employing  the  slowly  varying  envelope  approxima¬ 
tion, 

^=-(j0n0«f>  [p  slnty+fi+Qcosty+fi]  (2.115) 

and  Eq.  (2.112)  yields 


i^yC-2) sin<£-(££-p)cosc£  =  0 
(££-p)  sin<£+(£y£-.2)cos<£=0 


(2.116) 


Two  expressions  may  now  be  formed  for  tan  cp.  When  they  are  equated,  one  finds 


(W-/3)2.-W£-2>2 


(2.117) 


where  £  is  as  defined  in  Eq.  (2.77)  and 


y=  w0A/cr 


(2.118) 


Since  all  terms  are  real,  each  side  of  Eq.  (2.117)  must  be  zero  and  hence  one  obtains 


P--IC 

2--y',e 


(2.119) 


From  Eqs.  (2.119)  it  is  seen  that  y  represents  the  ratio  of  the  in-phase  component  of 
polarization  to  the  out  of  phase  component.  For  a  dispersionless  system  y  =  0  and 
hence  2=0.  The  effects  of  dispersion  may  therefore  be  considered  to  be  large 
when  y  >  1 . 


From  Eq.  (2.39) 


*  2.  1-A2 

1  '  C2  o  +  r2) 


(2.120 
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and  from  Eq.  (2.37a) 


dA —  — £L_  / 1  _/)  2  \ 
dt  "  5Ti+7 r  } 


(2.121) 


The  population  is  again  seen  to  be  inverted  by  the  pulse  for  the  solution  of  Eq.  (2.121) 
is 


A  =  -tanh  (t/rp) 


(2.122) 


where 


tp  =  £  (l  +  y2)/ft 

The  pulse  shape,  which  follows  from  Eq.  (2.120)  is 


(2.123) 


t  -  ~^r~  ssch  f/rp 


(2.124) 


From  Eqs.  (2.88),  (2.119)  and  (2.120), 


*±z-JLA 

at 


(2.125) 


and  the  "instantaneous"  frequency  is 


w 


inst 


d±  . 

w°~af ' 


■A) 


(2.126) 


Hence,  for  y  comparable  to  unity,  the  frequency  sweep  is  proportional  to  the  first 
power  of  the  ratio  of  optical  period  to  pulse  width. 

For  large  y  it  may  be  easily  seen  that  the  population  inversion  takes  place  by 
means  of  adiabatic  rapid  passage'"^-5'.  Introducing  a  position  vector  in  a  three 
dimensional  ft  ,  2,  A  space  according  to 


r-.\P+\2+VA  (2.127) 

and  a  vector  describing  the  electric  field  and  the  detuning  of  an  individual  two- 
level  system  by 


—  rs.  *  fl  ,  f  34* 


(2.128) 
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Eqs.  (2.37),  with  f  set  equal  to  zero,  may  be  written 

— -  exR 
df  ~  ~ 

For  large  y,  the  angle  between  R  and  E  is  given  by 

COS  8  =  -Arfr  "  I 

|R||E|  y  » I 


(2.129) 

(2.130) 


Hence  the  position  vector  remains  collinear  with 

two-level  system  and  goes  from  v/t  to  -y/r  the  1c 
^  of  p  p 


1  to  -1. 


E.  As  the  pulse  passes  a  given 
component  of  R  must  proceed  from 


It  is  to  be  expected  that  future  research  in  the  field  of  ultrashort  pulses  will 
place  increasing  emphasis  upon  the  phase  characteristics  of  such  pulses. 


2.5  Transformation  of  Equations  for  a  Two-Level  System 

It  has  been  found  useful  to  observe  that  Eq.  (2.37)  are  a  set  of  scalar  equations 
which  have  exactly  the  same  structure  as  the  Frenet-Serret  equations  of  differential 
geometry(36) .  j-j-  qs  known  that  the  solution  to  such  a  set  of  equations  is  equivalent 
to  the  solution  of  a  Riccati  equation''-'  To  show  this,  one  first  recalls  from 
Eq.  (2.39)  that  an  integral  of  Eqs.  (2.37)  is 


A2+  P*+2‘-- 1 

(2.131) 

Two  new  functions  may  now  be  introduced  by  writing 

A+iP  .  \+2 

1  -2  A+i/> 

(2.132a) 

A+i  P  i+j2 _ L  -  a,* 

1  -2  A  +  i  P  + 

Equations  (2.132)  may  be  inverted  to  yield 

(2.132b) 

*.  _  _  2Re<j> 

|<£|2-h 

(2.133a) 

p  _  j  i+M  .  2lm<£ 

0_v|/  |<^j2+l 

(2.133b) 

„  I0I-1 

J  ’  t-*  ’  |0|2+i 

(2.133c) 

Equations  governing  the  time  dependence  of  9  and  'If  are  readily  deduced  by  inserting 
Eqs.  (2.133)  into  Eqs.  (2.37).  It  is  found  that  9  satisfies  the  Riccati  equation 


JI  =  +  lif  +4£)(*2-l> 


(2.13^) 


and  that  i|r  satisfies  the  same  equation. 
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One  may  now  employ  the  usual  transformation  to  convert  this  Riccati  equation  to 
a  second-order  linear  equation.  If  the  phase  term  is  neglected  the  problem  is  reduced 
to  that  of  solving  the  equation 


W  +  3-  (f2+£2+2i£)  W  =  o  (2.135) 

where  the  dot  signifies  differentiation  with  respect  to  x-  The  new  dependent  variable 
w  is  related  to  <p  through  the  transformations  w  =u  exp(-^ J*  T  dr'  Z)  and  cp  =  (2i/f) 
d(lnu)/dj.  From  well-known  properties  of  such  second-order  differential  equations,  it 
follow.,  that  in  general  it  is  impossible  to  write  cp  or  i}r  explicitly  in  terms  of 
quadratures  of  £  . 

Equation  (2.135)  is  particularly  instructive  since  it  puts  power  broadening  in 
evidence  and  provides  an  immediate  contact  with  results  obtainable  from  the  well- 
known  vector  model  for  describing  the  response  of  a  two-level  system  to  an  external 
f ield^T) .  We  now  digress  briefly  to  consider  this  approach.  For  a  constant 
envelope  Z  =  £Q/Eq.  (2.135)  is  readily  solved  in  terms  of  the  functions 


w  =  exp  f±^V f’+*V  )  (2'136) 

If  the  population  is  initially  in  the  lower  level  and  the  pulse  is  turned  on  at 
t  =  0,  then  the  proper  initial  condition  for  w  is  readily  found  to  be  w(o)/w(o)  = 
i(f  -£Q)/f2  +£2)2  and  one  obtains 

/V  =  -i+2cosz a  sirf(fr/2)  (2.137) 


where 


a  =  tan-l(f/£0) 


(2.138) 


This  constant  field  result  agrees  with  that  obtained  from  the  geometric  model.  It 
should  be  emphasized  that  although  the  vector  model  itself  is  applied  for  arbitrarily 
short  pulses,  Eq.  (2.137)  is  only  applicable  if  the  pulse  envelope  varies  slowly 
on  the  time  scale  Q-'*'.  For  the  ultrashort  optical  pulses  under  consideration  here 
this  condition  is  violated.  An  example  of  a  time  dependent  pulse  profile  for  which 
Eq.  (2.134)  is  still  soluble  in  closed  form  is  the  steady  state  solution  for  self- 
induced  transparency,  namely 


Z  -  2  sechr  (2.139) 

The  solution  of  Eq.  (2.134)  when  £  has  the  form  given  in  Eq.  (2.  39)  and  which  takes 
on  the  velue  -1  as  t  -*>  -®  may  be  shown  to  be 
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f  -  i  exp(icr/2) 
f  -  i  exp  (- 1  a/2) 


4>  =  - 


(2.l4o) 


where 


a 


=  4  tan 


T 

e 


(2.l4l) 


Eqs .  (2.133)  then  yield  expressions  for  the  response  of  the  system  which  agree  with 
Eqs.  (2.56) 


2.6  Specialization  to  a  Soluble  Model 

Although  it  is  possible  to  obtain  a  fairly  complete  analytical  description  of 
steady  state  pulse  propagation  in  an  inhomogene ously  broadened  medium,  other  very 
interesting  features  of  pulse  propagation  have  not  thus  far  been  found  to  yield  to  any 
simple  analytical  treatment.  It  has  been  noticed,  however,  that  if  inhomogeneous 
broadening  is  neglected,  the  analysis  may  be  pursued  much  further.  Fortunately,  it 
has  been  found  that  results  predicted  on  the  basis  of  such  a  model  are  preserved  to 
a  considerable  extent  when  inhomogeneous  broadening  is  included  and  the  more  complete 
set  of  equations  is  investigated  by  numerical  computations'  Furthermore,  the 

model  is  not  without  physical  interest  in  its  own  right  since,  as  mentioned  above,  it 
may  be  used  as  an  approximate  description  of  optical  pulse  propagation  under  condi¬ 
tions  of  extreme  saturation  broadening. 

The  simplification  introduced  by  the  assumption  of  vanishing  bandwidth  is  im¬ 
mediately  evident  when  one  recognizes  that  Eq.  (2.134)  becomes  linear  for  f  =  0. 

The  solution  is  then 


<p--  ±  e'a 

where  o  is  given  by  Eq.  (2.l4l)  or  equivalently, 


(2.142) 


£  =  da-/  dr 


(2.143) 


The  choice  of  sign  in  Eq. 
ditions  A  (^  -  <*>)  =  ±1. 


(2.142)  is  again  related  to  the  two  relevant  initial  con- 
From  Eq.  (2.133)  there  follows 


A=  ±cosa 

( 2 . l44a) 

P  =±  sin  a 

(2.l44b) 

2-0 

(2.l44e) 
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and  Eq.  (2.35)  takes  the  form 


d2<r 

d(dr 


±sin<r 


(2.145) 


This  nonlinear  partial  differential  equation  is  the  fundamental  equation  of  the 
model.  Fortunately,  it  has  already  been  studied  extensively  since  it 
arose  long  ago  in  the  theory  of  pseudospherical  surfaces,  i.e.,  surfaces  of  constant 
negative  curvature(3o)  More  recently,  it  has  arisen  in  the  analysis  of  many 
physical  problems'^" 


The  general  solution  of  Eq.  (2.145)  is  unknown.  However,  a  variety  of  particular 
solutions  corresponding  to  specific  pseudospherical  surfaces,  have  been  discovered. 

One  rather  large  class  of  solutions  is  expressible  in  terms  of  the  variables 


U  =  0  T  +  £/0 

v  =  ar-£/a  (2.146) 


where  a  is  an  arbitrary  constant.  In  terms  of  these  independent  variables,  Eq.  (2.145) 
takes  the  form 


d2cr  d2a  .  . 

d7-± 


(2-147) 


The  above-mentioned  solutions  are  of  the  form 


<r  ( u,v)  =  4  tanH^F(u)/G(v^ 

Substitution  of  this  assumed  form  into  Eq.  (2.147) 
F(u)  and  G(v)  satisfy  the  equations 


(2.148) 

leads  to  the  requirement  that 


F'2=  -kF4+m2F2+n 

Q>'2  =  kG*  +(m2-i)G2-n  (2.149) 


where  k,  m  and  n  are  arbitrary  constants.  The  various  psuedospherical  surfaces  cor¬ 
responding  to  such  solutions  are  known  as  the  surfaces  of  Enneper  and  have  been 
exhaustively  catalogued  by  R.  Steuerwald^^^ . 


Among  such  solutions  are  to  be  found  analytical  expressions  that  describe  not 
only  the  steady  state  2tt  pulse  associated  with  self- induced  transparency  but  also 
solutions  that  correspond  to  a  4tt  pulse  as  well  as  pulse  envelopes  for  which  the  total 
pulse  area  is  zero,  so-called  zero  n  pulses.  As  noted  previously,  the  negative  part 
of  the  envelope  in  such  a  pulse  is  merely  the  way  in  which  the  present  model  accom¬ 
modates  a  phase  change  of  tt  that  could  take  place  in  a  more  complete  theory  in  which 
the  phase  term  cp  of  Eq.  (2.22)  were  retained'0'*4^) .  The  4rr  solution  exhibits  the 
pulse  breakup  phenomenon  that  has  been  observed  both  experimentally*'0''  and  in  more 
complete  numerical  computations ( ^ ' . 
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Another,  and  more  general,  method  of  obtaining  solutions  is  to  use  the  fact 
that  Eq.  (2.145)  is  an  example  of  an  equation  which  admits  of  a  Baecklund  transforma¬ 
tion^®#*^),  This  transformation  theory  may  be  interpreted  geometrically  as  the 
transformation  of  a  surface  which  represents  a  solution  of  a  given  partial  differen¬ 
tial  equation  into  another  surface  which  is  the  solution  to  the  same  or,  in  some 
cases,  another  known  equation. 

For  Eq.  (2.145)  the  transformation  equations  are^®>  39) 


(2.150a) 

7  -°b>  = 0  (2.150b) 

One  may  easily  show  that  both  c0  and  o-^  satisfy  Eq.  (2.145).  Hence,  from  a  given 
solution  aQ  one  may  obtain  a  new  solution  which  contains  not  only  the  constant  a 
but  also  an  arbitrary  constant  of  integration.  This  transformation  may  be  used 
repeatedly  to  generate  a  solution  02  ^rom  °i  etc.  For  extensive  calculations  of  this 
sort,  it  has  been  found  convenient  to  use  a  symbolic  representation  of  Eqs .  (2.150) 
in  which  a  transformation  from  a  solution  a ^  to  a  solution  c .  with  a  constant  a^  is 
represented  as  shown  in  Fig.  2.  As  a  first  usage  of  such  multiple  transformations, 
one  may  show  quite  readily  that  the  four  solutions  that  are  related  by  the  trans¬ 
formation  depicted  in  Fig.  3  will  satisfy 


-on  (-a**)  .  .an  (^)  (2.151) 

which,  quite  remarkably,  does  not  involve  any  quadratures.  A  simple  algebraic 
manipulation  of  the  eight  equations  implied  by  Fig.  3  leads  immediately  to  Eq.  (2.151). 
It  will  be  shown  subsequently  that  this  result  may  be  used  to  construct  a  4rr  pulse  as 
well  as  a  number  of  different  types  of  zero  n  pulses. 


If  a^  =  &2  an<*  the  integration  constants  in  o^  and  Oq  are  different,  this 
relation  merely  yields 


aiz<T0±7r  (2.152) 

When  the  integration  constants  are  the  same,  the  resulting  indeterminate  form  may  be 
evaluated  from  the  usual  Taylor  expansion  and  one  finds 


tan 


Qi  ( dc\  1 

*  /"  2  \da,  dy, 


d/A 

da,  / 


(2.153) 


where  y^  is  the  constant  of  integration  arising  from  the  solution  of  Eqs.  (2.150). 
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It  will  bo  shown  that  the  compounding  of  transformations  shown  in  Fig.  4  yields 
a  6tt  pulse.  An  obvious  generalization  to  obtain  a  2m  pulse  and  undoubtedly  more 
complicated  zero  n  pulse  suggests  itself  immediately  but  the  subject  has  not  been 
pursued  beyond  this  point. 

All  of  the  above  solutions  represent  modes  of  propagation  that  are  realizable 
in  an  attenuator.  It  has  been  found  that  a  solution  which  applies  to  a  lossless 
amplifier,  the  n  pulse,  is  contained  within  the  family  of  solutions  that  results  from 
the  specification  that  a  be  of  the  form  o(5t).  Such  solutions  have  also  been  con¬ 
sidered  within  the  context  of  differential  geometry' and  have  been  used  recently 
in  the  attenuator  case  to  describe  coherent  resonance  fluorescence^ ' . 


2.6.1  Specific  Pulse  Profiles 

We  now  turn  to  a  more  detailed  consideration  of  the  various  solutions  of 
Eq.  (2.11*5)  that  were  mentioned  above.  The  results  are,  of  course,  meagre  in  com¬ 
parison  with  the  complete  analytical  description  of  the  evolution  of  arbitrary  initial 
pulse  shapes  that  can  be  obtained  when  dealing  with  linear  initial  value  problems. 
However,  the  particular  solutions  that  have  been  found  do  exhibit  many  of  the  import¬ 
ant  and  interesting  features  of  optical  pulse  propagation  as  't  has  been  observed 
both  experimentally  and  an  output  from  numerical  computations. 


2. 6. 1.1  2tt  Pulse 


As  has  been  indicated  above,  a  large  number  of  pulse  profiles  may  be  obtained 
for  propagation  in  an  attenuator.  Perhaps  the  moot  widely  Known  solution  of  this 
type  is  the  one  related  to  self- induced  transparency.  It  my  be  obtained  in  a  number 
of  ways,  the  simplest  being  that  of  assuming  a  steady  state  solution  of  the  form 
ct( t  -  x/v).  Such  steady  state  solutions  will  be  discussed  subsequently.  The  solution 
may  also  be  obtained  by  noting  that  o  r  0  in  a  solution  of  Eq.  (2.ll»5).  This  solution 
may  then  be  used  as  aQ  in  the  Paccklund  transformation  given  by  Eqs.  (2.150).  Choosing 
the  lower  sign  in  the  rc-ond  Eqs.  (2.150),  nr.  is  required  for  propagation  in  an  atten¬ 
uator,  the  two  resulting  first-order  differential  equations  have  the  solution 


cr,  *  4  ton',e‘' 


(2.15'0 


where  v  is  as  defined  in  Eq.  (2.l!»6).0nc  eons  tant  of  integration  has  been  net  equal 
to  zero  since  it  merely  serves  to  translate  the  initial  location  of  the  solution  along 
the  v  axis.  The  corresponding  electric  field  follows  from  Eq.  (2.1U3)  and  one  finds 


C  *  20  sech  ux  20  sech 


[oHO-x/v)] 


(2.155) 


where 


(2.156) 
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From  Eq.  (2.155)  one  sees  that  (afi)  determines  the  width  of  the  pulse  envelope. 
Therefore,  setting  afi  =  tp"1  where  rp  is  the  pulse  width,  the  expression  for  the 
electric  field  envelope  becomes 

=  |-sech  (t-x/v)Tp  (2.3.57) 

where 

i=~C~+  a'TP  =  T(i+um/Uw)  (2.158) 

cy '  is  as  defined  in  Eq.  (2.30)  and  Uw  and  Um  are  the  energy  densities  of  wave  and 
medium  respectively.  The  last  form  for  the  velocity  given  in  Eq.  (2.158)  is 
particularly  instructive  and  has  been  derived(52)  on  simple  physical  ground  by 
merely  equating  the  average  energy  of  both  wave  field  and  medium,  Vrp(Uw  +  Um),  to 
cr X  the  amount  of  energy  that  flows  through  the  volume  Vt  at  the  light  velocity  c. 
Eqs.  (2.157)  and  (2.158)  agree  with  Eqs.  (2.54)  and  (2.44)  in  the  appropriate  limit, 
namely  g(Auj)  -♦6(Auj)* 


2.6. 1.2  4tt  Pulse 

It  has  been  observed  both  experimentally^  and  from  machine  computations ^  ^  ^ 
that  the  combination  of  field  strength  and  magnitude  of  dipole  moment  sufficient  to 
induce  two  inversions  in  the  population  of  the  two-level  system,  a  so-called  4rr 
pulse,  does  not  propagate  as  a  single  pulse  but  rather  separates  into  two  separate 
2tt  pulses.  Such  pulse  separations  are  also  exhibited  by  the  analytical  solutions 
(17,19)  s  q>he  4n  pulse  is  obtainable  as  the  function  in  Eq.  (2.15!)  when  one 
chooses  <j0  =  If  one  chooses  the  lower  sign  in  Eq.  (2.145),  as  is  appropriate  for 
the  attenuator,  then 


where 


»  =  a.  r-(/a 

1  1  1 


(2.159a) 


(2.159b) 


The  resulting  expression  for  may  be  put  in  the  form 


cr  3  *  4  tan 


/Ojl  02 

\a,  -  a2 


) 


sin  h^~(t/|-^2) 
cosh-g-  (i/,+i/2) 


(2.160) 


For  ai  >  0,  a0  <  0  the  function  in  Eq.  (2.l6o)  varies  from  -2n  to  2n  as  t  proceeds 
from  -®  to  ®.  Since 
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Q  =  dt£  =  CT(OC)-  (-08) 

n  J-co 


(2.161) 


one  may  expect  that  the  associated  electric  field  will  correspond  to  a  4tt  pulse. 
Setting  a-jQ  =  t-^1,  -a20  -  Tp1  using  Eq.  (2.14-3),  the  electric  field  is  found 

to  be 

Eo£  a  _  (2/r,)sech  x -f  (2/t2)sechY 
*  "h  *  A  i  -  B  ( tanh  x  tanh  Y-sechxsechY) 

(2.162) 

where  A  =  (r2z  -t,2)/(t22  +  r,2) 

B  =  2t,  t2/(t22  +  t,2) 
x  =  (t  -x/v,  )/r, 

Y  =  (t  -x/v2)/r2 

(2.163) 


and  the  velocities  V 


1 


and 


are 


ft2r,2e) 


+  a'r,2 


( 1  +  £l2  T22e  ) 


a'  r22 


(2.164) 


A  graph  of  Eq.  (2.162)  is  shown  in  Fig.  5-  As  the  pulses  become  completely  separated 
Eq.  (2.162)  reduces  to 


=  -f^sech (x  ±/3)  +  ^-sech(Y  ±  /3) 


(2.165) 


where  the  upper  sign  is  to  be  used  for  and  the  lower  sign  for  and  where 


(3  -  tanh-1  B  (2.166) 

In  order  to  obtain  a  pulse  envelope  that  begins  at  5  =  0  with  only  one  peak  as 
a  function  to  time,  one  must  impose  the  requirement  5^ £ <0at§=T=0.  This 
condition,  along  with  the  requirement  £  >  0  leads  to 


(|_  *)(|  +  4  2  -ll  )  >  0  ,  *  =  T,  /T2 


which  is  equivalent  to 

\  (3-/5  )<  r,  /r2  <  -^(3  +^5’) 


(2.167) 


(2.168a) 
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Figure  5  could,  of  course,  be  continued  back  to  negative  values  of  §  to  provide 
an  example  of  the  envelope  distortion  that  takes  place  when  an  ultrashort  pulse  over¬ 
takes  a  slower  pulse  and  passes  through  it. 


In  addition  to  the  above  results  for  specific  pulse  profiles,  certain  conclusions 
relating  to  4rr  pulses  of  arbitrary  initial  shape  may  also  be  given.  These  results 
are  obtained  by  employing  a  technique  previously  applied  to  the  Korteweg  de  Vries 
equation^)  which  has  solutions  with  many  of  the  same  pulse  like  properties  as 
Eq.  (2.145).  It  will  be  discussed  further  in  Section  2.6.3. 


One  begins  by  noting  that  Eq.  (2.145)  possesses  the  two  conservation  laws 


dr 


+ 


( 1  -  COS  cr) 


0 


(I  -  c°s  a )  + 


(2.168b) 


which  follows  from  Eq.  (2.l45)  after  multiplication  by  and  Srr/dT  respectively. 

Integrating  over  all  §  and  over  t  from  t  =  0  up  to  an  arbitrary  final  time  one 
finds 


4C, 


d£(i  -cosa) 


4CZ 


(2.168c) 


where  c^  and  c2  are  positive  constants  that  are  obtained  by  evaluating  Eq.  (2.l63c) 
for  a  given  initial  form  of  o.  After  the  pulses  have  completely  separated  one  may 
write 


cr 


=  X  4  tan-1  (exp  [ai  (t  - 

i  =  I » 2  1  L 


Tj>  -£/ai]} 


(2.l68d) 


Eqs.  (2.l68c)  then  yields  the  paid  of  algebraic  relations 


Oj-'+Of1  =  C, 

0  1  +  0  2  ”  C  g 


(2.168b) 


Introducing  the  positive  quantities  a  =  a^,  cy2  =  -a2,  one  finds  that  the  final 
amplitudes  are  given  by 


a 


1,2  = 


j(cz  ±  y^-ACz/C,  ) 


(2.l68f ) 


Also,  one  notes  that,  since  the  <y ^  are  real 
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C,  C2>4 


(2.i68g) 


(78 ) 

It  has  been  shown'  that  Korteweg  de  Vries  equation  satisfies  an  infinite 
number  of  conservation  laws.  For  that  equation,  the  above  technique  may  be  extended 
to  treat  pulses  which  evolve  into  any  number  of  isolated  pulses.  Thus  far,  no 
additional  conservation  laws  have  been  found  for  Eq.  (2.145)  although,  as  will  be 
shown  below,  it  is  known  to  have  solutions  that  break  up  into  more  than  two  pulses. 


2. 6.1. 3  Zero  rr  Pulses 

As  mentioned  previously,  Eq.  (2.145)  also  admits  of  solutions  for  which  the 
associated  electric  field  envelope  becomes  negative.  Such  solutions  cannot  be  dis¬ 
carded  on  any  physical  basis  since  a  negative  field  envelope  merely  corresponds  to 
a  phase  change  of  tt  in  the  carrier  wave. 

Three  different  types  of  zero  rr  pulses  have  been  constructed  from  the  solutions 
mentioned  above.  The  simplest  type  is  obtained  by  merely  choosing  a^  >  0  in  the 
previous  solution  for  the  4tt  pulse.  The  electric  field  envelope  in  dimensionless 
f orm  is 


t  = 


_ Oi  sech  I/,  -  og  sech  vz _ 

1  -  0  ( tanh  v,  tanh  uz  +  sech  1/,  sech  vz ) 


(2.169) 


where 


A  =  ( a |Z  -  a2z  )/(a,z  +  a2z ) 


(2.170) 


B  s£o,  az/(a,z  +  a2z) 


(2.171) 


An  example  of  this  result  is  shown  in  Fig.  6. 

In  the  limit  a^  =  a^,  Eq.  (2.151)  becomes  indeterminate.  In  this  case,  with 
again  set  equal  to  zero,  one  may  use  Eq.  (2.153)  and  obtain  as  the  second  form  for 
a  zero  tt  pulse, 


a  =  4 tan'1  ( **  sech  v) 


(2.172) 


This  yields  the  field  envelope 
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z 


-  4a  sech  vi— — fr  anLg— ' ) 
\i  +  uzsechzi/ ) 


A  graph  of  this  result  is  shown  in  Fig.  7* 


(2.173) 


The  third,  and  by  far  the  most  interesting  type  of  zero  tt  pulse  is  obtained 
by  allowing  the  parameters  a^  and  ag  in  Eq.  (2.151)  to  become  complex  and  requiring 


a,  =  c*  =  a  =  a  +  i/3 


One  then  finds 


a  - 


4tarf 


fa.  sing  \ 
\/3  cosh  p  j 


where 


p  =  a(r-£/lalz) 


and 


q  =  0(t+  £/lalz) 


(2.174) 


(2.175) 


(2.176) 


(2.177) 


The  electric  field  envelope  is 


z 


4  a  sech 


P 


cos  q  -  (a/P)  sin  q  tan  h  p  ' 
!+  (a//3)zsinzq  sechz  p 


A  graph  of  this  result  is  shown  in  Fig.  8. 


(2.178) 


Unlike  the  two  previous  types  of  zero  rr  pulses,  the  envelope  given  in  Eq. 
(2.174)  tends  to  remain  as  a  single  localized  disturbance.  It  provides  an  alternate 
and  more  flexible  form  of  self- induced  transparency.  It  has  been  found(53)  from 
numerical  computations  that  this  pulse  shape  is  remarkably  insensitive  to  variations 
in  inhomogeneous  broadening. 


2. 6. 1.4  6tt  Pulse 

As  mentioned  above,  the  6rr  pulse  is  obtained  from  the  sequence  of  transforma¬ 
tions  depicted  in  Fig.  4.  From  this  diagram,  the  corresponding  analytical  expressions 
are  easily  seen  to  be 
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where 


<7 


crz  *  4 ran’1  ^K|3tan 


o~a~  o~b 
4 


-I  V 2 

cr2  «  4tan  e 


crQ  =  4  tan 


-i 


sinh^  (»/,  -  i/2) 
_  12  cosh-jj  (j/(  +  v2) 


(2.179) 


and. 


-i 

4  tan 


sinh  ^  lv2" 
cosh^-  U/2  +  i/j) 


K|j  =  (aj+aj)  /(a;  -aj ) 


(2.180) 


(2.l8l) 


One  may  immediately  impose  a  number  of  constraints  upon  the  triad  of  constants 
a-^,  a2,  a^.  In  the  first  place,  for  the  envelope  function  corresponding  to  a2  to  be 
positive,  one  must  require  a2  >  0.  One  then  proceeds  to  make  cra  a  4rr  pulse  which 
requires  a-^  <  0.  Also,  is  made  a  zero  tt  pulse  which  requires  0  <  a^  <  a^.  The 

three  constants  a^  may  be  related  to  the  widths  of  the  three  pulses  when  complete 
separation  has  taken  place  by  setting  -a-jfl  =  Tq,  agfi  -  t 2>  a^Q  =  t^. 

As  with  the  4n  pulse,  one  must  impose  additional  restrictions  in  order  to 
assume  a  pulse  shape  that  consists  of  a  single  peak  at  §  =  0.  The  inequality  is  much 
more  complicated  in  this  case  than  in  Eqs.  (2.167)  and  has  not  been  analyzed  in 
detail.  By  a  trial-and-error  method  the  case  shown  in  Fig.  9  has  been  obtained. 


2. 6. 1.5  tt  Pulse 

To  describe  a  rr  pulse,  one  must  turn  to  the  other  completely  different  class 
of  solutions  of  Eq.  (2.145)  that  were  mentioned  above.  These  are  solutions  in  which 
the  independent  variables  occur  solely  in  the  product  form  §t.  When  one  sets  Z  = 
in  Eq.  (2.1^5)  it  reduces  to  the  ordinary  differential  equation' '  ' 


Z a"  +  a'-  sincr -  0  (2.182) 

where  the  prime  indicates  differentiation  with  respect  to  Z. 

A  new  dependent  variable,  w,  related  to  a  by  w  =  exp(to),  may  be  shown  to 
satisfy  a  special  case  of  the  equation  that  defines  the  third  Painleve  transcendent' '  ' 
Since  the  properties  of  such  functions  are  essentially  unknown  and  tables  appear  to  be 
unavailable,  it  seems  preferable  to  resort  to  a  direct  numerical  integration  of 
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Eq.  (2.182).  The  result  of  such  a  numerical  solution  is  shown  in  Fig.  10  which  also 
includes  the  result  for  a'  =  ( £ /§)  as  well  as  a  phase  plane  diagram  of  the  solution. 
The  'example  shown  in  Fig.  10  satisfies  the  condition  a(0)  =  0.1  as  well  as  <r(0)  = 
sino(o)  which  is  required  for  continuity  at  the  origin. 

Scaling  laws  for  tt  pulse  propagation  may  be  inferred  from  these  results. 

Since  the  abscissa  for  the  pulse  envelope  is  §t ,  the  actual  pulse  envelope  narrows 
linearly  with  increasing  distance  of  propagation.  Also,  since  £  =  §a',  the  amplitude 
of  the  envelope  increases  linearly  with  distance.  This  spatial  evolution  of  the  pulse 
shape  is  shown  explicitely  in  Fig.  11.  Similar  results  have  been  obtained  from 
direct  numerical  analysis^"^ 1  'of  the  partial  differential  equations  governing  optical 
palse  propagation  in  a  resonant  medium.  A  comparison  of  Fig.  11  with  results  present¬ 
ed  in  Ref.  16  shows  that  until  the  signal  becomes  so  large  that  the  linear  loss  is 
dominant,  neglect  of  the  loss  term  introduces  no  significant  change  in  propagation  in 
an  amplifier. 

The  fact  that  the  self-consistent  interaction  of  field  and  resonant  matter 
should  give  rise  to  ringing  is  not  unexpected  in  view  of  the  known '  response  of 
an  inverted  population  to  a  specified  spatial  mode  of  the  electric  field.  The  ring¬ 
ing  may  also  be  inferred  from  a  theorem  concerning  solutions  of  Eq.  (2.14-5).  It  may 
be  shown^9 /that  there  is  no  function  which  satisfies  Eq.  (2.145)  and  at  the  same  time 
remains  within  the  interval  0  <  0  <  tt 


2.6.2  Steady  State  Solutions 

An  example  of  a  steady  state  solution  has  already  been  given  with  the  discus¬ 
sion  of  the  2tt  pulse.  This  solution  is  actually  a  limiting  form  of  a  more  general 
oscillatory  solution  which  is  now  considered.  Similar  resultg  for  propagation  in  an 
inhomogeneous ly  broadened  medium  have  also  been  reported >  9  . 

Steady  state  solutions  will  be  functions  of  one  of  the  variables  defined  in 
Eq.  (2.145).  Choosing  the  variable  to  be  v,  one  readily  shows  that  the  conservation 
laws  given  by  Eqs .  (2.38)  and  (2.39)  take  the  form 

-  /+  A  -- Ai  (2.183a) 

P  +  A  :  1  (2.183b) 

where  A  is  a  constant  of  integration  and  the  constant  a  is  again  (Qt  )"^.  Allowing 
for  a  steady  state  solution  in  which  £  is  non-zero  when  the  entire  population  is  in 
the  ground  state,  one  sees  from  Eq.  (2.l83b)  that  the  constant  A  may  be  less  than 

-1. 


From  Eq.  (2.37a)  and  (2.183)  one  readily  obtains 
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4^-)  =  2(A-A,Hi-Az) 

dv;  (2.134) 

which  shows  that  A  may  be  expressed  in  terms  of  elliptic  function. 

If  -kA^l,  a  solution  for  which  the  population  varies  between  A.  and  1  is 
given  by 


which  leads 


I 


dA 


TV  V(/V-A,)( i -A2) 


A  =  1-2  k2sn2' 


iv-V0),^ 


where  k 


i(l  -A1)-  From  Eq.  (2.l83a), 


(2.185) 


(2.186) 


£=  2akcn[(v-v0),k]  (2.187) 

A  solution  for  which  -l<A<Aq  could  also  be  given,  but  it  requires  that  a 
be  imaginary.  This  implies  envelope  function  propagation  faster  than  the  light 
velocity.  In  the  limit  k  -*■  1  this  solution  goes  over  to  one  which  represents  2n 
pulse  propagation  in  an  amplifier  which  is  unstable. 


For  A 1  <  -1  it  is  seen  that  k  >  1. 


Using  the 


relations 


(56) 


sn(v,k)  =  kHsn(kv,  k-1) 
cn(v,k)=  dn(kv,  k"1) 


(2.188) 


one  finds  that  the  population  difference  and  field  envelope  may  be  written 


A  =  i-2sn[k'l(v-v0),k]  (2.189a) 

C-  2ak’'dn  [k"'(v-v0),  k]  (2.189b) 

p 

where  now  k  =  2/(l  -  A-,).  These  latter  forms  may,  of  course,  be  obtained  by  direct 
integration  of  Eq.  (2.1o4).  In  the  limit  A]_  -*■  -1  both  solutions  reduce  to  that  for 
the  2tt  pulse  in  an  attenuator. 

(57) 

It  has  been  conjectured  that  these  steady  state  solutions  may  be  realized 
in  self-pulsing  situations. 
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2.6.3  Relation  to  the  Korteweg-deVries  Equation 

(  Sfi'l 

It  has  been  noted'  '  that  the  hyperbolic  secant  solution  of  Eq.  (2. 145)  and 
the  decomposition  of  pulses  into  a  sequence  of  such  "solitary  waves"  is  very  similar 
to  results  obtained  in  recent  investigations  of  the  Korteweg-deVries  equations (58-63) t 
In  fact,  one  may  show  that,  for  steady  state  solutions,  the  square  of  the  envelope 
function  £  satisfies  the  corresponding  steady  state  Korteweg-deVries  equation.  This 
may  be  seen  quite  readily  by  writing  o(§,t)  in  the  form  a(v)  where  v  =  t  -  §•  A 
first  integral  of  the  resulting  equation  is  readily  obtained.  The  integral  satisfying 
£,  (-oo)  =  o(-oo)  =  0  is 


COS  a 


(2.190) 


When  this  result  is  solved  for  a  and  differentiated,  one  obtains 


(2.191) 


where  u  £  .  Two  derivatives  of  this  equation  yields 


The  substitutions 


—  3u  du 
dv 


1/2 


X 


(2.192) 


(2.193) 


transforms  Eq.  (2.192)  to 


-kf'+ff'+f"'=0  (2.!94) 

where  the  prime  indicates  differentiation  with  respect  to  x.  Eq.  (2.194)  is  the 
steady  state  form  of  the  Korteweg-deVries  equation.  The  solution  that  vanishes  for 
large  values  of  x  may  be  written  in  the  forrr/^) 


f  =  3ksech2  (2.195) 

which  is  readily  shown  to  be  equivalent  to  the  result  given  in  Eq.  (2.157)*  Periodic 
solutions  in  terms  of  elliptic  functions  have  also  been  given  .  They  are  related 
to  the  result  given  in  Eq.  (2.187)  and  are  referred  to  as  cnoidal  waves. 
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(77) 

Although  the  multisoliton  and  oscillatory  solutions  of  the  Korteweg-deVriec 
equation  are  similar  to  results  obtained  above  for  2nrr  and  tt  pulses  respectively,  and 
the  criterion  given  in  Eq.  (2.167  is  similar  to  oneft appearing  in  the  breakup  of  two 
soliton  solutions  of  the  Korteweg-deVries  equation^  ,  to  the  authors  knowledge  no 
quantitative  relation  between  Eq.  (2.145)  and  the  time  dependent  Korteweg-deVries 
equation  has  been  discovered. 


2.6.4  Stability  Considerations 


When  inhomogeneous  broadening  is  present,  the  stability  of  the  area  under  the 
electric  field  envelope  £  (x,t)  may  be  inferred'  *  from  the  solution  of  Eq.  (2.68). 
However,  an  integration  over  the  frequency  of  the  detuning  associated  with  inhomo¬ 
geneous  broadening  is  a  crucial  step  in  the  derivation  of  this  result.  In  the  model 
being  considered  here,  inhomogeneous  broadening  is  neglected  and  so  one  must  rely 
upon  other  considerations  to  infer  area  stability.  This  is  accomplished  quite 
readily  by  noting  that  Eq.  (2.145)  may  be  written 


d£  _  flef 


d  x 


P 


sincr 


(2.196) 


where 


er(x,t)=<r(x,-*  )+/*  dt'£(x,t') 


(2.197) 


For  a  system  initially  in  the  lower  level,  one  may  take  a(x,-®)  =  -tt  for  then 
A  (x,-®>)  =  coso(x, -'»)  =-l.  For  a  system  initially  in  the  upper  level  one  may  assume 
rr(x,  -oo )  =  0. 

For  the  hyperbolic  secant  pulse  envelope  given  in  Eq.  (2.155)  9  =  (lP/4  )j  ^dt^ 
=  2n  so  that  near  the  trailing  edge  of  this  pulse,  Eq.  (2.196)  goes  to 

jf  --  $-a'sin[(or)  +  Zlr]  (2.198) 

where  the  upper  choice  is  made  for  the  attenuator  and  the  lower  choice  for  the  ampli¬ 
fier.  Now  if  there  is  a  perturbation  in  £  such  that  the  total  area  0  is  greater  than 
2tt,  then  in  an  attenuator  9  <£/3x  ~  sin(n  +  e)  <0.  The  field  at  the  trailing  edge 
therefore  tends  to  decrease  to  recover  a  total  area  of  2n.  On  the  other  hand,  if  the 
perturbation  is  such  that  9  is  less  than  2tt,  then  d  <£7dx  >  0  and  the  field  at  the 
trailing  edge  increases.  The  total  area  of  such  a  pulse  therefore  tends  to  remain 
at  2rr.  In  the  amplifier,  the  inequalities  are  reversed  and  the  hyperbolic  secant  no 
longer  represents  a  stable  pulse  envelope.  These  results  are  in  agreement  with  those 
previously  obtained''*0'  for  the  case  in  which  inhomogeneous  broadening  is  included. 
The  above  considerations  predict  only  the  area  stability  and  leave  open  the  question 
of  perturbations  in  which  the  total  area  remains  unchanged.  We  now  take  up  this  topic 
and  show,  by  exhibiting  a  Liapunov  functional  with  vanishing  derivative,  that  in  the  at¬ 
tenuator  the  pulse  shape  is  stable  but  not  asymptotically  stable,  i.e.,  perturbations 
remain  finite. 
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(64) 

An  appropriate  Liapunov  functional 


F(u) 


is  given  by 


F'ul=0''[(-fv),  +  (^)2+  *"-“H  (2.l99) 

which  is  proportional  to  the  total  energy  residing  in  field  and  medium.  Differentiation 
with  respect  to  u  and  a  subsequent  partial  integration  yields 


=  />£( 


a2cr 

au2 


+  siner  + 


da_ 

du 


CO 


-co 


(2.200) 


which  vanishes  by  virtue  of  Eq.  (2. 147) (in  which  the  lower  sign  has  been  chosen  as  is 
appropriate  for  an  attenuator)  and  the  boundary  condition  that  a  represent  a  pulse  and 
hence  da/d u  and  da/dv  must  vanish  at  v  =  iso. 


Since  dF/du  is  merely  zero  rather  than  negative  definite,  it  is  not  unexpected 
that  a  first-order  perturbation  analysis  of  Eq.  (2.147)  will  contain  a  zero  eigenvalue. 
This  is  readily  seen  to  be  the  case.  Setting  a  =  a  (v)  +  a  ^  (u,v),  one  finds 
that  satisfies 


<32  cr(l) 
<3u2 


d2*{,) 

~d^~ 


-(l-2sech2  v)cr(l) 


0 


Expressing  a 


(1) 


in  the  form 


(2.201) 


cr(l)(u,v)  =  V  (v)esu 

(2.202) 


V(v)  is  found  to  satisfy  the  Schroedinger  equation 


V"  +  (X-2sechzv)V=  0 


(2.203) 


where  X  =  -(S^  +l).  For  X  =  -1  (  and  hence  S  =  0)  one  readily  finds V(v)  =  sechv 
which  is  the  solution  corresponding  to  the  expected  zero  eigenvalue. 


2.7  Inhomogeneous  Broadening  and  Photon  Echo 

A  very  extensive  analysis  of  the  effects  on  pulse  shapes  of  both  homogeneous  and 
inhomogeneous  broadening  has  been  carried  out  numerically'  ’  ' .  Thus  far  only  the 
steady  state  pulse  shape  has  been  described  analytically  when  inhomogeneous  broadening 
is  present.  However,  if  one  foregoes  consideration  of  the  actual  pulse  shapes  and 
confines  attention  to  the  area  under  the  envelope  of  the  pulse,  then  further  progress 
may  be  made.  In  particular,  it  has  been  shown  (^)  that,  a  very  simple  description  may 
be  given  of  the  space-time  evolution  of  the  photon  echoes  that  may  appear  behind  two 
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optical  pulses  no  they  propagate  through  a  resonant  inhomogeneounly  brondenec  medium. 
Thin  was  carried  out  by  noting  that  the  area  theorem,  Eq.  (2.68)  ic  still  satisfied 
if  the  pulses  are  assumed  to  be  infinitely  narrow,  i.o.,  of  the  form 


CUj)  --  eu)  8(r) 


(2.20»*) 


where  6(t)  in  a  delta  function.  The  assumption  of  propagation  at  the  light  velocity 
is  consistent  with  that  of  zero  pulse  width  according  to  Eq.  (2.57n) 


The  apparent  inconsistency  of  using  a  delta  function  in  the  slowly  varying 
envelope  f*  (?,t)  does  no  violence  to  the  theory.  It  merely  provldeo  a  convenient 
device  for  obtaining  nolutionc  to  Eqn.  (2.35)  to  (2.37)  in  the  ehort  pulce  limit. 
We  now  give  a  derivation  of  the  area  theorem  for  delta  function  pulocn  nnd  show 
how  some  of  its  implications  may  be  readily  explored. 


The  response  of  the  system  ic  governed  by  the  Riccati  equation  given  in  Eq. 
(2.13*4).  Introducing  the  new  complex  function  u  ■  u  *  Ui,  by  the  definition 
®  •  cxp((u)  Eqn.  (2.13M  and  (2.133)  becoec 


if  sin/i  =  C 

(2.205) 

P  *  sech  /x,  sin/** 

(2.206) 

2  e  -tonh/x, 

(2.207) 

Af  s  S€Ch/i|  cos/x* 

(2.208) 

Using  the  firm  for/*  given  in  Eq.  (2.20*»)  and  integrating  Eq.  (2.205)  ncrons  the 
singularity  ni  t  «  0,  one  1  «ndn 


=  9(C) 


(2.209) 


Glnce  «  ic  real,  u,  ir.  continuous  across  the  pulse  and  according  to  Eq.  (2.207 ),2 
ir  also  continuous  across  the  pulse.  The  change  in  population  ic 


afiisN>_N<  s  sech  (cos/x>  -  cos /x<) 


(2.210) 


which  may  be  vrltten  as 


Ar5,  s  A<cos8  -fe  sin  e 

1*3 


(2.211) 
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Similarly,  Eq.  (2.206)  yields 


P>  =  P<  cos 9  +N<sin  6 


(2.212) 


Eqn.  (2.211),  (2.212)  and  the  continuity  of  .2  across  the  pulse  may  be  summarized  in 
the  vector  form 


0 

I 

0 


(2.213) 


The  3x3  mtrix  represents  a  rotation  about  the  .2 axis  by  an  angle  q  and  may  be 
represented  symbolically  by  R ^2  (0 )  *  Considering  00  the  three  components  of 

a  vector  P,  Eq.  (2.213)  may  be  written 


(2. 214) 


While  the  pulse  is  not  acting,  the  system  evolves  according  to  the  homogeneous 
counterpart  for  Eq.  (2.205)  which  has  the  solution 


e,M  =  icot  Tp(fr  +  a) 


(2.215) 


where  a  •  a  ■*  Ur  is  a  constant  of  integration.  It  is  now  a  simple  matter  to  show 
that  r  1 


(p(r)  *  secha,  sin(fr  ♦  a4)  =  p(r0)cosf  (r -r0) +  £(10) sinf (r -r0) 


(2.216a) 


j?(r)  s  secha,cos(fr  ♦  a*)  =5fro)cosf  (r-r0)  p  (r0)sinf  (r- r0) 


(2.216b) 


t){t)  *  tonh  a,  =  77(t0) 

(2.216c) 

The  matrix  representation  of  this  result  in 

fP[r)\  I  cost (t-t0)  smf(r.r0)o\  fP(To)\ 

=  1-sinf  (r-r0)cosf(r-r0)o|  I  ^(r0)  j 

\i?(t)/  \  0  01/  V7(fo)/  (2.217) 


which  represents  a  rotation  through  an  angle  T(t  -  t  )  about  the  h  axis  and  may  be 
written 
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If  a  transverse  relaxation  time  were  retained  in  the  analysis,  so  that  between 
pulses  /^and  2  satisfied 


d2 


then  Eqs.  (2.218)  would  be  replaced  by 

p(T)  =  e_(t,‘r°)/IiT2  Rt?(T-T0)p(r0) 


II 

(2.219a) 

nVr-2  =  ° 

(2.219b) 

(2.220) 


At  a  time  r  after  interaction  with  a  pulse  of  area  0,  the  state  of  a  system  that 
was  initially  in  the  lower  level  is  given  by 


(-sin®  cosfr\ 

Sine  sinfr 

-cose  /  (2.221) 

Eq.  (2.35),  with  g(Au>)  ■  g(0)  to  accommodate  all  spectral  components  of  the  delta 
function,  takes  the  form 


5  -wg(0)8(r)sin® 

04 


which  yields  the  area  theorem  in  the  fona 


(2 .222) 


=  -Sin®  ,4'=Trg(0)C  (2.223) 

Thin  scheme  may  now  be  used  repeatedly  to  describe  the  response  of  the  medium 
to  a  sequence  of  pulses.  Th*'  response  due  to  two  pulses  of  area  and  Qg  a  time  T 
apart  (in  units  of  is  found  to  have  n  contribution  at  t  *  2T.  Evaluating  P 

Just  beyond  thin  time  t  ■  2T,  one  finds 


PIT)  =  R,,(t-2T)R.  (®»)Rt,(T)R  (®a)R,iT)R..J®,)P(0) 


(2.22k) 


Carrying  out  the  indicated  multiplications  for  a  system  that  in  initially  in  the 
ground  state,  one  obtains 

P(C,t)  s  -sin®,  cost t -sin®, cos®, cost (t-T) 

-  [sin®, cos G,  cos  ®,  -  cos®,  sin*(®,  /2 ) sin®,]cos  f  (r-  2T)  ( 


k$ 
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If  a  relaxation  time  T2  were  included,  this  entire  expression  would  merely  be 


multiplied  by  e 


-2T 


Ul 


Substituting  this  result  into  Eq.  (2.35)  yields 

|fr  =  -sinfi, 

=  -  sin02  cos0, 

=  -sinfl,  cos 02 cos 0,  +  cos0,  sin2(02/2) sin0, 

Eqs .  (2.226)  and  (2.227)  have  solutions 

e 


(2.226) 

(2.227) 

(2.228) 


tan-^r  =  exp(-£'  +  a) 


(2.229) 


where 


fan  -  /3sech (£'-  a )  =  >3 sin  0, 


(2.230) 


e° 


.  0  =  tan  -i2~  esc  0,  (0) 


(2.231) 


How,  if  9^(0)  *  tt/2,  92(0)  n  n,  the  optimum  cace  for  photon  echo  experiments, 
then  |B|  »  1  and  from  the  solution  for  02  one  5002  that^g  remains  nearly  equal  to 
Its  initiul  value  until  0^  decays  to  a  value  equal  to  B-1.  Until  this  final  state 
In  the  pulse  evolution  is  reached,  one  may  set  92  «  n  in  Eq.  (2.228).  The  resulting 
equation  may  then  be  transformed  to 


+  cos0,  +  cote,  sine, s  0 

0C7| 


(2.232) 


Upon  substituting  y  ■  tan(9^/2)  this  becomes  a  Riccati  equation  which  may  be  converted 
to  a  second-order  linear  equation  by  the  substitution  y  =  -2(du/do)/u.  Setting 
k  *  r.in(91/2)  one  finally  obUiins 


k(k*-i)u"  +  ( 3k 2  —  i)u'  +  ku  s  0 


(2.233) 


where  the  prime  indicates  differentiation  with  respect  to  k.  Eq.  (2.233)  has  the 
solution'** ) 


u  r  0K(k)  +  bK(k') 


h6 


(2.234) 
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where  K(k)  is  the  complete  elliptic  integral  of  modulus  k  and  k'  is  the  complementary 
modulus.  Finally,  the  solution  for  may  be  written 


63  _  (kVk)B(k')  —  (k /k')B(K) 
2  "  K(k)  +  K(k’) 


(2.235) 


where  B(k) 


is  a  tabulated  function 


(67) 


related  to  the  complete  elliptic  integral  by 


k  dk 


(2.236) 
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2* 


Figure  12  contains  a  graph  of  Eq.  (2.235)  as  well  as  the  variations  in  0^  and 
It  is  seen  that  the  amplitude  of  the  echo  increases  at  a  rate  approximately 

This  has  been  observed  experi- 


equal  to  that  at  which  the  first  pulse  decreases. 
mentally'1^ . 


2.8  Level  Degeneracy 

(8.18) 

It  has  been  pointed  out  that  pulse  propagation  under  conditions  which 

prevail  experimentally  may  lead  to  considerably  different  results  from  those  predict¬ 
ed  here  since  level  degeneracy  may  be  present.  It  has  also  been  shown  that  level 
degeneracy  has  a  marked  effect  on  the  direction  of  polarization  of  the  electric  field 
vector  of  the  echo  pulse  in  a  photon  echo  experiment'2^.  Thus  far,  however,  only 
the  source  term  for  the  echo  pulse  has  been  calculated  when  degeneracy  is  present. 

No  consideration  has  as  yet  been  given  to  the  complete  problem  in  which  the  spatial 
evolution  of  the  photon  echo  is  followed  in  the  presence  of  level  degeneracy.  Hence 
this  latter  topic  will  not  be  pursued  here. 

To  avoid  detailed  consideration  of  specif  ic  molecular  modes,  level  degeneracy 
will  merely  be  expressed  in  terms  of  a  simple  Jm  scheme.  The  two  states  previously 
denoted  by  a  and  b  arc  now  characterized  by  angular  momentum  quantum  numbers  J'm' 
and  Jm  respectively.  Additional  quantum  numbers  that  would  be  associated  with  molecular 
vibration  will  commute  with  angular  momentum  operators  and  may  be  ignored.  Each 
element  in  the  2  x  2  matrices  of  Eq.(2.10)  now  becomes  a  (2J  +  l)  x  ( 2J *  +  l)  submatrix 
itself  with  elements  <J'm'  |jm>.  As  ic  well  known'°°’  transitions  in  J  are 

restricted  to  AJ  =  J'  -  J  *-1,0,1,  the  three  alternatives  frequently  being  referred  to 
as  P,  Q  and  R  branch  transitions  respectively.  In  addition,  if  the  quantization  axis 
is  alligned  parallel  to  the  electric  field  polarization  vector  then  only  the  IP  „ 
matrix  elements  need  be  calculated.  All  such  matrix  elements  vanish  unless  m*  =  m 
One  then  finds  that  nj  -  O’"  1*1  J”>  -  where  ^  ic  the  largest  value  of 
in  each  of  the  three  cases  and 


Km  =  yi2-rr^j  ,  Aj  =  -1 


(2.237a) 


Km  =  m/j  ,  Aj  =  0 


(2.237b) 
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Km  =  y(j+i)2  -  m2/(  j  +  i)  ,Aj  =  i 


(2.237c) 


Since  the  submatrices  of  1/  =  -E  *  IP  =  -E  fP  are  now  diagonal  in  m,  the  various  pairs 
of  levels  designated  by  different  m  values  are  not  coupled  by  the  interaction  and  may 
be  treated  separately. 


Hence,  for  each  value  of  m  one  may  write 


pm  +  o,ob2  pm  T)m  e  |<  jm/pz  /  j 'm  > |‘ 


(2.238) 


V 


2E 


Pm 


1Wob  rm  (2.239) 

Assuming  that  all  sublevels  of  the  lower  state  are  equally  populated  initially. 


^m  = 


^0  /  ^  .» 


2  j  + 1 


(<  jmlP!  jm>  -  <  jm|P|  jm>) 


(2.24o) 


The  electric  field  is  governed  by 


if-  ♦  ■  ^/>«pT 


(2.241) 


where  J<  refers  to  the  lesser  of  J  and  J'.  When  this  relation  is  integrated  over  the 
duration  of  the  pulse,  one  obtains,  in  analogy  with  the  derivation  of  Eq.  (2.68) 


dx  =  2(2 jQ+  I )  sin(k"0)  (2.242) 

where 

0lx)=  dt'£(x,t') 

(2.243) 

For  transparency  to  tukc  place  it  is  necessary  that  the  right-hani  side  of  Eq. 
(2.243)  vanish.  For  Q  branch  transitions  this  will  be  possible  for  9  =  2nn  Just  as 
in  the  nondegenerate  c,  so  since  the  various  Km  are  integrally  related.  For  P  and  R 
branch  transitions,  however,  the  irrational  ratios  of  the  various  prevent  a 
simultaneous  vanishing  of  all  K  except  in  the  few  cases  in  which  there  is  only  one 
nonvanishing  value  of  K^.  For  !P  and  Q  branch  transitions  this  takes  place  for 

J  -  0,  l. 
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However,  it^has  been  noted.'  that  the  right-hand  side  of  Eq.  (2.243)  will 
also  vanish  if  f  dt'/T  =  0  •  Such  zero  tt  pulses  should  exhibit  transparency  indepen- 

J -X 

dently  of  the  values  of  the  K,n .  Although  profiles  of  zero  tt  pulses  have  been  described 
in  Section  2.6.1,  it  should  be  emphasized  that  they  have  been  obtained  for  a  simple 
nondegenerate  system  and  are  not  directly  applicable  to  the  present  situation. 


For  large  values  of  j  the  summation  may  be  approximated  by  an  integration. 
Setting  m  =  j  coney,  the  results  quoted  in  Eqs.  (2.237)  may  be  replaced  by  continuous 
variable  K  given  by 

f  sina  ,  Aj  =  ±  l 

K  M 

(.cos  a  ,  Aj  =  0  (2.244) 

For  A.'1  =  0 

y-yj-E  Km  sin(Km0)  — •>  y  J” da  sin  a  cos  a  sin(0cosa) 


=  (sin 9  -  9cos9)/92  =  j ,(0) 


where  J  (q)  is  a 


spherical  Bessel  function 


(72) 


(2.24J,) 


ypy  Z  Km  sin(Km0)  — ►  y  J  wda  sinza  sin(0sin  a) 
=  y[Ho(0) -H,(0)/0] 

(69) 


where  the  H  (h)  are  Struve  function; 
n 


(2.246) 


The  two  forms  for  the  urea  theorem  in  the  presence  of  deg^racy  for  large  j  are 
now  obtained  by  combining  Eq.  (2.243)  with  Eqs.  (2.246)  and  (2.247).  It  is  again 
evident  from  these  results  that  a  zero  tt  pulse  should  exhibit  transparency. 


For  small,  h,  the  area  theorem  for  Q  branch  transitions  reduces  to  the  smaLL  0 
form  of  Eq.  (2.68)  when  orientational  averaging  in  included. 

As  in  Section  2.6,  pulse  chapes  may  be  obtained  in  the  limit  of  extreme  saturation 
breadening.  Setting  g(Au>)  =  f  ( Atu) ,  Eq.  (2.14^)  then  goes  over  to 


d2a 
dx  at 1 


a1 

2]TT 


lKmsin(Kmcr) 


^2.247) 


where  t'  =  t  -  x/c  and  y'  is  as  defined  in  Eq.  (2.30). 

/  a  Q\ 

Examples  of  steady  state  pulse  profiles  have  been  obtained  numerically'1-  .  For 
Q-branch  transition  with  j  =  2,  the  result  may  be  given  in  a  simple  closed  form 
One  obtains 
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a  =  -4tan-,(‘/5  csch  w/r) 


(2.248) 


where  w  =  t  -  x/v  and 


(2.249) 


The  electric  field  is 


IPEoC  _  4  -/5~  sech  W/r 

1)  “  ru  +  4  sech 2  w/r) 


(2.250) 


which  is  shown  in  Fig.  13. 
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Section  3 


CERENKOV  TYPE  RADIATION  FROM  LIGHT  PULSES 


3 • 1  Intr oduct ion 

Optical  rectification  in  non-linear  dielectric  media  produces  a  polarization  that 
follows  in  time  the  intensity  of  the  optical  pulse.  The  so-called  dc  polarization  has 
frequency  components  in  the  millimeter  and  submillimeter  region  of  the  spectrum,  and  in 
general  has  vector  components  both  longitudinal  and  transverse  to  the  direction  of 
light  propagation.  Submillimeter  radiation  can  therefore  be  generated  from  picosecond 
light  pulses.  A  general  theory  of  microwave  generation  via  the  dc  effect  does  not 
exist. 


The  most  important  feature  of  the  generation  of  the  microwave  field  from  its 
source  is  the  fact  that  the  source  polarization  moves  through  the  crystal  at  the  group 
velocity  of  the  light  pulse  ,  which  is  approximately  c/n,  whereas  the  phase  velocity 
of  the  free  microwave  field  in  the  same  crystal  is  usually  lower,  c/m,  where  n  and  m 
are  the  optical  and  microwave  indices  of  refraction.  Because  of  the  infrared  resonances 
between  the  optical  and  microwave  frequencies,  m  is  usually  larger  than  n,  so  the 
microwave  field  moves  slower  than  its  source,  in  close  analogy  with  the  Cerenkov 
radiation  from  a  relativistic  electron.  The  analogy  with  Cerenkov  radiation  will  be 
shown  to  be  very  close  for  the  microwave  field  generated  by  the  longitudinal  component 
of  the  polarization.  Our  understanding  of  the  microwave  generation  from  light  pulses 
in  non-linear  media  is  greatly  increased  by  studying  the  Cerenkov  analog  of  the  problem, 
viz.  the  ooupling  of  the  source  to  the  free  fields  of  the  dielectric.  For  example, 
one  learns  to  what  extent"phase-matching"  is  important.  A  comparison  of  the  induction 
fields  and  the  radiation  fields  shows  the  possibility  of  greatly  increasing  the  micro- 
wave  power  generation  by  changing  the  boundary  conditions  on  the  problem  in  such  a 
way  as  to  couple  to  the  induction  fields. 

This  work  deals  with  the  radiation  from  longitudinal  and  transverse  polarizations 
in  optically  isotropic  and  uniaxial  crystals.  The  physical  picture  that  one  associates 
with  these  phenomena  is  a  light  packet  accompanied  by  a  dc  polarization  radiating 
microwave  power  in  Cerenkov  cones,  one  cone  system  at  the  leading  edge  of  the  pulse 
and  another  at  the  trailing  edge.  The  interference  between  these  Cerenkov  cones  is 
completely  accounted  for  by  Fourier  analysing  the  polarization  pulse  into  its 
frequency  components.  The  radiated  fields  at  any  frequency  are  proportional  to  the 
corresponding  Fburier  component  amplitudes  of  the  equivalent  source  current.  In 
the  same  way,  the  interference  between  waves  originating  from  different  regions  in 
the  cross  section  is  exhibited  in  the  Fourier-Bessel  transforms  of  the  traverse 
distribution  of  source  polarization. 

The  first  problem  that  we  treat  is  that  of  radiation  in  an  optically  isotropic 
medium.  Both  longitudinal  and  transverse  polarizations  are  treated.  The  theoretical 
technique  is  to  solve  the  inhomogeneous  differential  equations  for  the  field  potentials, 
then  derive  the  actual  fields  by  differentiation  of  the  potentials.  The  second 
problem  is  that  of  a  uniaxial  crystal,  in  which  the  tensor  nature  of  the  dielectric 
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constant  so  much  complicates  the  problem  that  it  is  easier  to  deal  directly  with  the 
field  quantities  rather  than  the  potentials.  The  method  of  analysis  here  is  to 
resolve  the  source  polarization  into  a  continuous  spectrum  of  plane  waves,  find  the 
electric  field  components  for  the  ordinary  and  extraordinary  plane  waves  of  microwave 
field,  then  sum  them  all  through  integration  to  find  the  complete  ordinary  and 
extraordinary  wave  fields. 


3.2  Optically  Isotropic  Media 

The  microwave  field  is  derived  from  the  vector  and  scalar  potentials  A  and 


H  =  curl  A 


(3.1a) 


_  1  SA 

E  =  -  -  _=  -  grad 

c  at 


(3.1b) 


The  polarization  field  P  is  split  into  its  homogeneous  and  inhomogeneous  parts: 


(3.2) 


where 


E 


> 


(3.3) 


and  Pi,  the  inhomogeneous  part  is  the  source  polarization  caused  by  the  dc  effect  in 
the  non-linear  medium.  The  potential  field  equations  are 


2  m2  S2A 

V  A  -  — 

c  at 

C  dt 

-Hoj 

<1 

“e|” 

1 

~v> 

CM 

> 

= 

“ |  div  P. 

c  at 

m 

2 

div  A  +  — 

2i 

=  0. 

—  C 

dt 

(3A) 


(3.5) 


(3.6) 


This  last  is  the  gauge  condition.  Only  two  of  the  three  equations  are  independent. 
We  use  the  first  to  derive  A  from  Pj  ,  the  third  to  derive  i|r  from  A,  then  find  E  and  H 
by  differentiation  of  these  potentials.  The  source  of  A  is  a  current  density.  For 
Cerenkov  radiation  from  electrons  the  source  current  is  the  actual  current  due  to 
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motion  of  the  electron, 


and  the  right  hand  side  of  the  wave  equation  for  A  is 


-4lT.il 

c 


The  comparison  between  the  Cerenkov  radiation  from  an  electron  of  charge  e-^  and 
that  from  a  longitudinal  polarization  front  of  length  A  and  face  charge  Q  is  easily 
exhibited.  In  both  cases,  let  z  be  the  direction  of  the  motion  of  the  source.  For 
the  electron,  the  charge  density  is  given  by 


p  =  e15(x)8(y)  5(z  -  vt) 

where  v  is  the  electron  speed.  The  current  density  is 


(3-7) 


J  =  p v  =  e  v6(x) 6  (y )6  ( z  -  vt)  z 


(3.8) 


We  use  the  time  convention  exp  (-iajt)  throughout,  and  thus  have  the  Fourier  component 
of  current  at  frequency  m  as 


j  =  zj  =  ~  6(x)6(y)e1U)Z//v 

'll)  113  2tt 


(3.9) 


The  vector  potential  equation  is  therefore 


o  22  ei  •  / 

2.  m  id  ,  1  /  \  1U3Z/V 

7  A  +  — A  =  —  6  p  e  ' 

U3  q2  U)  TTCp 


(3.10) 


where  A  =  (0,  0,  A  )  and  where  §(x)6(y)  has  been  replaced  by  6(p)/2np  in  cylindrical 
coordinates.  The  polarization  front  is  represented  by 


P_  =  Q[|  -  1/tt  tan  1f(z  -  vt)/A}]  2  6 (x)6 (y) 


(3.Ha) 


Its  Fourier  component  equivalent  current  is 

dP 


.  _  a  Q  -A  si)  /v  iidz/v  6  ( p ) 
,  =  -  io)P,  =  z  —  e  11  e 

fit  /  — Iuj  2n  2rrp 

(13 


(3.Hb) 


Thus  the  vector  potential  ’.rave  equation  is 

2  2 


I  2  m  uj  \ 

V  +-r) 


,  4ttiu) 

A  =  - 

0)  c 


ttCo 


~I(13 

(3.12a) 

6(p)elu,z/v 

(3.12b) 
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where  again  A  =  (O.O.A  ).  The  energy  radiated  per  cm  of  path  per  unit  frequency 
—0)  (!) 

interval  by  tne  electron  is 


d2w 

dcudz 


Ji 

c 


( 1  — tt)  w 

\  m2v2/ 


(3.13) 


Hence,  by  comparison  of  the  two  vector  wave  equations,  the  radiation  from  the  advanc¬ 
ing  polarization  front  must  be  given  by 


For  a  single  electron,  e-^  =  ^.8  x  10  esu,  and  the  energy  radiated  per  cm  path  per 

cm--1-  bandwidth  is  ~5  X  10"28  erg,  For  a  picosecond  pulse  focussed  to  10"1-^watts/cm^ 

in  KDP  we  have  PT  =  8  esu  in  say  a  square  millimeter  section,  so  that  Q  =  8  x  10-^  esu 
51  11  -1 

=  2  x  10°e-^ .  The  energy  radiated  here  is  now  ~  2  x  10“  erg  per  cm  path  per  cm 

andwidth  at  wavelengths  long  compared  to  A  and  to  the  cross  section  dimensions. 

Admittedly,  the  energy  output  in  this  case  is  small,  but  techniques  for  increased 

coupling  can  be  found, and  the  important  point  is  that  the  equivalent  of  10®  relativistic 

electrons  in  a  cubic  millimeter  volume  is  easy  to  achieve  with  picosecond  pulses  and 

NLO.  The  defocussing  of  light  is  less  serious  a  problem  than  Coulomb  repulsion  and 

beam  collimation  in  megavolt  electronics. 


We  now  proceed  with  accurate  calculations  of  the  field  quantities  for  this 
Cerenkov  type  interaction  -  the  weak  coupling  situation. 


Calculation  of  the  Electromagnetic  Fields 


We  want  to  solve 


v2  ♦  ^  A*  •  P. 


\Z  iw 


(3.15) 


and 


„  .  iu/m2  ,  „ 

V#£W--c—  +»  -  0 


(3.16) 


for  the  general  polarization  pulse; 


Pilr,n  =  p  F,(x,y)F2(z  -  -^) 


(3.17) 


which  is  derived  from  a  collimated  picosecond  pulse.  Fourier  analysis  of  gives 


f2 ( z-vt)  =  [g2 (ou)einztu/c]  e-'"*  du 


(3.18) 
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The  following  quite  general  form  of  F  (x,  y),  expressed  in  cylindrical  coordinates, 
will  take  care  of  most  polarization  forms  derived  from  rectification  of  the  mode 
fields  of  a  laser  output: 


F,U.y)-F, (/>,<*>)  «  £F i(p)e't* 


(3.19) 


The  vector  p  will  have  longitudinal  and  transverse  components.  Because  of  the 
linearity  of  the  problem,  we  may  consider  either  one  component  alone  and  also 
select  one  1-component  of  the  transverse  polarization  distribution.  Equation  (3-15) 
now  becomes 


(»*  ♦ 


Acu  s 


47TIOJ 


Ft{p)e't4>  G2(<u)einwZ/c 


(3.15a) 


where  A  stands  for  either  of  the  rectangular  components  A  •  X  or  A  •  2  correspond¬ 
ing  to  transverse  or  longitudinal  polarizations  respectively.  The  "otvious  method 
of  solving  Eq.  (3.15a)  is  to  substitute 


AW*  u (p)e't*  einwz/c 

which  yields  for  the  scalar  u(p)  the  inhomogeneous  Bessel  equation 


(3.20) 


B'USE  -p  ippw*  (s*-$)u 


(3.21) 


where  S2  =  (m2  -  n2)  is  the  square  of  the  transverse  k-vector. 


The  solution  u  is  readily  found  and  will  be  written  down  later.  For  the  moment 
it  is  instructive  to  look  at  the  structure  of  the  solution  being  forced  on  this  problem. 
The  only  boundary  condition  is  the  radiation  condition.  Equation  (3* 20)  is  consistent 
with  a  linear  loss  (diffraction)  balancing  the  energy  generation.  The  exponential 
factor  einujz/c  in  Eq.  (3-20)  is  an  expression  of  the  translational  invariance  along 
z  of  the  problem.  There  is  no  possibility  of  exponential  growth  of  the  microwave 
field  along  z  here.  The  magnitude  of  F  (p)  in  Eq.  (3.15a)  is  independent  of  z,  which 
is  an  approximation  that  neglects  diffraction  of  the  picosecond  light  pulse  along  the 
path  of  interest,  and  neglects  absorption  of  the  light  pulse.  Of  course,  the  micro- 
waves  rob  the  light  pulse  of  some  of  its  energy,  and  they  take  their  momentum  from 
the  crystal  directly,  and  some  momentum  from  the  light  pulse  indirectly  via  the 
crystal . 
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Fields  Generated  by  Longitudinal  Polarization 

When  p  =  2,  the  solution  of^Eqs .  ^  3  •  ^0)  and  (3«2l)  is  A  =  (0,0,uei^eirUUZ^C^ 
whereupon  Eq.  (3.16)  gives  i|r  =  —  ue  ™e  .  From  Eq.  ”T?.l)  we  therefore  obtain 

W  m 


(3.22a) 

The  energy  radiated  is  found  by  integration  of  the  Poynting  vector 

(3.22b) 

S'*'*  5V-5-Re{E"»xH'"} 

(3.23) 

JZ& (,ld* *  fRe{/I(^x3<*')d"} 

(3.24) 

Taking  into  account  the  frequencies  m  and  -u)  simultaneously  we  calculate  the  energy 
radiated  per  unit  frequency  interval  through  a  cylinder  of  1  cm  length  coaxial  with 

2  by  integrating  the  co-derivative  of  Eq.  (3.24)  over  the  cylinder  surface.  The 

result  is 


dwdz  '  fo  4  Re{(EwXHw)/>  +  (E-wx  H-w)/»}/>d<^  (3.25) 

When  p  is  large  enough  for  the  cylinder  to  enclose  the  sources  completely,  Eq.  (3.25) 
is  independent  of  p.  From  Eq.  (3.22), 


du 


dP 


(3.2 6) 


and  Eq.  (3.25)  becomes 


_d 

d 


l2w  .  oi  /,  n2  \  f  du  n  du 

TJdT  ’  Tl,_7^)27ri^lud^  UW, 


(3.27) 


This  formula  is  for  F  (x,  y)  ~  e1  ^  =  cos  ftp  +  i  sin  ?.cp  wherein  the  cosine  and  sine 

terms  radiate  independently.  If  only  one  is  present  (e.g.  cos  j£cp)  we  have  to  divide 

Eq«(3.27)by  2.  Note  that  i[2tt  2  i  ,  i^-cp  -iLcp, ,  except  when  L  =  0. 

2^  J0  cos  f  cpdcp  =  f  Re[e  T  e  7 
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Fields  Generated  by  Transverse  Polarization 

When  p  =  x,  the  solution  of  Eqs .  (3*20)  and  (3*21)  is 

Ao,  =  (uelJ*elfWZ/c  ,o  ,o)  —  (u  cos  ^einwz/c,o ,o) 
s  u  cos  i<£cos<£  eina*/c /8 -u  cosi<£  sin^e'^2^  (3.28) 

where  the  cp-dependence  has  been  written  more  explicitly.  Equation  (3.l6)  gives 

^  [djr  “s^+cos^  +  f-  Sin^^sin^]  e,nw!/c 

Use  of  Eqs.  (3-1)  and  (3-2)  gives  the  fields  E  and  H  from  which  one  readily  obtains 

—01  —01 

Re  {(EwXHto)^  -L-^cos^.j  |u|S--u^_)(sin^+  -^-cos2^,)  (3>3o) 

jL 

We  have  dropped  terms  in  E  and  H  that  fall  off  faster  than  l/p2  because  they  will  not 
contribute  to  the  radiation  field.  These  induction  fields  are  capable  of  coupling 
energy  out  of  the  system  through  microwave  structures  that  come  in  to  within  a  couple 
of  wavelengths  of  the  source.  The  presence  of  such  a  structure  changes  the  boundary 
conditions  on  the  problem,  so  that  our  equations  no  longer  give  the  correct  expressions 
for  these  fields.  A  simple  analogy  may  be  helpful  here.  A  radio  frequency  choke 
is  a  poor  radiator  because  of  the  small  overlap  of  its  fields  on  the  free  fields  of 
the  vacuum.  The  same  choke  coil  when  UE'vi  as  the  primary  winding  of  a  transformer 
is  capable  of  coupling  strongly  to  a  transmission  line  connected  to  a  secondary 
winding. 

The  formula  for  the  energy  radiated  by  the  transverse  polarization  as  calculated 
from  Eq.  (3-25)  depends  on  whether  £  =  0,  C  =  1,  or  £  >  1.  The  result  is: 


d2W 

dwdz 


i  =  o 


d2w 

dwdz 


d2w 

doudz 


■) 

(3.31a) 

du  \ 

dp  1 

(3.31b) 

du  \ 

dp  ) 

(3.31c) 

The  Eqs.  (3.27) (longitudinal  polarization)  and  (3. 31) (transverse  polarization) 
both  contain  the  expression  .  ,  du  -  du.  in  which  u  is  the  function  that  solves 
Eq.  (3.21)  1<’(uarUd^ 
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.  du  -  duv 
Evaluation  of  ip(u  —  -  u  — ) 

dp  dp 


The  solution  of  Eq.  ( 3 * 2l)  can  be  found  most  readily  from  expression 


u ip)  --  -  l<1Jc2(a;)  So  <*l{P/Po)n{Po)fh*Po 

where  g  (pip  )  is  the  Green's  function  of  the  problem 
e  1  o 


(3-32) 


Bj  w/p, )  =  - 

The  Green's  function  for  Eq.  ( 3 •  33)  is 

qJt(p/p0)2^\  Jj,(sp)Hj,($Po)  P<P0 

and 

q/P/Po)2rr*  i J^(sA,)H2(sp)  P>P0 


(3-33) 


(3- 34a) 


(3.34b) 


Both  expressions  (3*34)  must  be  used  to  express  the  field  inside  the  source.  We  need 
only  Eq.  (3* 34b)  for  the  radiation  fields.  Combination  of  Eqs.  (3-32)  and  (3.34b) 
gives 


u(/»)  '  2 " * “  ,W 1  [  fopt (ftj  J„( Sf,„ ) p0 d/>„]  ( sp) 


27T2Cl)  G2  (CD) 
C 


G*(S)  H^S p) 


(3.35) 


where  we  use  G^(s)  to  denote  the  Fourier- Bessel  transform  (the  expression  in  square 
brackets)  of  the  source  distribution. 

Thus 


'/‘[‘lip -*%-)-  4’r<#lGH'  h'/Wi-c.c.] 

G2(0l»)  |2  [gA(S  )]  2 


5  I67T3-^- 


(3.36) 


Energy  Radiated  in  an  Isotropic  Medium 


Equation  (3*36)  is  now  substituted  into  Eqs.  (3-27)  and  (3.31)  with  the  follow¬ 
ing  results 
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Longitudinal  Polarization: 


where 


and 


G2(OJ)|2|V(S)]2€* 

=  1  for  1  =  0 

=  \  for  i=  1,  2,  3,  .  .  .  . 


crw 

dzdu 


:  87T4 


os’ 

“c* 


-5 


Transverse  Polarization: 


d2  W 
dzdcu 


l- 1 


G2(cj)|2[V(S)]Z 


(3.37) 


(3.37a) 


(3.37b) 


(3.38a) 


d2w 

dzdcu 


/  9*  I 


47T4 


Cl)3 

C2 


■('♦&) 


|GZ(W) 


(3.38b) 


In  the  above  expressions  is  the  transverse  component  of  the  k- vector  of  the  field. 


For  any  actual  polarization  distribution,  one  uses  the  decomposition  of  Eq. 
(3.19)2and  sums  the  Eqs.  (3*37)>  (3*38)  over  f.  Equation  (3-37)  contains  the  factor 
(1  -  £^)  in  close  analogy  with  Cerenkov  radiation.  The  transverse  polarization 
case  cfoes  not  have  this  factor,  because  the  radiated  power  does  not  vanish  when 
m  =  n.  This  is  in  sharp  distinction  to  the  case  of  Cerenkov  radiation  from  a  trans¬ 
verse  electric  dipole  on  a  relativistic  particle.  The  distinction  arises  because 
for  the  particle  the  equivalent  current  sources  are  always  along  the  direction  of 
motion,  and  furthermore,  the  transverse  dipole  partly  transforms  to  a  magnetic  dipole 
at  right  angles  to  both  the  electric  dipole  and  the  direction  of  motion.  Such  a 
Lorentz  transformation  does  not  enter  the  NLO  case,  where  the  equivalent  current 
source  may  actually  be  transverse  to  the  direction  of  motion. 


The  case  m  =  n  represents  a  match  between  the  speeds  of  the  optical  pulse  and 
the  microwaves.  In  the  case  of  SHG  this  phase  matching  causes  a  large  increase  in 
power.  The  present  situation  is  different,  because  of  the  essential  role  played 
by  diffraction  of  the  microwaves.  Since  (l  ±  n^/m  )  is  slowly  varying  in  n/m, 
nothing  spectacular  happens  at  or  near  n  =  m.  Transverse  and  longitudinal  polariza¬ 
tions  radiate  about  equally  well  according  to  our  formulae. 


Next  we  examine  the  case  of  an  anisotropic  medium  to  see  what  new  features  are 
introduced. 
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3-3  Optically  Anisotropic  Media 

This  method  of  analysis  is  more  powerful  than  the  one  used  for  isotropic  media. 
In  essence,  we  resolve  the  polarization  pulse  into  its  plane  wave  components,  find 
the  electromagnetic  plane  wave  for  each  polarization  wave  component,  then  add  all 
the  resulting  plane  waves  to  obtain  the  total  field. 


Again  we  start  with 


vt. 


Pj  (£,t)  =  pF1  (x,y)  F2  (z  -  — ) 


(3-17) 


The  Fourier  transform  of  F^  is  given  by 


F,(«.»!  =//"Gl(k.,ky)eilk“*k'»,dk. 


dk» 


(3.39a) 


or 


F,(^)=^*/02’de.G,(x,e)ei'",c'>s8-+' 


(3.39b) 


where  x  =  p  cos  cp,  y  =  p  sin  cp,  k  =  K  cos  9,  k  =  K  sin  0.  The  inverse  of  Eq.  (3*39b) 
gives  ^ 

(J  40a) 

Again  assuming  that  F^(p,cp)  =  EjF^(p)e^"^ 

1  (3.40b) 


J — rreit(8-w/2  )  rl 


2v  Z 

t 


G  *(K) 


(3- toe) 


where  G  (x)  is  the  Fourier-Bessel  transform  defined  previously  in  Eq.  (3-35 )>  and 


/•»  /*2  7T 

F|  (p,4>)--  z  Jo  dKKj0  68  ~  e'*ie~v,z)  Gt{K)e'*pcosie-<t,)*ilt 


(3.39c) 


Using  Eq.  (3«l8)  to  transform  F2  now  gives  the  combined  result 


,W)r&p~2V" /dwG2(w)j[d»cKGi(K)/d0{eiw(,-ni/#  d) 


(3.4i) 
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The  quantity  in  curly  brackets  defines  the  frequency,  wavelength,  direction  and 

starting  phase  of  an  elementary  plane  wave  component  of  polarization  with  amplitude 

P  .do), Kdx, do  where 
(uK 


Po,k‘  jLre-i^GjIuOG'l*)  (3.1,10) 

If  such  a  plane  wave  gives  rise  to  an^eligctromagne^c  wave, with  one  £^.eld  component 
(say  x)  of  amplitude  and  phase  g^  e  e  D  '  ™  dutfCdKdQ, 

the  total  field  component  due  to  tne  totality  of  elementary  plane  waves  is  given  by 
the  integral 

£,<£.')  ‘  %  f&r11*  *  I*"68  GaUG<l*>  (3.42) 

Our  problem  reduces  to  finding  the  ratios  g  /P  ,  g  /p  and  g  /p  (where  e.g. 

x  uK  y  mK  z  uK 

g  s  g  )  and  integrating  Eqs.  (3.42)  to  obtain  the  electromagnetic  field.  These 
ratios  are  found  in  the  general  case  by  solving  the  inhomogeneous  wave  equation. 

The  general  solutions  can  then  be  specialized  to  the  required  uniaxial  crystal  or 
biaxial  crystal,  etc. 


Solution  of  the  Inhomogeneous  Wave  Equation 

In  the  anisotropic  medium,  the  displacement  vector  is  related  to  the  field 
intensity  and  the  inhomogeneous  part  of  the  polarization  by 


D  -*  €  •  E  +  4  7T  PT 


The  wave  equation  is  now 


I  a*E  -4tt  d*Pl 
VkV  K  E  +  -ZT  €  ■  -r-4-  *  ■— T-  . 

~  c‘  at  cz  a  tz 


Resolved  into  Fourier  components.  Equation  (3.44a)  becomes 

z*s\{ zrf  k2{e  -$($•«)}-  €  $|  s  4fpJ 


(3.^3) 


(3.44a) 


(3.44b) 


~k 


where  g,  e  =  E  ,  the  Fourier  transform  of  E,  and  k  =  ks  •  A  subsidiary  condition, 
coming  fron  div  D,  =0  is 

0--  div  D«  *  k«(€»Ek+  4  tt^) 

•  A  A  .  _  _  A  A  /  p  (3.45) 

•••s»€.e!-4irpw(fp.s/^K 
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(84) 

The  inhomogeneous  wave  equation  (3.44b)  is  solved  by  the  method  used  by  Kleinman' 
in  his  study  on  SHG. 


E,.  S2E2.  SC, 


(3.1-6) 


where  S  is  the  wave  normal  (k  =  kS) ,  e  is  the  direction  of  polarization  of  the  E 
field  for  a  free  ordinary  wave  and  e^  is  the  direction  of  polarization  for  the  free 
extraordinary  wave,  e  ,  e ,  and  S  are  not  orthogonal.  When  Eq.  (3.44b)  is  resolved 
along  these  three  directions  we  obtain 

(^f{-k2S(S«,)E,-k%«^E2.k^E,.KS2E2} 


-  (€•#,  E,  +  €  •  ^2E2)  -€«tE3=  A7T  (3.44c) 

The  quantities  e  *  e  and  c  *  e2  in  Eq.  (3.44c)  are  found  by  solving  the  homogeneous 
wave  equations  for  the  ordinary  and  extraordinary  waves  respectively.  These  equations 
are 


Substituting  e  •  e  into  Eq.  (3.44c)  we  find 

(i-)2  [(kf-kz)  {m4,)  -«,}£,♦  (kz-kz)  {S(s.ez)  -^}eJ 

The  scalar  product  of  Eq.  (3-44d)  with  S  gives 


(3.47) 


(3.44d) 


E3  -  -  47rp  ($•§)/ (s»e»§)  (3.48a) 

3  UIK 

which  is  now  substituted  into  Eq.  (3*44d)  to  give 

(s)!(kf-k!)  {»(*•«,) -«,}e,  *  (c-flk^-k2)  {StS.Sjl-it.jE, 

=  4Trpw^[^  -  (e  *s)  (p.s)  /  (s#e.  §)J  (3. He) 

In  anisotropic  media  the  D  vectors  of  the  ordinary  and  extraordinary  waves  are  polar¬ 
ized  in  the  directions  5  and  9.  respectively,  which  are  both  perpendicular  to  the 
wave  normal  S,  and  each  is  perpendicular  to  the  non- corresponding  e  vector 


A  A  A  A 

d|  »e2=  o  =  d8*e, 


(3.49) 


The 


scalar  products  of  Eq.  (3-44e)  with  d  and  d^  yield  the  remaining  field  components: 
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E(=47r(^f-)  Pw<^i#P  “  (^i*€*s )($•§)/ (s»«*§)]  /(8|«e,)(k2-k2,) 

E2=  47T'(t')  Pu>«  [$2*P  _(^2*€*sHp*s)/(S«€«s)/(^2  *e2)(k2-k2  ) 

In  equations  (3*47)  and  (3.48),  k  ,  denote  the  wave  vectors  for  ordinary  and 
extraordinary  waves  and  are  given  in  terms  of  the  phase  velocities  Vq  and  v  of  the 
ordinary  and  extraordinary  waves  belonging  to  the  wave  normal  S;  i.e.,  k  =  oj/v  , 

k2 - »At. 

These  ratios  like  £  /p  needed  for  Eq.  (3.42)  come  from  Eqs.  (3.46)  and  (3*48) 
when  e  ,  e  ,  and  S  are  projected  on  to  the  x,  y,  and  z  axes.  To  progress  from 
Eqs.  (3.48;,  we  specialize  to  the  case  of  uniaxial  crystals. 


(3.48b) 

(3.48c) 


Uniaxial  Crystals 

Before  we  can  use  Eq.  (3.42),  we  have  to  evaluate  the  various  scalar  product 

terms  that  appear  in  Eq.  (3-48).  Assume  that  the  motion  of  the  source  polarization 

is  in  the  direction  of  a  principal  axis  of  the  dielectric  tensor  e.  The  direction 

of  motion  is  z  •  x  and  y  are  also  principal  axis  directions.  For  a  uniaxial 

crystal,  e  =  e  makes  2  the  optic  axis  while  e  =  e  makes  x  the  optic  axis.  In 
xy  yz 

the  first  case  the  source  moves  along  the  optic  axis  with  either  transverse  or 

longitudinal  polarization.  The  latter  of  these  two  possibilities  can  be  solved 
almost  by  inspection  by  looking  at  the  solution  for  Cerenkov  radiation  from  an  electron 
moving  along  the  optic  axis  and  adjusting  the  coefficients  in  a  way  suggested  by 
comparison  of  Eqs.  (3.13)  and  (3-37).  We  will  concentrate  on  the  transverse  polariza¬ 
tion  case,  for  motion  along  the  optic  axis.  This  case  is  easier  to  treat,  though 
of  perhaps  less  interest  than  the  case  where  the  motion  is  perpendicular  to  the  optic 
axis  with  the  polarization  along  the  optic  axis.  In  this  latter  case  the  Cerenkov 
cones  are  not  circular  and  the  integration  over  0  in  Eq.  (3*42)  is  best  omitted  so 
that  one  can  consider  the  azimuthal  distribution  of  radiated  energy. 


Figures  14  and  15 depict  the  ellipsoids  of  wave  normals  for  uniaxial  crystals. 


The  relevent  phase  velocities  v  corresponding  to  the  wave  normal  S  =  (S  , 
are  derived  from  the  equation  fSr  the  normal  surface 


Sz) 


+ 


*  0 


(3.50) 


For  the  source  moving  along  the  optic  axis  (Figure  l),  v  =  v 
vz  =  vg  =  c//sz.  The  solutions  for  v  are  then:  J 

Vp  =  V0  =  c/v^T 


c//«x  and 

(ordinary) 

(3.51a) 
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Vp=  Vc 


=  (V  cos2/3  +  -^-sin2/3) 


(extraordinary) 

(3.51b) 


Equations  (3*51a)  and  (3.51b)  correspond  to  a  sphere  and  an  ovaloid  respectively. 

Here  0  is  the  angle  between  S  and  z.  For  the  source  moving  transverse  to  the  optic 

axis  (Figure  15).  v  =  v  =  cZ/’e  .  and  v  =v  =v  =c //$  .  The  phase  velocities  are 
x  l  x  yzO  z 

Vp  =v0  =c/y/7^  (ordinary) 

(3.52a) 


Vp  =  Ve=^-^“  sin2/? COS20  +  -^-  (sin2/3sin20-tCOS2/3)j  (extraordinary) 


(3.52b) 


Here  0  is  the  azimuthal  angle  for  the  wave  normal  S. 

We  now  have  sufficient  information  to  evaluate  the  various  unit  vectors  and 
the  scalar  products  in  Eqs.  (3.^8).  The  results  for  the  source  moving  along  the 
optic  are  as  follows: 

=  (sin  0,  -cos  0,  0) 

3.^  =  (cos  8  cos  0,  cos  0  sin  0,  -sin  0) 

S  =  (sin  0  cos  0,  sin  0  sin  0,  cos  0) 

e-S  =  e  (T|  sin  0  cos  e,  7]  sin  0  sin0,  cos  0);  7]  =  e  /e 

Z  X  z 

^1  =  ^1  2  2  2  — 

e  =  (cos  0  +  T)  sin  0)2  (cos  0  cos  0,  cos  0  sin  0,-T]  sin  0) 

P  =  * 

cL*p  =  sin  0;  >p  =  cos  0  cos  0;  $  •  S  =  sin  8  cos  0;  (d^  ^  e 1)  -  1 
d.  -  (cos2  0  +  Tj2  sin2  0)"2  (cos2  0  +  7)  sin2  0);(d  •  e  •  S)  =  0 

^  *  a  a  2X  2 

(cl  .e-S)=  e  (Tj  -  l)  sin  0  cos  0;  (S  •  e  •  S)  =  (e  cos  0  +  7]  sin  0). 

(3.53) 

rT1he  results  for  the  source  moving  transverse  to  the  optic  axis  with  the 
polarization  along  the  optic  axis  are 


d1,  =(i-sin2/3cos20rl/2  (0,  -  cos/3, sin/9sin@) 

d2=  (i-  sin2/9cos2@rl/2  (i-sin2/3cos20,  -sin2/3sin0cos0,-sin/3cos/3cos0) 


s  =  (sin/3cos0,  sin/3sin0,cos/3) 


A  4 

e,  =  d. 


=  (i-sin2/9cos2©rl/?  [i  +  (i72-i)sin2/9cos20]’l/2  (i-sin2/3cos20,-i7sin2/?sin0cos^-^sin/3cos/3cos0) 
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€  *  s  =€z  (17  sin  /9cos0,  sin/3sin0, cos/3) 

A  A 

P  =x 

d, •  p  =  0 ;  d2*p  =(1- sin2/3cos20)l/2  ;  p*  s  =  sin/3cos0,  d,  •  e(  =  0  ; 

d2*  e2  =  [|  +  (i72-0  sin2/Sc°s26l]  1/2  [i  +  (i?-i)sin2/3cos20j ;  (d‘|*€*s)  =  0; 

(d2‘€*s)  =  (17-1)  sin/3cos0(i  - sin2/3cos20)l/2  ; 

(s*f*s)=  ez+€z  (17 -i)sin2£cos20  (3.54) 

We  are  now  in  a  position  to  combine  Eqs.  (3*46)  and  (3-48).  The  case  of  motion 
along  the  optic  axis  is  treated  here.  The  following  substitutions  are  first  made: 


k2  =  *2+n2($2 
k2-k,2=(^)2(T2+n2-*?«z) 
sin/3  =  r /(r2+n  2  }l/2 


2  r2 


cos)S=n  /  (r2+n2)l/2 


(3.55) 


The  field  components  to  be  substituted  into  Eq.  (3*42)  are  thus 


£x/47Tp( 


sin 

‘  r2+ni-‘ 


Z9 


n2cos  z9 


j?cq£l 


Tz+nz-77tz  (77r2+n2)(T7T2+n2-T7€z)  €z(77r2+n2) 


a  /A „ _  _  -sinflcosfl  .  n2cqsflsinfl _  ,  r2cos0sinfl 

V4  Pojk-  r2+n2-T7€z  (i7T2+n2)(i7T2+n2-77€z)  €z (i7T2+n2) 


-9  ,A  _ _ -qnrcosfl 

^Z/47rP  "  .  _2\/ _ 


nr  cos 0 


{77^2+n2)(i7r2+n2-'i7€z)  ^(^z+nZ) 

The  integration  over  9  is  performed  first.  These  integrals  are  of  the  form 

^  da 

2  2 

where  C(e)is  the  sin  9,  cos  9,  sin  9  cos  9,  or  cos  9.  These  all  integrate  into 
derivatives  of  Bessel  functions.  For  example,  consider  the  easiest  one  with 
C(9 )  =  cos  9-  We  substitute  p  cos  9  =  x  and  p  sin  9  =  y.  This  integral  is  now 


(3.56a) 


(3.56b) 


(3.56c) 


±  Jocose  e 


de 
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Using  the  well  known  integral  representation  of  the  Bessel  function,  J  (£)  = 

-l  .  ^ 

■5 —  <j)  e^  C0S  ^  +  in^d£p,  the  above  integral  is  obviously  a  derivative  with  respect 

to  x  of  a  Bessel  function,  and  in  fact  is  -  e^  -  J'  (xp).  The  other  integrals  are 
only  slightly  more  complicated.  ^ 

The  Eqs.  (3.4-2)  for  E  ,  E  and  E  in  terms  of  the  remaining  integrals  are  now 
as  given  below  in  Eqs.  (3-fj)  ^  7 

E„  =  l4ireifTdo,  f-""-n:,Cl(tfGzM/VrG'ff){ 

- (v-TKn2^;^—-)  l^*3^  +  COS^J.Wcl] 

[Sin2*  +  COS*J,Wc)]} 


(3.57a) 


fz(VT~S~+nZ)  L  T’  (WTp/ c) 

(£)%,(«,£“< It  r  0W){ 

{  [  cos</.sin^.Ji*“r^  -  cos^sin*  tfW/rtl] 


(t) 


Tqr„-^i+^.^j  [  °os  ^  sinf  -  cos  ^  sin  ^  If“(iu  rp  /c)]  <3-57b) 

[cos<t  sin  0  ^  'ojrP/c*  _  cos^  sir^^'twrp/fc)]  J 


£z(^r2+nZ)  l (u»r/wfc) 

+  ..(,Xn*)  [-*  co*^^'  CunpA*]} 


{  (T,T*+n»)(T,T»+n»-i,.p  [ 1  C0S *  *'  ‘‘"H 


(3.57c) 


We  now  proceed  with  the  integrals  over  t«  Quite  a  few  of  the  integrals  can  be 
disposed  of  immediately.  Those  /that  involve  a  Bessel  function  divided  by  p  have 
asymptotic  forms  that  go  as  p  '  and  contribute  nothing  to  the  radiation  field. 

They  are  induction  field  terms.  Consider  now  an  integral  like 

jf ”d  r  r  G‘  (^)  •  ■  J("  ($  TP) 

Set  J/  (*cf  TP' :  y  H£ 1  (^f  Tp) +  "2  h£  (~  tp)  .  Replace  the  Integral 

along  the  positive  real  axis  by  the  sum  of  the  two  contour  integrals  shown  in  Figure  16 
where  „(i)  n  is  integrated  along  the  upper  contour  and  h(^)H  along  the  lower. 


,(1) 


;(2) 


H  and  H  vanish  along  their  respective  contours  at  infinity.  The  vertical 
contributions  on  the  two  contours  cancel  as  can  be  shown  by  using  the  relation 
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iH^  (ix)  =  (-l),f+1  (-ix).  Thus  the  sum  of  the  contour  integrals  equals  the 

*1/  K* 

wanted  integral.  The  poles  are  located  at  t  =  V"Tl(e  +  1a)  -n2  .  These  are  in  the 
first  and  third  quadrants.  The  lower  contour  integral  encloses  no  poles,  thus  it 
vanishes .  The  upper  contour  integral  is  2rri  times  the  residue  at  t  -  T)e  “  and 

thus  2T'G<[  T  (7)€z~r)2^IZ .  ni'*  {  i:  ^V€z  ~n2)l/2^}  .  The  same  technique  is 

used  for  the  other  integrals.  The  contributions  from  the  pole  at  t  =  -i(n  +  iA ' ) //T| 
involve  the  function  H'  ^(-i  ^  -ip)  which  falls  off  exponentially  in  p,  and  gives 

no  radiation  field.  The  surviving  integrals  are  now: 

E„  =-2,r2ipif*£dw«''i“'('"'’I/Cl  Iff  G2(u)  [G«(kT0)H;""(kr°/>)  Sin2*#. 

+  G2(k*)H'""(kTV)COS2*] 

Ey  =  -  2  ir2  i  £  e'!*f  due’1"11  ni/cl  (M-)2  g2(uj]  [  G*(k°)  n111'  (kt0/5)cos*#>  Sin<(i 


Gf(kTe)  h)i1// (kjp)  cos<£sin<#>] 


E7  =  -27r2i  Z  e...  (3-58) 

0  e  1 

Here  and  k~  are  the  transverse  k  vector  components  for  ordinary  and  extraordinary 
waves 


k°  ■  f  (V,  -n2)1'2  k,*  *  ?  («:- 


(3.59) 


To  justify  the  labelling  in  Eq.  (3*59)  we  compute  the  corresponding  phase  velocities. 
For  the  ordinary  wave  we  have 


<,*•(£?  ="2(?f  +  <KT°)2  -  iff’, 


(3.60) 


Therefore  vQ  =  which  is  correct  for  the  ordinary  wave.  For  the  extraordinary 

wave  we  have 


*<?  ■  (£)*■  n2(f  f*  (Kr*  )2  -  (ft  [«,  *n2  (,- A.)] 


can  be  solved  for  n  . 


N*  -E  C0S*/3» 

i.(iJ^)cos2/9« 

Substituting  Eq.  (3.62)  into  Eq.  (j.6l)  we  get 


(3.61) 


A  _  _ 2/3  n 

Let  B  be  the  angle  between  k  and  2.  Then  COS  — —  -  /. - rr  ,  which  can 

e  ~e  *2  4,02  V  ~  T)l 


(3.6 2) 


(3.63a) 
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or 


ve2  _  cos2/3e  sin2/3e 
c2  *x 


(3.63b) 


This  is  identical  to  Eq.  (3. 51b),  so  our  identification  of  this  wave  vector  with  the 
extraordinary  wave  is  correct. 


The  Eqs .  ( 3 • 58)  are  of  the  form 

Ex  =  Lx{  E x(0’(w)  +E*Cel(u)}  e"iwtdoj 


(3.64) 


with  similar  expressions  for  E  and  E  .  The  superscripts  o  and  e  refer  to  the  ordinary 
and  extraordinary  fields  and  tXe  tvm  fields  can  be  treated  separately.  The 
asymptotic  expressions  for  E  ,  E  ,  and  E  can  be  written  down  by  using  the  asymptotic 
relations:  x  y  2 


-  HptlV/(z)  a  -  i »!%) »  Hf(l,(z) *  e  n/*] 


(3.65) 


The  results  are: 


*•%)*!  2  7T2i(cu/c)2G2UGf(kT0)|-^j^-)  c  ' (hT  ^ J  sin2 


(0) 

y 

(0) 

3 


(3.66a) 


«•»  r*  2  n2  /u>\2  t,  (  2  V,Z  •)»>>•) 

4"  <“>  •  p  (f )  <*«■»  «'<*?>  (7^)  * 

T- .  f  ^  ,ff  (^f 


cos*sm*  e'‘* 


•/* 


F**‘(wl  *  £  -2tt2. 
( 


If) 

i(*  H*nv/C 

COS  ♦  6 


W*  °2W  e'<K? >  l^)' 


(3.67a) 


The  cylindrical  components  arc  readily  found  from  the  rectangular  components.  The 
angular  factor:?  are  as  follows* 


Ep°\u;)  •*»  s«n$  Sin2$ 
E^’(w)~sin*C0S2*  #'** 


(3.66b) 
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e pe\w)~  cos<£  e'1* 

Ef\u>)  =  0 

The  magnetic  field  components  are  obtained  from 


(2.67b) 


H  M  =  -i  (£)  curl  E  (w) 


(3.68) 


Hp0)(u>)  = 


■n  e£*M 


»  E{p0\u>)  Hj°W) 


Ck*  to) 

U)~ 


(w)  (3.69) 


H^*\cu)  *  0  »n  E  *lw)  -  E2 


Hz(#)(w)  «  0 


(3.70) 


Evaluation  of  the  Poynting  vectors  gives 


[E(®M  X  H(®m]  =  E  ^E^0’*  8 ir*($f  IG^2  [G*(kT°)]2  y  Sin^COS2 2^(3- 71a) 


and 


(3.7ib) 


Integration  of  the  Poynting  vector  gives  the  energy  radiated.  For  the  ordinary 
waves  ve  find 


(3.72) 


where  tf  ■  1  for  1*0 

•  1/8  for  cos  3/6  for  sin  tp 
■  $  for  <1  >  1, 

while  for  the  extraordinary  waves  the  result  is 

5^=4**  $  mm 

3.4  Susaary 

The  interesting  case  that  we  still  have  not  worked  out  in  detail  is  that  where 


**1 


(3.73) 
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the  source  motion  is  transverse  to  the  optic  axis,  with  the  polarization  along  the 
optic  axis.  This  will  probably  be  examined  later.  The  case  of  motion  along  the 
optic  axis  with  longitudinal  polarization  is  less  interesting,  but  the  problem 
can  be  solved  by  inspection  as  noted  earlier. 

The  coupling  problem  solved  here  between  the  NLO  source  polarization  and  a 
microwave  field  is  probably  the  only  one  that  has  been  treated  anywhere  to  date.  We 
were  able  to  treat  the  problem  fairly  exactly  because  the  microwave  field  we  considered 
is  unguided.  The  case  with  different  boundary  conditions  on  the  microwave  field 
will  be  more  relevant  to  the  submillimeter  wave  generation  problem,  and  we  feel 
better  prepared  to  look  at  such  cases  now  that  we  understand  the  most  elementary 
one. 
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Section  k 

AN  EXTENSION  OF  TPF  TECHNIQUES  TO  SIMULTANEOUS  POWER  SPECTRUM  AND  PHASE 

STRUCTURE  MEASUREMENTS 


There  Is  considerable  interest  in  non-linear  optical  phenomena  that  modify  the 
amplitude  and  phase  structure  of  short  optical  pulses.  Some  of  these  phenomena  have 
been  discussed  extensively  in  the  literature  already,  such  as  the  self -steepening 
of  light  pulsea^nropagating  in  Kerr-active  liquids  and  the  associated  self -phase 
modulation* 05 '  '  Indirect  evidence  of  the  phase -modulation  has  been  inferred  from 

spectral  broadening  measurements}  '/  which  shew  a  characteristic  spacing'  between 
the  zeros  of  the  power  spectrum  of  the  output  pulses.  More  direct  evidence  has  been 
obtained  by  studying  the  output  from  a  grating  pair  compressor  with  and  without  the 
Kerr-activu  llauld  in  th*  b«»*m.(°9) 


Evidence  of  self -steepening  processes  occurring 


in  a  node-locked  laser  has  been  obtained  previously  in  this  laboratory. 


(90) 


More 


recently,  interest  has  developed  In  non-linear  modifications  to  the  propagating 


pulse  in  the  presenc 

rectification,^1)  a 

power  optical  pulses. 


nee  of  the  propagating  microwave  signal  generated  by  optical 
in  almultaneous  front  and  back  end  shock  formation  in  high 


Much  of  this  Interest  centers  arouqd  the  possibility  of  geneiating  optical  pulses 
of  ~10’1^  second  or  shorter  duratlon.(^)  The  various  proposals  for  10"^  sec  pulses 
involve  the  addition  of  phase  modulation  to  the  pulse  and  the  subsequent  compression 
of  those  parts  of  the  pulse  that  have  the  appropriate  frequency  sweep  characteristics, 
using  pulse  compression  techniques  similar  to  those  initiated  in  this  laboratory. ( ^' 

In  this  type  of  investigation,  it  would  be  desirable  to  obtain  simultaneous  power 
spectrum,  phase  structure  and  intensity  profile  lnfornatior  on  a  single  pulse  before 
trying  to  compress  it.  We  have  investigated  the  extent  to  which  such  informtlon  Is 
obtainable  using  the  well  established  principles  of  spectroscopy  and  the  more  recent 
two-photon  fluorescence  (TIT)  technique'  '  and  we  are  in  the  process  of  setting 
up  the  equipment  to  make  measurements.  We  believe  that  the  technique  we  are  developing 
will  come  close  to  extracting  the  maximum  available  Information  from  a  single  pulse 
measurement.  A  simplified  and  very  approximate  theory  of  the  measurement  technique 
and  a  description  of  the  experimental  equipment  follow.  When  experimental  results 
are  obtained,  a  more  precise  theory  will  be  developed. 

At  a  fixed  point,  the  optical  fields  of  a  pulse  vary  In  time  as  A(t)e,o(t)  where 
the  principal  term  in  <p(t)  ismQt.  The  rate  of  change  of  the  phase  ls<P(t),  which 
varies  in  time  if  there  is  phase  modulation.  One  could  detect  the  phase  modulation 
by  taking  a  set  of  filters  In  the  same  region  of  space  tuned  to  different  center 
frequencies  "l*  *2’ *3’  *tc ■  the  same  bandwidth  *•,  end  observing  the  sequential 
responses  of  these  filters  as  the  pulse  goes  by.  This  Is  Illustrated  in  Fig.  17  for 
a  complex  phase  modulation  in  a  pulse  with  large  time -bandwidth  product.  The  optimum 
filter  bandwidths  would  depend  on  the  phase  and  amplitude  structure  of  the  pulse,  but 
a  feeling  for  the  nagnitudes  involved  can  be  obtained  simply  by  considering  a  smooth 
pulse  envelope  (e.g.,  Gaussian)  with  a  linear  frequency  chirp.  The  product  of  the 
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pulse  duration  T  and  the  pulse  bandwidth  (Aoj/Stt)  is  denoted  by  N,  i.e.  T^»jj/2tt  *  N. 
Therefore  (T//N)(a<i>/2tt/N)  =  1-  A  filter  of  bandwidth  (W2n/N)  has  a  response  time 
duration  of  approximately  its  reciprocal,  which  in  this  case  is  T//N.  Hence,  by 
dividing  the  bandwidth  into  /N  parts,  the  original  pulse  duration  can  be  resolved 
into  ,/N  filter  responses.  The  curve  cp(t)  versus  t  can  in  this  case  be  drawn  with 
resolvable  points.  At  the  same  time,  the  power  spectrum  versus  tu  'an  be  sketched 
out  with  /N  resolvable  points .  We  do  not  expect  to  be  able  to  do  better  then  this 
with  a  single  pulse.  If  m  identical  pulses  could  be  generated,  the  effective  time- 
bandwidth  product  becomes  mN,  and  the  phase  and  power  spectrum  curves  could  be  drawn 
with  /(mN)  resolvable  points. 

The  construction  of  /Ti  separate  filtering  devices  that  could  be  located  in  the 
same  region  of  space  looks  difficult  if  one  thinks  in  terms  of  multilayer  film  deyices, 
etc.  Our  approach  to  the  problem  is  indicated  in  Figure  (  18).  The  laser  pulse  is 
analysed  in  a  low  rcs.olution  spectrometer,  and  the  relative  arrival  times  of  different 
spectral  components  in  the  focal  plane  is  measured  by  "two-photon  fluorescence" 
techniques.  The  incident  bean,  is  exfxinded  in  a  Galilean  telescope,  then  is  diffracted 
from  a  ruled  groting  with  1200  rulings  per  millimeter.  The  energy  in  the  light  pulse^ 
which  was  distributed  in  a  disk  at  right  angles  to  the  propagation  direction  is 

now  distributed  in  a  volume  that  is  canted  relative  to  t  he  general  propagation 

direction.  The  resolving  power  at  this  point,  for  first  order  diffraction,  equals 
the  number  of  rulings  intercepted  by  the  incident  beam.  The  ctepped  mirror  again 
distributer,  the  pulse  energy  into  a  disk  at  right  angles  to  the  beam  direction, 

although  the  energy  distribution  is  canted  in  any  local  region.  (The  angle  at  which 
the  pulse  is  canted  is  determined  by  the  angular  dispersion  caused  by  the  grating.) 

The  resolving  power  of  the  spectrometer  has  now  been  degraded  to  the  number  of  rulings 
in  the  intercept  of  Figure  (  19),  while  the  angular  dispersion  ic  left  unchanged. 

We  have  constructed  an  apparatus  like  that  shown  in  Figure  (i8  ).  The  stepped 
mirror  has  about  60  elements,  and  wan  made  by  annealing,  grinding,  polishing  and 
tilting  a  stack  of  microscope  slides.  The  filter  bandwidth  of  this  system  at  any 
point  in  the  focal  plane  is  about  20  cm"1-,  and  the  angular  dispersion  in  the  focal 
plan?  is  about  l4  cm"1  per  millimeter. 

At  each  point  in  the  focal  plane  there  is  a  filtering  property  with  center  fre¬ 
quency  that  varies  linearly  across  the  plane  and  with  filter  bandwidth  determined  by 
th?  degraded  resolution.  The  focal  plane  is  in  the  center  of  a  Rhodaraine  60  fluorescent 
dye  coli^  r/  .  One  team  ic  inverted  so  that  the  ^  versus  x  characteristic  is  reversed. 
For  a  linearly  chirped  pulse,  the  wnvepackcts  would  be  slewed  in  the  manner  illustrated 
by  Figure  (  )  as  they  npproach  the  focal  plane.  The  TPF  display  of  such  a  chirped 

pulse  would  be  slewed  as  illustrated  by  the  shaded  region  in  Figure  (20).  Denoting 
th*'  time  of  arrival  of  the  ;pectml  energy  centered  around  frequency  id  by  t(«)  the 
TPF  displays  the  odd  function 


T(a»)  ■  t(».  *  bm)  -  t(«  -  bm)  •  -T(-a») 


(4.1) 
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for  the  phase  structure  and  the  spectral  function 

S(a«j)  =  J1^  +  A®Jt')  l(»^  -  AwJt'Jdt' 


(4.2) 


where  the  simulaneous  functional  dependence  on  frequency  and  time  is  constrained  by 
the  time-bandwidth  product  considerations  mentioned  earlier. 


This  new  technique  for  using  TPF  has  the  great  advantage  of  being  able  to  present 
phase  structure  information  over  a  frequency  band  that  exceeds  the  bandwidth  of  the 
two-photon  absorption.  By  adjusting  the  relative  horizontal  positions  of  the  two 
beams  the  frequency  can  be  placed  in  the  middle  of  the  two- photon  absorption  band. 
Since  we  correlate  oe tween  l(«h  +  d<*»)  an<*  -  4<u)  the  sum  frequency  involved  in 
the  TPF  is  always  (uvj  +  A®)  +  -  A«)  =  2^  as  indicated  in  Figure  (20).  So  even 

when  a®  greatly  exceeds  the  two-photon  absorption  bandwidth,  the  TPF  display  is  still  visi¬ 
ble  and  in  fsct  will  show  an  increasing  contrast  against  the  background  as  a®  increases. 

Figure  (  21 )  shows  the  TPF  display  produced  by  directing  the  output  train  of 
pulses  from  a  mode-locked  neodymium  glass  laser  through  the  apparatus.  The  similarity 
between  this  Figure  and  Figure  (  20 )  Is  preliminary  evidence  that  the  technique 
works.  Our  aim  is  to  optimize  the  technique  and  use  it  for  investigation  of  non¬ 
linear  optical  effects  that  change  the  phase  structure  of  picosecond  pulses. 
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Section  5 

STIMULATED  SCATTERING  AND  NONLINEAR  PROPAGATION 
5-1  Transient  Stimulated  Raman  Scattering 

During  the  semi-annual  report  period,  a  number  of  experiments  on  stimulated 
Raman  scattering  with  picosecond  laser  pulses  have  been  carried  out  in  collaboration 
with  Dr.  R.  L.  Carman,  Dr.  F.  Shimizu,  J.  Reintjes  and  Professor  N.  Bloembergen  of 
Harvard  University*. 

The  essential  difference  between  stimulated  Raman  scattering  with  nanosecond 
duration  pulses  from  a  Q-switched  laser  and  that  produced  with  picosecond  pulses  from 
a  mode-locked  laser  lies  in  the  fact  that  phonons  generated  in  the  scattering  process 
require  time  to  build  up.  The  characteristic  time  associated  with  this  build  up  is 
the  reciprocal  of  the  spontaneous  Raman  linewidth.  In  liquids,  for  those  lines 
which  have  been  stimulated,  the  linewidths  range  from  about  0.1  cm"1  to  20  cm"1, 
giving  a  phonon  build  up  time  of  0.5  to  100  psec.  Consequently,  while  the  "steady 
state"  can  be  achieved  in  stimulated  scattering  with  Q-switched  pulses,  the  scatter¬ 
ing  produced  by  mode-locked  pulses  of  picosecond  duration  must  be  transient. 

More  quantitatively,  the  optical  electric  field,  E(z,t),  produces  in  the  molecular 
vibrators  a  dipole  moment,  p,  proportional  to  the  molecular  polarizability,  a,  that 
is, 

P  *  (5.1) 

In  general,  the  molecular  polarizability  is  a  function  of  the  normal  coordinates, 
q^,  of  the  various  normal  vibrational  modes  of  the  molecules.  Taking  for  simplicity 
only  a  single  mode  with  normal  coordinate,  q,  the  polarizability  may  be  expanded  to 
give 


P 


cr  E  + 
0 


(5-2) 


The  effect  of  this  polarization  on  the  electric  field  can  be  found  by  substituting 
Eq.  (5.2)  into  Maxwell's  equations,  remembering  that  the  bulk  polarization  is  just 
N  times  the  above,  where  N  is  the  molecular  density.  In  a  similar  fashion,  the 
electric  field  also  affects  the  molecular  vibration.  Let  us  assume 


E  -  EL(z,t)et(kl-'-  '  ■*>  .  Es(,,t)el(kS=  *  U»>  .  0. 


(5.3) 


and 


*  The  contributions  of  the  Harvard  group  supported  under  NASA  Grant  NGR22-007-117  and 
Joint  Services  Electronics  Program  Contract  N00014-67-A-0299-0006. 
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=  Qrz,t)et(kPhZ  ‘  ubt)  +  c. 


(5.4) 


where  and  Eg  represent  the  laser  and  Stokes  fields  respectively,  og  =  -  u)Q  and 

kph  =  k^  -  kg.  If  laser  depletion  is  neglected,  Maxwells'  equations  yield 


-3ES 


(5.5) 


where 


_  2ttNuk 

1  ’  A  0 

c  Ks 

Similarly  from  the  equation  of  motion  for  the  molecular  vibration 


(5.6) 


fr +  w  -  Wl*  (5-7) 

where  t  is  now  the  retarded  time,  r  is  the  Raman  linewidth  and 

Y2  =  (BqUq  09  0  (5.8) 

Here  ^  is  the  reduced  mass  for  the  vibration. 

The  coupled  Eqs.  (5*5)  and  (5-7)  can  be  treated  conveniently  in  two  limits.  The 
first  is  the  so  called  steady  state  limit  where  dQ*  Roughly  speaking  this 

ir<<rQ*- 

corresponds  to  the  case  that  the  laser  pulse  duration  is  long  compared  to  the  vibra¬ 
tional  damping  time,  l/l\  Consequently,  the  molecular  vibrations  are  able  to  follow 
the  changes  in  the  optical  fields.  In  this  limit 

Eg(z,t)  ce  exp(-p-^  |Ej2  z)  (5-9) 

The  Stokes  wave  experiences  exponential  gain  with  a  gain  coefficient 

VnYo  o 

s  -  -p2  Iel!  (5.10) 

Since  y^y  «  o,  the  integrated  Raman  cross  section,  only  lines  with  both  a  large 
cross  section  and  a  narrow  Raman  linewidth  will  exhibit  high  gain.  It  is  for  this 
reason  that  nearly  all  of  the  lines  stimulated  with  Q-switched  lasers  represent 
symmetric  stretching  modes. 
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The  transient  limit,  t— -  »  rQ*>  corresponds  approximately  to  the  case  where 

CLX 

the  laser  pulse  is  short  compared  to  the  reciprocal  linewidth.  In  this  case  the  molecular 
vibrations  cannot  follow  the  envelope  of  the  optical  pulses.  The  general  solutions 
to  Eqs.  (5*5)  and  (5. 7)  were  first  given  by  Wang'9°'.  They  have  since  been  discussed 
in  some  detail  by  Carman,  Shimizu,  Wang  and  Bloembergen'99' .  In  the  transient  limit 
it  is  found  that  near  the  peak  of  the  laser  pulse 


Eg(z,t)  *  exp  (y^z 


(5-11) 


The  Stokes  wave  again  experiences  gain  but  the  gain  is  considerably  less  than  the 
under  steady  state  conditions.  This  is  simply  a  reflection  of  the  fact  that  the 
phonon  wave  has  insufficient  time  to  build  up  to  its  steady  state  amplitude. 

Of  particular  importance  is  the  fact  that  in  contrast  to  the  steady  state  result, 
Eq.(5*H)  shows  the  transient  gain  to  be  independent  of  the  Raman  linewidth.  To 
exhibit  high  gain  a  Raman  line  need  only  have  a  large  cross  section.  On  this  basis 
stimulated  Raman  scattering  with  picosecond  laser  pulses  would  be  expected  to  occur 
in  different  lines  and  in  a  wider  variety  of  materials  than  with  Q-switched  laser 
pulses.  This  prediction  has  found  verification  in  our  most  recent  experimental 

work '100,101) t 

Also  of  interest  is  the  unusual  time  and  distance  dependence  indicated  in 
Eq.  (5.1l).  The  Stokes  gain  is  seen  to  reach  a  maximum  after  the  peak  of  the  laser 
pulse.  This  is  a  consequence  of  the  continuing  build  up  of  phonon  wave  amplitude. 

This  continuing  build  up  results  not  only  in  a  delay  in  the  Stokes  pulse  emission  but 
also  results  in  a  Stokes  pulse  width  which  is  critically  dependent  on  the  detailed 
time  evolution  of  the  laser  pulse.  If  the  laser  pulse  has  sharp  leading  and  trailing 
edges,  such  as  for  example,  in  a  square  pulse,  the  Stokes  pulse  will  be  dramatically 
sharpened  as  the  result  of  the  high  gain  occurring  at  the  trailing  edge  of  the  laser 
pulse.  If  the  laser  pulse  has  a  long  precursor  or  a  long  trailing  edge,  the  Stokes 
pulse  can  be  considerably  longer  than  the  laser  pulse.  For  a  Gausi'ian  laser  pulse 
shape,  the  laser  and  Stokes  pulses  should  have  about  the  same  duration^99J . 


5.2  Experimental  Results  -  Liquids 

Transient  stimulated  Raman  scattering  has  been  observed  at  United  Aircraft  in 
both  liquids  and  gases  using  a  mode-locked  ruby  laser.  Attempts  have  also  been  made 
at  the  McKay  Laboratory  by  the  Harvard  group  to  observe  stimulated  Raman  scattering 
with  a  mode-locked  neodymium  laser.  They  have  been  able  to  reproduce  many  of  the 
results  with  the  liquids  using  a  frequency  doubled  beam  but  have  had  no  success  with¬ 
out  frequency  doubling  and  have  been  unable  to  reproduce  any  of  the  results  in  the 
gases . 
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With  the  mode-locked  ruby  laser  stimulated  Raman  scattering  was  observed  in 
all  the  liquids  tried.  These  included  liquids  like  carbon  tetrachloride,  water, 
methanol,  ethanol,  propanol  and  acetic  acid  where  stimulated  scattering  is  obtained 
only  with  difficulty  with  Q-svitched  pulses.  Stimulated  scattering  was  also  obser¬ 
ved  in  the  more  usual  Raman  materials  like  benzene,  chlorobenzene,  carbon  disulfide 
and  acetone.  In  all  cases,  the  scattering  could  be  observed  with  a  collimated  input 
beam.  Sample  cell  lengths  of  25,  50  and  75  cm  were  used. 

Figure  (22)shows  a  laser  pulse  train  and  the  corresponding  Stokes  pulse  train 
f~r  carbon  tetrachloride.  It  is  not  yet  understood  why  the  Stokes  pulse  train 
peaks  early  in  the  laser  train  and  falls  off  rapidly  toward  the  middle  and  end.  The 
forward  to  backward  Stokes  emission  ratio  was  measured  and  was  found  to  be  greater 
than  1000:1,  the  limit  of  measurability  in  that  experiment.  This  confirms  the  travel¬ 
ing  wave  nature  of  the  scattering.  In  addition,  the  gain  in  the  backward  direction 
was  measured  by  reflecting  the  forward  emitted  Stokes  radiation  back  through  the 
liquid  cell.  There  was  negligible  gain.  This  result  indicates  that  any  background 
intensity  upon  which  the  mode-locked  pulses  are  superimposed  must  be  at  least  1000 
times  less  intense  than  the  picosecond  pulses  themselves.  This  information  is  quite 
important  in  studying  traveling  wave  interactions,  such  as  the  present  one. 

Of  all  the  liquids  tested,  carbon  tetrachloride  showed  the  greatest  energy 
conversion.  Up  to  20$  of  the  0.5  J  energy  in  the  laser  pulse  train  was  converted  to 
Stokes  energy.  The  relatively  high  conversion  allowed  a  Stokes  pulse  width  measure¬ 
ment  by  the  two  photon  absorption-fluorescence  technique.  Rhodamine  6C  in  ethanol 
was  found  to  be  a  suitable  material  for  that  purpose.  Figure (22)  shows  the  result. 

It  can  be  seen  that  the  Stokes  pulse  duration  is  comparable  to  the  laser  pulse  dura¬ 
tion  (5-10  psec).  The  fact  that  the  two  pulse  widths  are  abojtgComparable  lends 
support  to  the  notion  that  the  laser  pulse  shape  is  Gaussian^  ' .  The  spectral  width 
of  the  first  Stokes  and  of  the  laser  line  were  also  comparable. 

Spectra  were  taken  of  the  outputs  of  a  number  of  liquids.  All  showed  considerable 
generation  of  Stokes  in  several  orders  but,  curiously,  very  little  anti-Stokes  genera¬ 
tion.  The  Stokes  line  in  carbon  tetrachloride  and  chlorobenzene  was  found  to  be  a  doublet, 
in  contrast  with  the  results  with  Q-switched  lasers  where  a  single  shift  of  459  cm" 
is  reported.  The  doublet  may  arise  because  of  the  two  isotopes  of  the  chlorine.  In 
methanol,  two  separate  lines  are  observed,  one  at  2833  cm"1  corresponding  to  a  C-H 
Stretching  mode  and  one  at  2941  cm’1  corresponding  a  symmetrical  CH_  bending  mode. 

Figure  (23)  illustrates  both  the  spontaneous  and  stimulated  scattering  with  bending 
mode.  The  transient  nature  of  the  scattering  process  contributes  heavily  toward 
allowing  stimulated  scattering  in  this  rather  broad  line  to  be  observed. 

Near  and  far  field  patterns  have  also  been  taken  with  many  of  the  liquids.  The 
far  field  patterns  indicate  a  Stokes  beam  divergence  of  from  3  to  10  mrad  depending 
on  the  liquid  used.  This  is  to  be  compared  with  a  laser  beam  divergence  of  about  1 
mrad.  The  anti-Stokes  far  field  patterns  show  an  emission  principally  in  the  forward 
direction  but  also  in  phase  matched  cones.  The  near  field  photographs  indicate  Stokes 
emission  both  from  the  whole  beam  and  from  small  filaments  in  the  case  of  long  eells. 

With  short  cells  emission  is  observed  only  from  the  filaments,  indicating  that  as  In 
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the  case  with  Q- switched  lasers  the  threshold  for  Ranan  scattering  in  the  liquids  is 
determined  by  self-trapping  and  not  by  the  Ranan  effect  itself.  The  sharp  threshold 
for  the  trapping  and  Stokes  generation  is  shown  very  clearly  in  Figure  (2k)  . 

Dielectric  reflectors  were  used  to  vary  the  laser  intensity  without  changing  the  near 
field  pattern. 

Figure  (25)  shows  self-trapping  in  methanol  as  observed  at  the  laser  wavelength, 
at  the  Stokes  wavelength  and  at  the  anti-Stokes  wavelength.  Although  the  filaaents 
arc  very  pronounced  in  the  Stokes  and  anti-Stokes  photographs  they  do  not  show  up  in 
the  laser  photographs.  Instead  one  observes  a  number  of  small  dark  areas  about  the 
sane  size  as  the  Stokes  filaments,  suggesting  that  once  the  laser  light  is  trapped 
it  is  nearly  completely  converted  to  Stokes  and  anti-Stokes.  The  size  of  the  8tokes 
filaments  varies  widely  depending  on  the  material.  They  range  from  about  a  5  u 
diameter  (resolution  limited)  in  carbon  disulfide,  to  about  a  100  u  diameter  in  water. 
There  is  also  a  considerable  range  in  size  in  any  given  eatcrial.  In  eethanol,  for 
example,  the  Stokes  filaments  run  from  5  -  ko  u.  Curiously,  the  anti-Stokes  fl laments 
vary  much  less  in  size  and  arc  generally  quite  small.  In  eethanol  nearly  all  have 
the  resolution  limited  5  u  diameter.  More  curious  still  is  the  fact  that  each  of  the 
anti-Stokes  filaments  is  surrounded  by  from  1  to  6  concentric  and  nearly  equally 
spaced  rings  of  anti-Stokes  light.  The  diameter  of  the  first  ring  varies  from  fila¬ 
ment  to  filament.  The  formation  of  these  rings  is  believed  to  be  related  to  the  so 
called. class  II  antl-S+okcs  observed  in  eelf- trapping  liquids  with  Q-switched 
lasers1 103>. 

One  interesting  and  practical  result  cf  the  s«!f-xocuclng  study  is  the  fact  that 
focusing  the  laser  beam  inhibits  self- trapping.  This  is  shown  in  Figure  (26)  where 
the  results  with  a  collimated  beam  are  contrasted  with  those  with  the  laser  beam 
focused  by  a  50  cm  focal  length  lens  Into  tie  center  of  the  cell.  The  near  field 
photograph  taken  at  the  end  of  the  cell  shown  little  evidence  of  trapping  In  the  latter 
case.  This  is  true  not  only  for  acetic  acid  but  for  all  of  the  lew  Kerr  constant 
liquids.  Recent  experimental  results  indicate  that  some  filamentation  still  does 
occur  but  it  only  exists  over  a  chort  length  near  the  focus  of  the  lens.  Consequently, 
if  a  high  quality  Stokes  beam  is  required  it  is  best  obtained  by  focusing  the  laser 
into  the  liquid  cell  and  recollinating  afterward.  The  reason  trapping  is  inhibited 
in  the  focused  beam  la  not  completely  clear  at  the  present  time,  but  it  is  likely 
associated  with  the  decrease  in  laser  intensity  per  unit  solid  angle.  Similar  . 
inhibition  of  seif- focus  lng  in  focused  bea-s  has  been  observed  by  other  workers' 


5.3  Experimental  Results  -  Oases 

Recently  a  considerable  effort  has  been  devoted  to  transient  ntimlaied  Raman 
scattering  in  gases.  Transient  conditions  arc  easily  achieved  ty  virtue  of  the  fact 
that  the  linewidthc  are  relatively  narrow,  even  at  the  high  pressures  required  to 
reach  threshold.  To  observe  stimulated  scattering,  a  test  ceil  5**  cm  long  was 
pressurized  with  the  sample  gas.  The  output  from  the  mode-locked  ruby  laser  was  then 
focused  into  the  center  of  the  cell  with  a  50  cm  focal  length  lens.  The  spectra 
of  the  stimulated  scattering  from  gases  in  the  coll  were  recorded  on  Kodak  1R  plates 
using  a  Jorrell-Ash  75-152  Spectrograph.  The  experimental  results  for  stimulated 
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vibrational  scattering  are  summarised  la  Table  I,  which  iuu  all  the  feats  tried. 

The  nu-.be r  and  variety  of  gases  in  which  stimulated  mean  scatterlj*  has  been  observed 
stands  in  strong  contrast  to  previous  work  with  Q- switched  ^asers.  where  successful, 
results  were  obtained  only  with  methane,  hydrogen  and  deuterlua^^**  ®'.  in  other 
gases  competition  with  Brlllouln  scattering  prevented  the  build  up  of  stimulated  Pa  nr 
scattering^  -  •  Transient  excitation  with  picosecond  laser  pulses  not  only  allows 
the  stimulation  of  large  crocs  section,  broad  lines  but  also  very  effectively  dis¬ 
criminates  against  the  more  slowly  responding  Brlllouln  scmtterlig;. 


Of  all  the  gases  examined,  the  largest  conversion  efficiencies  occurred  for  C8^ 
and  SF^.  In  both  cares,  more  than  twenty  prerent  of  the  total  Incident  larer  energy 
was  converted  to  Stokes  energy.  In  at  stipulated  scattering  could  be  produced  without 
focusing.  Stimulated  vibrational  scattering  in  a  collimated  beam  could  also  be 
generated  In  SF^,  I,  and  00,  by  tclescopli*  down  the  Incident  ruby  laser  beam  with  a 
2  power  Galilean  telescope.  Nearly  all  of  the  conversion  lakes  place  In  the  early 
pulses  of  the  train,  as  Is  shown  In  Figures  (27a)  and  (27b)  for  SFf .  Up  to  seventy 
percent  of  the  energy  in  these  Initial  pulses  could  be  converted  to  Stokes  energy  ns 
verified  by  the  depletion  of  the  incident  User  in  Hgurc(27a)  and  by  direct  energy 
measurements.  The  reason  for  t*  decrease  In  conversion  towards  the  end  of  'he  train 
le  not  clear,  but  it  nay  represent  a  saturation  of  the  vibrational  population.  Another 
possibility  is  tha  guild  up  of  slower  thermal  or  acoustic  perturbations  such  pc 
discussed  by  Fohi(1  ,  and  observed  recently  by  one  of  the  present  authors* 1  . 

In  the  cats#  of  C0?,  0^  and  Kn0,  stimulated  pure  rotational  and  rotational -vibra¬ 
tional  ftaean  scattering  could  also  be  produced  if  the  Incident  ruby  beam  Is  circularly 
polar!  :sed  k  .  With  Q-rv  itched.  Infers  .only  hydrogen  and  deuterium  have  shown 

stimulated  rotational  scattering***’  *  In  N.O,  which  has  the  largest  rotational 
cross  section  of  the  three  gnsee,  both  first  and  second  order  rotational  and  rotational- 
vibrational  Stokes  lines  arc  excited.  Fvidence  of  rotational  scattering  le  also  teen 
In  propylene.  Figure  (26)  shews  pure  rotational  ftwan  spectra  and  corresponding  d«*c Ho¬ 
me  ter  traces  for  CO^  and  Cl.  In  00„  the  lines  3(12)  through  S(2b)  have  been  stimulated, 
while  In  Ol,  8(5)  through  if  13)  have  been  stimulated.  The  line  shifts  obrerved  agree 
within  experimental  error  with  those  reported  for  the  spontaneous  spectra**1  t  In  the 
case  of  K^O,  the  narrowness  of  the  line  epaclf*  in  comparison  to  the  llnevldihs  pre¬ 
cludes  an  a  ecu  rated  identification  of  the  stimulated  lines.  The  observed  rotational 
scatterlrg  spans  the  region  from  8(13)  to  8(20)  in  both  orders. 


Often  the  rotational  lines  observed  are  broadened  and  skewed  fra  line  center  or 
even  split  into  several  components.  This  Is  believed  to  be  the  result  of  an  optical 
Stark  effect.  Regiectli^  any  remanent  molecular  dipole  moments,  the  Stark  shift  for 
circularly  polarlred  light  Is  to  lowest  order  in  perturbation  theory 


gW. 


V  vl  J(J  ♦  1)  -  3n2 
3  c  {2J  -  1;(2J  •  3) 


(5.12) 


where  the  s-axle  has  been  chosen  In  the  direction  of  propagation.  Here  *  Is  the  anisotropic 
part  of  the  tolarl-lUIUy  tensor  and  I  is  the  incident  intensity.  Uslng(  ),  ore  -ay 
readily  confirm  that  th*  shifts  are  an  appreciable  fraction  of  the  total  rotvtlonal 
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energy  for  the  J- levels  involved,  and  that  the  splitting  la  the  rotational  shifts 
for  the  different  at- states  Is  of  the  carder  of  a  fev  reciprocal  cer.tiaeters,  la  quali¬ 
tative  nj*f»ct  vlth  the  splitting  observed.  la  cases,  where  rotational  freedom  Is 
preserved,  the  fitark  shift  a  are  also  responsible  for  the  observed  Kerr  effect.  For 
V  >  0  In  Eq.  (5*12),  the  lowest  energy  states  are  a  •  £J,  corrvi;  ponding  to  molecular 
rotation  In  the  plant  of  the  electric  field,  Thlc  Is  equivalent  to  classical  allcment. 
further  evidence  of  the  Importance  of  the  optical  Start  effect  In  the  transient 
regime  Is  provided  by  the  observation  of  self-trapping  In  a  coi  11  sated  bean  In  the  most 
anisotropic  esses,  CO  and  X^O.  This  Is  clearly  seen  in  the  near  field  patterns  shown 
in  Figure  (J?9). 

If  linearly  polarised  light  Is  used  In  an  attempt  to  observe  stimulated  rotational 
n-nn  scattering,  one  tees  Instead  a  broad  wing  shaded  to  the  Stokes  side  on  both  the 
ruby  line  and  on  the  vibration  Stakes  line.  Similar  results  are  found  In  the  case  of 
f*yl«*lgh  vli^  gciu«rliq(ll‘*  .  Explanations  for  the  difference  in  behavior  with  circular 
and  linearly  polarised  light  have  been  given  by  Herman  4 •  ''an!  by  Chlao  and  Coding 
The  seme  of  circular  polarisation  of  the  stimulated  first  rotational  an!  rotat local - 
vltraltacal  lines  Is  always  reversed  from  that  of  the  Incident  laser  light.  The  sane 
situation  Is  found  lr.  the  stimulated  rotational  scattering  from  hydrogen^ and  In 
Xaylelgh  wing  scattering^1*  .  These  experimental  results  are  in  agreement  with 
theoretical  predictions  giving  a  forward  gain  ratio  for  the  opposite  to  saw  sense 
of  polarisation  of  6*1,  If  the  anti -Stokes  coupl  leg  is  Ignored*  108'and  6:0,  If  this 
coupling  Is  taken  Into  account  iL  K  The  second  rota  i local  ani  rotational-vibrational 
ft  tees  lines  observed  In  K^O  are  polarised  in  the  same  sense  as  the  laser  light,  ani, 
therefore,  In  the  opposite*" sense  as  the  first  rotatl  ml  Stokes  lines,  as  expected 
fnw  the  sequential  nature  of  the  higher  ord^r  Stokes  production  (Figure  30). 


5.1  Vibrational  Decay  Tine  Xufurm.ts 

Curing  stimulated  ftiman  scattering  not  only  Is  the  Stokes  wave  amplified  but  so 
also  Is  the  phonon  wave.  For  every  Stakes  photon  produced,  an  excited  vibrational 
quantum  Is  also  created.  Stimulated  Nran  scattering  Is,  therefore,  a  convenient  way 
of  generating  a  non-equlllbrlun  population  distribution.  TV'  short  lim  scale  of  the 
present  experiments  makes  transient  stimulated  ham  an  scattering  a  useful  tool  for  the 
study  of  gas  kinetics.  A  major  area  of  application  Is  In  chemical  laser  diagnostics. 
TV  re  standard  techniques  often  carrot  give  the  desired  relncatlon  rates  under 
conditions  pertinent  to  chemical  laser  operating  conditions.  Shock  tube  analysis  Is 
often  useful  but  high  te  praturee  and  pressures  arc  required.  With  electrical 
excitation  Ionised  species  ard  electrons  are  necessarily  generated.  In  etlisulated 
h  su  scattering  only  the  excited  hamar  active  vltrailoml  stale  Is  populated. 
Experiments  can  be  run  at  row  le-qe-miure  or  below  if  desired.  High  pressure  is 
required  for  the  Pamn  oscillate*’  .*  the  relaxation  rate  Is  Utvsr  with  pressure, 
then  the  results  *ny  be  extrapo  »  -ed.  If  not  linear,  then  a  lew  pressure  fiamar. 
amplifier  mny  be  used  in  conjunction  with  a  high  pressure  Pawn  oscillator.  The  gas 
properties  would  then  be  studied  In  the  amplifier. 


J9£{*79-21 


The  simplest  type  of  Measurement  would  be  a  fluorescence  measurement.  Such  a 
■murrafut  could  be  used  to  determine  the  vibrational  decay  rate  in  HF  as  tfc  vn  in 
Fl£ure{  U) •  Other  candidates  are  the  terminal  laser  levels  in  00.  and  JLO  and  the  first 
vibrational  level  in  CO.  Preliminary  experiments  are  new  under  way  in  nF. 

A  sore  general  though  none  difficult  technique  for  determining  vibrational  decay 
rate  is  shown  in  Flgure($2).  Here  me  monitors  the  excited  state  population  by  detest* 
ir*  spontaneous  anti-Stokes  scattering  produced  by  a  second  high  intensity  light 
source,  possibly  another  laser.  Prior  to  the  arrival  of  the  mode-locked  pulse  train 
and  the  generation  of  stimulated  Stokes  light  little  or  no  spontaneous  anil-Stckee 
light  would  be  observed  due  to  the  absence  of  excited  state  population.  The  anti- 
Stokes  Intensity  vould  build  up  with  the  creation  of  vibrational  population  during  the 
node-locked  pulse  train  and  would  decay  afterward  with  the  decay  rate  to  be  dele  mined. 
This  dwble  Faman  technique  has  been  used  with  a  u- switched  rub v  laser  to  determine 
the  vlbratlocal  decay  time  in  hydrogen  gas  at  room  temperature' '•  with  transient 
ct  inula  tel  Paean  scattering,  this  technique  can  be  used  to  sake  Measurements  In  a 
host  of  different  gases,  including  all  of  those  listed  in  Table  1. 
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Section  6 
OAGAIXC  DTI  lASDtS 

6.1  r  mple,  High  Intensity,  Short  Pulse  FlMkkf 

Electrical  discharge  flsshlsspe  eiptblt  of  witting  intense  light  pulses  of  wry 
short  duration  are  required  in  a  rnitty  of  application^,  inciting  photochemical 
e*  idles  and  laser  excitation.  Special  annular  leapt '“6fil?,i  15)  ^YC  ftXina 
particularly  uttful  in  such  applications  but  often  require  considerable  care  in 
construction  and  are  not  readily  adapted  irradiation  of  Irregularly  shaped  staples. 
Setae  perforaanee  characteristics  art  reported  here  for  high  Intensity  short  pulse 
leaps  which  art  easily  constructed  and  can  be  formed  to  efficiently  irradiate  various 
saaple  configurations. 

A  rnafcer  of  linear  cylindrical  ablating  wall^1^1*  )u*rc  vere  constructed  as 
shown  In  rig j(33aX  using  quarts  tubing  of  1  ■  wall  thickness  for  all  leap  envelopes. 
The  seal  •flexible  plastic  (or  rubber  !  tubes  shown  absorbed  aschanlcal  shocks  associated 
with  the  leap  discharge  and  effectively  prevented  dosage  at  the  ends  of  the  leap. 

The  wire  electrodes  were  soldered  or  spot  welded  to  the  inside  wells  of  the  at  tel 
sleeves,  which  provided  the  neeeeaary  external  electrical  contacts.  Tungsten,  copper, 
alralrca,  end  stainless  steel  electrodes  were  tested  end  ell  were  found  to  give 
slallar  results.  Performance  was  also  essentially  unchanged  when  the  leaps  were  bent 
Into  various  fores  for  irradiation  of  Irregularly  shaped  staples,  ss  In  Fig. (36b) 

To  enhance  the  light  pulse  rlsetias  end  to  minimise  the  Internal  explosive 
shock  associated  with  the  Intense  short  duration  discharges,  the  leaps  were  continu¬ 
ously  evacuated.  A  —11  continuous  sir  leak  was  provided  In  the  pimping  line  to 
salntaln  res  1  dim  1  ma  pressure  high  enough  (typically  1  to  3  torn'  to  ensure  rapid 
reproducible  discharge  Initiation,  but  low  enough  to  minimise  explosive  lamp  failures, 
'fo  significant  changes  were  observed  when  gases  such  ss  argon,  hell  us,  tenon,  nitrogen, 
and  COf  were  substituted  for  sir  In  the  leaps.  This  Is  consistent  with  the  idee 
Uat  the  discharge  was  carried  prlsarlly  be  ablation  products  from  the  quarts  walls 
while  the  res  Id -.ml  low  pressure  gas  served  only  to  facilitate  discharge  Initiation. 

electrical  energy  storage  for  the  lamp  was  provided  by  low  inductance  capacitors 
(Maxwell  laboratories.  Series  M)  which  were  connected  to  the  leaps  through  a  triggered 
sparh  rap  (£33,  y  del  rtp-ike)  as  shown  In  *ig.(3b)  To  minimise  pulaa  rlsetimes,  all 
electrical  connections  In  the  discharge  circuit  were  nde  as  short  ss  possible.  In 
sane  rases  several  leaps  were  connected  in  parallel  with  another,  and  all  were 
discharged  simultaneously  (measured  delay  between  lamps,  less  than  0.05  a  esc)  through 
the  spark  sap. 

inch  lamp  waa  tested  to  destruction  to  determine  the  explosion  energy,  l.e. 

■he  discharge  energy  above  which  the  leap  would  fail  within  one  or  a  few  shuts.  The 
failure  mode  In  these  teste  was  slways  that  of  violent  explosion  near  the  center  of 
the  leap  envelope. 
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After  determination  of  the  explosion  energy  for  a  given  lamp,  an  identical 
lamp  was  fired  at  approximately  seventy  percent  of  the  explosion  energy  and  the 
light  pulse  rise tins,  duration,  and  peak  intensity  ears  recorded,  using  an  8-h 
photodiode  and  an  oscilloscope  with  a  total  lnstnasntal  response  time  shorter  than 
2.5  nsec. 

Operation  of  a  leap  at  energies  higher  than  seventy  percent  of  the  explosion 
energy  generally  resulted  in  increased  light  pulse  fall  times  without  substantially 
increasing  the  peak  intensity.  Kleetlmee  ears  somewhat  shorter  at  energies  near 
the  explosion  limit,  but  lamp  life  «aa  greatly  reduced  from  the  several  hundred 
shots  often  obtained  at  seventy  percent  of  the  explosion  energy. 

htrformanc*  characteristics  typical  of  the  lamps  tested  are  presented  in  Table  11. 
As  expected,  the  light  pulse  rise  times  (lot  to  90£)  were  shortest  for  the  smallest 
capacitors  used,  due  to  their  relatively  low  inductance,  ftleetlmts  were  also 
dependent  to  some  extent  on  lamp  geometry,  generally  increasing  signifleently  for 
long  thin  lamps  whose  Inductance  contributed  excessively  to  overall  circuit  Inductance, 
and  for  large  diameter  lamps  in  which  wall  ablation  apparently  developed  slowly  due 
to  the  large  free  channel  available  to  the  discharge.  Data  for  lamp#  exhibiting 
rlee times  in  excess  of  s  microsecond  ere  not  included  In  Table  1  since  conventional 
gee  filled  lamps  can  often  be  used  in  applications  where  longer  risetimee  are 
permissible. 

Pulse  widths  (between  half  intensity  points)  were  generally  greatest  for  lamps 
of  large  energy  handling  capability  and  tended  to  be  smallest  for  leaps  of  Hall 
diameter  in  irfiich  emitting  atoms  were  apparently  quickly  deexcited  by  wall  collisions, 
task  light  intensity  and  energy  handling  capability  were  found  to  Increase  with 
increasing  lamp  dimensions  and  in  feet  were  quits  simple  related  to  the  projected 
area  (length  times  disaster)  of  the  leaps,  as  shown  in  Pig<(35).  The  simplicity  of 
the  approximate  sapirleal  relations  of  Plg.(35)  makes  it  possible  to  quickly  compare 
energy  handling  capabilities  and  aaxiaua  light  intensities  for  various  leap  configura¬ 
tions.  It  should  be  recognised  however  that  extrapolation  of  these  results  to  lamps 
of  very  large  disaster  or  length  ay  not  be  Justified. 

The  two  righthand  columns  in  Table  I  represent  qantitiss  roughly  indicative 
of  the  total  light  emitted  (the  product  of  half  power  pulse  width  by  peek  light 
intensity)  and  of  the  imte  of  rlee  of  light  intensity  (the  ratio  of  peak  light 
intensity  to  light  pulse  risetime).  The  latter  quantity  any  be  useful  in  characteris¬ 
ing  lamps  designed  for  certain  applications,  such  as  ot^nlc  dye  laser  excitation, 
ytm re  feet  rising  intense  light  pulses  ere  required.  71 

Measurements  using  a  ballistic  thermopile  in  conjunction  with  various  filters 
indicated  ♦hat  the  lamps  used  in  these  tests  radiated  mote  than  twenty  percent  of  the 
electrical  Input  energy  into  the  2000  to  6000  ft  spectral  region,  with  more  then  three 
fourths  of  the  radiation  felling  in  the  2000  to  hOOO  A  range*  These  mrasuremsnts, 
combined  with  the  results  of  several  tests  performed  to  estlsate  leap  emisslvity. 
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•uggest  that  within  this  spectral  range  the  lwp<  salt  essentially  aa  blackbody 
rad  In  tort  with  t  afiratur—  In  the  20,000  to  30,000  °K  ran# e.  The  corresponding 
ste  radiancy  at  peak  intensity  in  tha  2000  to  **000  ft  spectral  ran#*  would  thus  ba 
typically  a  faw  hundred  wattle*2.  ataradlan.  A). 

Aa  a  practical  taat  for  lamps  of  tha  type  described  above,  five  leaps  200  aa 
lory  with  3  aa  bores  wars  placed  alongside  a  10  m  die  by  100  an  Ion#  quarts  call 
containing  a  10**  aolar  solution  of  either  rhodaaine  65  or  7-diethyiamino-l*-aethyl 
co'aarln  In  ethanol.  Broadband  mirrors  of  100  percent  and  33  percent  reflectivity, 
respectively,  were  placed  at  either  end  of  the  cell  and  aliened  perpendicular  to  the 
cell  axis.  When  all  leaps  were  simultaneously  fired  in  parallel  through  the  spark 
rap,  using  a  total  stored  electrical  energy  of  about  700  J,  laser  oscillation  ensued 
at  the  expected  wavelengths  of  about  5700  A  for  the  rhodaaine  and  1*600  A  for  the 
counerln  dye.  Typical  oscillographs  of  the  leap  pulse  and  of  the  counarln  laser  out¬ 
put  are  shown  in  Flg.(j6>  The  laser  output  energy  emitted  through  the  33  percent 
mirror  was  about  one  joule  In  each  case,  being  somewhat  higher  for  the  rhodaaine  dye. 
This  represents  one  of  the  highest  output  energies  reported  to  date  for  organic  dye 
lasers  of  this  type  end  demonstrates  one  practical  application  of  the  above  described 
lamps. 


6.2  Ultra  fast  Flashlampe  for  tye  lasers 

Since  the  development  of  dyt  lasers  two  different  approaches  to  the  problem  of 
pumping  the  dyes  have  been  pursued.  laser  maplng  wss  first  demonstrated  In  Infrared 
emitting  dyes  using  a  Q- switched  ruby  laser*2*  .With  frequency  doubling  this  technique 
wee  then  extended  to  visible  and  near  *Jf  emitting  dyesf  -)  K  re  recently,  the  pulsed 
nitrogen  laser  has  been  used  with  great  success  to  pu*p  %  wide  vs  rifely  of  lyes  over 
the  W,  visible  and  near  1ft  regions  of  the  spvctrvsi  *  .  7tx  fast  rls*tlv»e 

achievable  with  laser  rasping  avoid  the  problem  presented  by  singlet  to  triplet 
crossover  and  reault  In  a  relatively  high  conversion  efficiency.  Fieehlaatp  pumping 
was  first  demonstrated  in  19S7  by  Sorokin  and  laniard  .  The  most  successful 
flashlampe  to  date  ere  those  of  Furomoto  end  'tccon,^*  who  have  achieved  lasing 
in  nny  of  the  d/ee  which  can  be  pumped  with  laser  pumping.  Flashlampe  offer  the 
adrantare  of  broadband  pumping  and  generally  higher  ojtput  energies.  Considering 
that  the  ptstping  la  a  one  step  process,  the  overall  efficiency  la  often  higher. 

The  principal  disadvantage  of  flashlampe  is  that  feet  risetimes  have  not  been 
achieved  so  that  triplet  quenching  limits  the  number  of  dyes  which  can  be  aade  to 
last.  A  ttashlamp  with  a  rise  time  of  5-10  nanoseconds  end  an  energy  handling 
capability  of  a  few  Joules  would  offer  all  the  advantages  of  both  methods  of  pumping. 
Curing  the  seml-ann'jml  report  period  an  experimental  effort  to  devlss  such  a  lamp 
wis  initiated. 

For  the  initial  efforts  storage  In  high  voltage,  barium  tltirmte  "doorknob" 
capacitors  woe  tried.  These  ere  came  only  available  in  a  593  picofarad  capacitance 
end  a  20  or  30  kV  operating  voltage.  The  33  kV  capacitors  have  e  self-resonant 
frtj4»p*y  of  about  50  Ho,  while  the  20  kV  capmcilrs  have  a  60  He  self- resonant 
frequency.  Two  types  of  discharges  wtre  tried,  an  unconfined  discharge  in  air  or 
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nitrogen,  and  a  capillary  discharge  in  nitrogen.  The  results  with  a  teat  fixture 
bolding  three  (3)  of  the  30  kV  capacitors  are  shown  in  Pig.  (37).  These  photographs 
represent  an  input  energy  of  about  3/4  Joule  at  33  kV  in  each  case.  For  the  open 
arc  1.2  ca  disaster  ball  electrodes,  separated  by  1.9  ca  were  used.  The  light  output 
has  a  10  nsec  risetias  and  a  width  at  helf-aaxlmum  of  35  nsec.  This  is  shout  the 
risetlas  expected  considering  the  30  Me  self-resonant  frequency  for  the  test  fixture. 
With  the  capillary  discharge  the  light  output  had  a  risetias  of  18  nsec  and  a  width 
at  half-aaxiatai  of  70  nsec.  The  peak  light  output  is  about  the  saae.  The  arc  vis 
1.9  ca  long  and  confined  to  the  inside  of  a  1  aa  bore  quarts  capillary  tube  with  a 
0.9  >rs  wall  thlckneaa.  A  pressure  of  3  ata  of  Kg  was  necessary  to  hold  off  the  33  kV 
firing  voltage;  however,  the  light  output  and  risetias  were  found  to  be  lhirly 
insensitive  to  pressure  over  a  range  of  fraa  1  to  4  ata.  A  parallel  trigger  was 
used  to  fire  the  laap.  The  slower  risetias  of  the  capillary  lamp  is  believed  due 
to  the  foraatlon  of  a  hotter,  aillsr  dlaaster  and,  therefore,  higher  Inductance  arc. 
The  greater  Integrated  light  output  in  the  capillary  laap  attributes  to  the  hotter 
discharge.  Similar  results  have  been  found  by  other  workers  with  somewhat  slower 
discharges 

The  open  air  arc  was  coupled  with  a  dye  cell  in  an  effort  to  produce  lasing  In 
an  ethanol  solution  of  rhodaalne  6c.  In  fhet,  lasing  was  obtained  even  with  only  a 
39  percent  output  reflector.  However,  the  lasing  was  vary  erratic  due  principally 
to  the  very  erratic  behavior  of  the  open  arc  in  the  presence  of  a  near  by  ground  plane. 
Atteapta  to  correct  this  problem  were  unsuccessful. 

Aa  a  result  of  those  experiments  a  aylar  insulated  30  kV  strip  line  capacitor 
storing  two  (2)  joules  is  under  construction.  The  erratic  behavior  of  the  open 
arc  discharge  in  the  presence  of  ths  dye  cell  mandates  the  use  of  capillary  laaps. 

To  reduce  ths  circuit  Inductance  several  of  these  will  be  used  in  parallel  in  a  scheme 
similar  to  that  in  the  700  Joule  flaahlamp  array  described  in  Section  6.1. 

6.3  Mode- Locked,  Flaahlamp  Pumped  Co-jmarln  Eye  laser  at  4600  £ 

Experlaants  during  this  period  have  resulted  in  what  is  believed  to  be  the  first 
mode- locking  of  a  flaahlamp  pushed  dye  laser  operating  in  the  blue  region  of  the 
•pec t run.  Eyes  suitable  for  passive  node-locking  of  such  lasers  are  difficult  to 
find,  due  to  the  extremely  large  absorption  cross  sections  required.  However,  the  n 
technique  of  intracavity  adulation  at  the  longitudinal  node  difference  frequency 
is  readily  applied  to  node-locking  the  relatively  long  duration,  high  intensity  output 
pulses  typical  of  these  lasers. 

The  iastr  used  in  these  experiments  consisted  of  100  percent  and  33  percent 
reflectors  spaced  by  160  cm,  with  a  20  ca  by  1  ca  die  quarts  cell  centered  between 
the  reflectors  aid  filled  with  a  10*1*  molar  solution  of  7-diethylaalno-4 -methyl 
couaarin  in  ethanol.  The  dye  solution  flowed  through  the  cell  at  a  rate  of  about  one 
liter  per  alnute  and  was  excited  by  aaana  of  five  specially  constructed  linear  flash- 
lam? a  of  the  type  described  in  Section  6.1,  symetricelly  placed  alongside  the  cell. 

The  laaps  were  fired  simultaneously  through  a  triggered  spark  *ap  from  a  3.6  microfarad 
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capacitor  charged  to  20  kV.  The  resulting  flash  was  characterised  by  a  half ‘intensity 
duration  of  a  few  microseconds  and  a  risetime  of  several  hundred  nanoseconds. 

An  acoustically  driven  fused  quarts  block  was  located  in  tbs  laser  cavity  close 
to  the  100  percent  reflector  and  periodically  diffracted  a  few  percent  of  the  laser 
radiation  out  of  the  cavity  by  virtue  of  the  acoustic  standing  wave  in  the  block.  ^ 
When  the  periodicity  of  the  diffraction  was  made  equal  to  the  11  nsec  round  trip 
transit  time  of  radiation  in  the  laser  cavity,  the  4600  ft  laser  output  seen  through 
the  33  percent  reflector  exhibited  the  mode-locked  characteristics  shown  in  Fig«(38). 


6.4  Superradiant  Traveling  Wave  Dye  Laser 

In  recent  years  a  number  of  dye  solutions  have  been  made  to  lase  with  either 
flashlaap  or  laser  pimping.  One  striking  characteristic  of  many  of  these  dyes  is  the 
very  high  gain  which  can  easily  be  obtained.  Neumann  and  Hercher^^have  recently 
demonstrated  lasing  in  a  laser  pumped  rhodamine  6g  solution  a  millimeter  thick  with 
cavity  mirrors  having  a  reflectivity  of  only  a  few  percent.  Experiments  such  as  this 
suggest  that  with  a  solution  thickness  of  a  few  centimeters  and  with  a  somewhat  higher 
pimiping  density,  superradiant  emission  should  be  observed. 


Of  particular  interest  is  the  case  of  pumping  with  picosecond  duration  pulses 
from  a  mode-locked  laser.  The  dye  Inversion  is  then  a  wave  traveling  with  a  velocity 
equal  to  that  for  a  pump  pulse.  This  situation  parallels  that  in  Shipman's  traveling 
wave  nitrogen  laser.  &30)  Superradiant  traveling  wave  emission  has  been  observed  from 
several  polymethlne  cyanine  dyes  pusped  by  a  mode-locked  ruby  laser.  This  technique 
offers  greater  simplicity  and  ease  than  the  more  usual  mode- locking  technique 3 l/f or 
generating  short  dye  laser  pulses. 


The  experimental  apparatus  for  observing  traveling  wave  dye  laser  emission  is 
shown  In  tig^39u).  Excitation  of  the  dye  was  accomplished  with  a  mode-locked  ruby 
laser ^  which  produces  a  train  of  pulses  each  having  a  duration  of  from  2  to  5  psec 
and  a  peak  power  of  up  to  about  5  GW.  A  beam  divergence  of  1  mrad  is  typical  for  this 
laser.  The  test  cell  containing  the  dye  is  2  cm  thick  and  is  wedged  so  that  the  windows 
make  an  angle  of  about  10°  to  each  other.  The  length  of  the  cell  was  chosen  to  be 
sufficiently  long  that  stimulated  emission  could  be  achieved  with  a  moderate  dye 
concentration,  yet  the  cell  was  short  enough  that  stimulated  Raman  scattering  in  the 
solvent  (methanol)  was  negligible.  For  photoelectric  detection  and  for  the  far  field 
photographs,  Coming  CS-7-69  filters  were  used  to  separate  the  dye  emission  from  the 
ruby  pump  light.  At  full  intensity  a  slight  leakage  of  ruby  light  through  these 
filters  does  occur  but  this  is  small  compared  to  the  dye  laser  emission.  The  additional 
filtering  provided  by  the  presence  of  the  dye  in  the  test  cell  reduces  this  leakage 
to  below  detectability. 


A  superradiant  traveling  wave  emission  was  observed  for  three  dyes,  cryptocyanine, 
PDI  (1,  1’ -diethyl-2, 2 '-dicarbocyan me  iodide)  and  DTTC  (3,3’-Diethylthiatricarbocyanine 
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iodide)  etch  dissolved  in  methanol.  The  results  for  DTTC  are  repreaentative.  jifith 
thia  dye  the  emlaaion  occurs  in  a  130  to  180  K  wide  band  centered  at  from  7920  A  to 
8080  %  depending  on  the  dye  concentration.  Maximum  output  occurs  with  a  concentration 
giving  a  band  center  at  7980  £.  At  this  concentration  the  low  level  transmission  at 
6943  A  through  the  2  cm  test  cell  corresponds  to  an  optical  density  of  6.0.  More 
than  90  percent  of  the  incident  pump  light  is  absorbed. 

The  output  from  the  dye  (Fig.  39b  )  occurs  as  a  train  of  pulses  which  follows 
more  or  less  the  overall  shape  of  the  incident  pump  pulse  train.  Far  field  photo¬ 
graphs  such  as  that  in  Fig^39c)ahow  that  the  emission  takes  place  in  a  narrow  beam 
with  an  angular  divergence  of  about  19  mrad.  As  in  the  case  of  pumping  with  a  Q- 
switched  laser,  0-33)  the  direction  of  polarization  of  the  output  beam  is  the  same  as 
that  in  the  pump  beam.  With  1  joule  of  pump  energy  the  total  energy  emitted  in  the 
forward  direction  is  generally  from  10  to  30  mJ.  The  forward  to  backward  emission 
ratio  was  measured  photoelec trically  and  was  found  to  be  about  100:1,  thus  confirming 
the  traveling  wave  nature  of  the  device. 

An  attempt  was  made  to  measure  the  pulse  width  directly  by  the  two  photon 
absorption-fluorescence  technique.  This  was  not  successful.  However,  the  pulse 
width  can  be  inferred  from  other  measurements.  The  high  forward  to  backward  emission 
ratio  indicates  that  the  fluorescence  decay  time  and,  therefore,  the  pulse  duration, 
must  be  less  than  the  90  psec  transit  time  through  the  cell.  One  indication  of  how 
much  less  is  given  by  the  fact  that,  if  the  entrance  window  of  the  dye  cell  is  sufficiently 
thin,  an  emission  occurs  not  only  in  the  forward  direction  but  alos  in  the  direction 
normal  to  the  outside  ffcce  of  this  window.  This  emission  arises  due  to  the  amplifica¬ 
tion  of  fluorescence  radiation  generated  near  the  window,  which  travels  through  the 
window  and  is  reflected  back  at  the  glass  to  air  interface.  The  resultant  emission 
is  shown  in  Fig.(39d), where  a  window  1.9  mm  thick  was  used.  With  a  window  twice  as 
thick,  as  was  the  case  in  Fig.(39d)pnly  a  faint  emission  in  the  direction  of  the 
window  normal  is  observed  even  with  the  entrance  face  nearly  perpendicular  to  the 
punp  beam.  This  indicates  an  appreciable  decay  of  gain  in  a  time  equal  to  the  difference 
in  the  round  trip  transit  times  for  the  two  windows.  That  is  to  say,  the  pulse  duration 
must  be  on  the  order  of  19  psec  or  less. 

That  such  short  pulses  are  to  be  expected  can  easily  be  shown.  In  the  frame  of 
reference  moving  with  the  emitted  pulses  the  growth  of  the  energy  density,  W,  per  unit 
frequency  interval  per  unit  solid  angle  follows  the  relation'  134' 

^£W(z,t)  -  hvS(v)^AN(z,t)  *  hvS(v)AN2(z,t)  (6.1) 

where  S(v)  is  the  lineshape  function,  A  and  3  are  the  Einstein  A  and  B  coefficients 
per  unit  solid  angle  for  given  polarization  and  N2  and  AN  are  the  excited  state 
population  and  inversion  in  the  dye,  respectively.  Saturation  effects  are  ignored  in 
Bq.  (6.1).  In  general  Ng  and  AN  are  functions  both  of  position,  z,  in  the  moving 
frame  and  of  time,  t.  However,  if  the  difference  in  group  velocities  of  the  emitted 
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dye  laser  pulse  and  the  Incident  pump  pulse  is  ignored,  Ng  and  AN  are  not  explicitly 
dependent  on  time  and  Eg.  (6.1)  can  be  readily  integrated. 

The  results  are  shown  in  Fig.(4o)  for  an  excited  state  population  of  the  fora 


N  (tQ)  =  N(e-to/Tf  -e‘to/Tr)  (6.2) 

Here  the  position,  z,  has  been  replaced  by  t0  =  z/v,  where  v  is  the  common  group 
velocity.  For  the  polymethine  cyanine  dyes,  a  fluorescence  decay  time  of  Tf  «  50 psec 
would  seem  to  be  a  reasonable  estimate  from  the  present  experiments.  The  population 
risetime  would  be  the  longer  of  the  pump  pulse  duration  or  the  Franck-Condon  time  ^35) 
for  the  dye.  A  value  of  Tr  *  2  psec,  corresponding  roughly  to  the  pump  pulse  duration, 
was  chosen. 

From  Fig.  21  it  can  be  seen  that  a  fairly  modest  gain  will  give  a  substantial 
pulse  sharpening  initially.  However,  once  the  dye  laser  pulse  duration  becomes 
comparable  to  the  pumping  pulse  duration,  or  more  accurately,  to  the  inversion  risetime, 
a  considerable  increase  in  gain  is  required  to  achieve  even  a  slight  decrease  in 
pulse  duration.  If  the  dye  laser  emission  builds  up  to  the  point  that  the  gain 
saturates,  then  there  may  be  an  additional  pulse  sharpening  due  to  the  nonlinear 
amplification.  0-36) 

Although  this  technique  has  been  demonstrated  only  for  a  single  class  of  dyes 
all  emitting  in  the  infrared  region  of  the  spectrum,  it  should  be  possible  to  obtain 
visible  region  emission  in  other  dyes  by  pumping  with  the  second  harmonic  of  either 
a  mode-locked  ruby  or  neodymium  laser.  However,  the  longer  lifetimes  of  most  of  the 
other  common  laser  dyes  would  necessitate  a  very  high  gain  or  possibly  strong  satura¬ 
tion  in  order  to  obtain  picosecond  pulse  durations. 
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Section  7 


OTHER  AREAS  OF  INVESTIGATION 


7<1  Stimulated  EJnission  from  laser  Produced  Plasmas 


The  possibility  of  obtaining  coherent  radiation  from  laser  produced  plasmas  is 
of  great  Interest  and  seme  preliminary  experimental  and  theoretical  investigations  of 
this  possibility  have  been  carried  out.  The  production  of  ultraviolet  radiation  by 
laser  action  is  hampered  by  the  very  fast  spontaneous  population  decay  which  typifies 
this  region  of  the  spectrum.  For  efficient  excitation,  the  pumping  source  must  have 
a  rise  time  comparable  to  ?r  less  than  the  spontaneous  decay  time.  Conventional 
optical  and  electrical  pumping  sources  are  capable  of  rise  times  in  the  nanosecond 
region.  A  thousandfold  improvement  in  this  figure  can  be  achieved  through  the  use  of 
picosecond  laser  pulses.  The  most  general  scheme  for  pumping  would  be  the  production 
of  a  plasma  by  laser  Induced  breakdown  in  the  material  which  is  to  lase.  The  inversion 
could  then  be  achieved  by  the  same  mechanisms  as  in  the  more  usual  electrical  discharge 
lasers,  although  the  electronic  temperatures  are  much  higher  in  the  present  case. 

Some  general  considerationsQare  treated  here  and  applied  to  the  possibility 
of  obtaining  emission  frem  the  3371  A  line  in  nitrogen  and  from  metal  vapors. 

The  consideration  can  be  extended  to  shorter  wave  lengths.  Sane  experimental 
observations  of  fluorescence  from  metal  vapors  are  also  presented. 


The  gain  coefficient  for  an  inverted  lossless  medium  is 

a  =  huB6teu)  AN 

c  l  »• 

where  B  is  the  Einstein  B  coefficient,  g(&u)  is  a  linewidth  function  and  AN,  the 
inversion.  It  is  convenient  to  define  a  gain  cross  section,  a,  such  that 


«  =  9L  =  huBg(Au) 
Afl  - 5 - 


(7.2) 


In  many  of  the  low  pressure,  high  temperature  superradiant  lasers  which  have  been 
made  to  operate,  doppler  broadening  is  the  most  important  line  broadening  mechanism. 
For  the  purposes  of  discussion,  let  us  temporarily  assume  that  this  is  so  in  our 
case  also.  Then  at  line  center  ^ 

g(4u  .  0)  -  2 jij  (7.3) 

where  AuD  is  the  doppler  width  (FWHM).  This  gives 

(7--> 
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vhere  A  Is  tbs  Einstein  A  coefficient.  This  yields 

Sines  A  is  the  reciprocal  of  tbs  radiative  decay  tiae,  all  of  the  quantities  in 
the  above  relation  are  Measurable. 


In  a  high  density  plasaa  of  the  kind  produced  by  laser  induced  breakdown, 
the  predominant  optical  losses  are  due  to  free-free  absorption  t  37). 

The  loss  coefficient  is  given  by  tbs  celebrated  f omuls 


K  -  3-1  a  10“*1  , 

Tjflw  - 


(7.6) 


vhere  Z  is  the  degree  of  ionisation,  W  is  the  ion  density,  T  is  the  plasaa  taapeta- 
ture,  hu  is  the  photon  energy  in  ev,  and 


g  -  0.55  In  C-Vt0^] 

iV 


(7.7) 


The  total  power  gain  coefficient  is 


•  o-K  -  Hflr  -  3.1  *  10 -31_|Vg. 


T*(h u) 


(7*8) 


vhere  f  -  $■ 


Laser  generated  plasaas  in  pies  typically  have  high  plasaa  tsar  statures 
and  a  high  degree  of  ionisation.  The  values  T  •  5  *  10^  °K  and  Z  •  5  are  represen¬ 
tative  of  the  conditions  which  sight  be  achieved.  The  free-free  losses  for  such 
a  plasaa  are  Indicated  in  FlgJUi)  as  a  function  of  wavelength  and  of  pressure.  It 
has  been  assuaed  that  the  ion  density  is  equal  to  the  particle  density. 

For  the  3371  A  line,  the  doppler  width  is  Aup  **  8.5  x  10*C  sec  **  at 
5  x  105  ^he  radiative  decay  tiae  for  this  line  is  kO  nsec,  giving  A  -  2.5  x 
107  sec*  . '  •  These  values  give  for  the  line  center  91  in  cross  section 

0  -  1.25  x  lO*11  ca2.  (7.9) 
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Tram  this,  tad  Um  ftwftw  iktontion  given  i«  FlgiM),  the  pm  can 

be  calculated  u  s  function  of  the  pimw  end  of  Um  fractional  inversion,  f. 
the  result  1«  shown  la  Figi*2).  Far  l>  10* 3  supenedlsnt  eUsslan  vlU  occur  la 
ft  1  ca  tttii  length.  b*  the  otlar  hud,  if  the  fractional  i eversion  is  less  than 
f  *  l(f  ,  tlw  act  fftla  will  alwa ya  be  negative  so  that  oscillation  win  be 
laposslble.  It  ftiutU  to  noted  that  ftm  our  estiaated  cross  ftcctloa  sad  tto 
■assured  gain  coaffi elect,  tto  fractional  inversion  achieved  la  a  conventional 
nitrogen  disctorie  laser  is  calculated  to  to  f  *  6  a  10’*. 

to  estlsate  gala  for  a  W  lias  we  asstae  that  tto  dipole  a  want  of  tto 
3371  A  lias  la  nitrogen  is  representative.  If  tto  principal  broaden  lag  aectoalan 
Is  doppiec  broadening  (probably  not  a  very  goodLassuapt ion),  than  tto  gala  cross 
section  Is  tto  mm,  naasly  •  •  1.1)  a  10’*  cm  .  Tto  aost  laportaat  effect  of 
decrees  lag  tto  wavelength  is  to  decrease  tto  free-free  atoorptloa  (see  Fig.  hi ). 
tto  aetgda  coefficient  vs  wavelength  and  pressure  Is  shown  la  Fig. (bj) for 
f  •  10  .  tto  peak  gala  la  seen  to  laersase  aad  save  to  higher  pressure  with 
decreasing  wavelength. 


la  recent  years  superradlaat  aad  salon  has  been  observed  la  aetal  vapors 
such  as  those  of  lead,  copper,  tine,  eadalta  aad  aaagaaeaa.  Excited  import  of 
auch  astals  can  to  produced  by  ftouaalag  a  laser  been  onto  s  surface  of  the  solid. 
To  our  knowledge  s  definitive  ate (dr  of  tto  eoadltlona  esl sting  In  tueh.a  plaaaa 
has  base  nade  only  for  carbon.  1  Tto  pertinent  paiaastere  are  T  a  10?  °K, 

If  a  l**9  sod  3  <  2  <  6.  For  want  of  bettor  data,  we  ehaU  take  T  •  KT  °K, 

II  •  lGr9,  2  •  5~for" the  aetols.  Tto  overall  pin  coefficient  Is  then 


0  •  1019  lb  -  1.1  a  lO1^*.  <7.10) 

For  copper  vapor  A  •  1.3  a  1C6  sec*1  and  g  •  ft  i  1C*1* m?  at  10*  °K.  Similarly 
for  lead  A  •  2.2  x  1C r  esc*1  and  r  •  7.2  *  lO*1"  at  l(r  °K.  tto  net  pita  coefficient 
for  copper  sad  leed  are  plotted  as  s  function  of  f  aad  1  in  FigaU%*)sndfl»^>  Again  a 
fraction  Inversion  of  10**  to  10**  Is  needed  for  eupenadlant  auction. 

Tto  calculations  outlined  shove  art  very  rough.  nevertheless,  tto  results 
are  sufficiently  encouraging  that  an  experiaental  effort  has  been  undertaken  to 
observe  stipulated  suasion  la  laser  produced  p Issue,  tto  experiasntal  apparatus 
la  shown  la  Flg4*6).  2s  tto  Initial  expertaants,  tto  node -locked  r toy  laser  baaa  was 
focussed  with  e  cylindrical  lens  onto  tto  target  end  no  resonator  was  used,  tto 
uaplea  Included  copper,  lead,  tine  aad  eadhtita.  tto  visible  region  aUssion  was 
observed  with  e  tost  S-h  photodiode  looking  down  the  axis  of  tto  cylindrical  plasu. 
It  was  hoped  that  the  ^ln  Ught  be  hl0>  to  observe  superradisnt  enissian.  Althoi^h 
such  was  not  the  case,  a  maker  of  interesting  results  were  obtained. 
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A  ItfOH  fluOrUMMM  «M  Obaorvud  fOT  *U  Of  U«  —pill.  Dtt  fMUlU  fOT 

fipr  art  represaeUtiva  (rtf .  A?  ).  Tbt  fluortattaot  tnhaa  plana  in  a  wa*ar  of  llnaa 
including  tba  5105  I  and  57*  t  liaaa  which  art  oapanudlaat  la  a  pulaad  dlactorfo. 

The  ovtmU  fltortMtaot  talaoloa  bacim  naar  tba  ptak  of  tbo  and*  Inched  n*y  lattr 
NlM  train  and  build*  op  until  tbt  tad  of  tba  train  at  vhJeh  tint  it  tboa  dtcajra 
with  an  AO  note  dotty  tint.  Two  okitnatlont  art  of  particular  laportaatt.  Tbt  flrtt 
to  that  the  fluortootnot  rlaao  with  a  o^boaaototoad  (dlodo  llal tod)  tint  constant 
after  each  laotr  pul**,  ladloaUaf  Hot  tbt  fhat  puaplag  rot to  can  to  actuated  by 
Uilo  technique.  Tto  atoned  oboonatlon  la  ttot  overall  flaorttctact  oalaoton  lncrtaaoo 
with  Increasing  power  danalty  at  tto  taryet  with  a  fined  ruby  later  output.  This 
Indloatta  ttot  tto  rdby  later  ebould  to  foci  oood  into  at  Mali  a  volant  at  poet  I  bit. 

It  1*  planned  ttot  further  ext*  river  a  in  this  ir»o  will  bt  confect**.  ttot* 
experiment*  have  U^o  deferred  awaiting  tto  conniption  of  a  10  Joult,  autnaao—cood 
out**  rtodyalua  laeer  facility  that  U  under  construction  at  Uil*  laboratory. 


7.  V  m  i rear  Profanation  EfTtcto 

Wn  on  Intense  lifht  wove  propagate  through  a  Kerr  active  liquid,  ihr 
Ir’cnslty  dependent  index  of  refraction  can  lead  to  a  variety  of  effect*  ouch  a* 
re  f.frbiirr,  **tf-ptof*  modulation  and  steepening  and  alteration  of  tto  iuu 
Priori  ’niton  of  tto  wave.  Burin*  tbl*  period*  <  oe  Id*  rut  Ion  to*  toen  given  to 
latter  effect  and  *everal  conclusion*  haw  torn  reactod.  It  tot  been  xn  wn  for 
tl«e  that  vton  an  elllptleally  polarised  wav*  propagate*  through  each  a  liquid* 
orientation  of  tto  axlt  of  tto  ellipse  are  routed  by  an  avaunt  proportional  to 
product  of  tto  field*  alo*  the  vnjor  and  tto  minor  ant*  and  tto  d I* lance  traveled 
(Inference  1  7)  .Tbl*  effect  to*  Wn  u*ed  to  «ea«rur-  •  m*  of  tto  nail  near  coefficient* 
of  tto  liquid  (deference  1JS).A  xcnrvnat  nor*  general  approach  to  Ul*  effect  to* 
torn  carried  out  ini  It  ha*  been  conclud'd  ttot  llrsar  polar  I  ration  It  an  on*  table 
«ode  of  propagation  In  nucti  a  liquid  while  circular  poiarl'utl  a  I*  n  neutrally 
« t*bi*  •ole.  Tto  natyiU  to*  al»o  be**,  extended  to  tto  case  wtor*  an  addition* 
linear  Hrefrlndetvce,  such  a*  could  be  proiu -*d  by  an  ent^mal  electric  field,  I* 
pre-ent.  Til*  !nw* ligation  I*  *tlll  In  prflgrrszj  an  outlln*  of  tto  method  of 
a (e roach  and  nv  of  tto  aaln  rveuit*  1*  given  tele*.  A  vote  deulled  dl*cu*«lon 
will  be  reported  upon  completion  of  tto  saaiyt!*. 

V*  consider  tto  profanation  of  an  optical  pul**  in  n  liquid  having  an 
anisotropic  molecular  polarizability.  For  tto  present*  v*  will  assume  that  tto 
molecular  reorlenUlton  time,  v*  it  *h  rt  ctmfnred  to  tto  duration  of  tto  |Ml*e. 

F>r  typical  liquids  such  a*  CS_  and  CJLW)  ,  v  1?  a  few  plcooeeond*  to  about  dO 
pleoeecord  .  Xn  this  case  we  rind  ttot\to  polar!  niton  In  tto  liquid  can  be 
written  a* 

f  .  Kill  •  \  •  *IB  •  %L  ,I*>U« 


*2 
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for?  i*se  first  vent  rfprttfut*  Ur  Tittmrjr  Umr  fOiarUabillty,  U*  second  t*rrt 
at  intercity  dependent  bait  Uotrople  poUrlrxfcil  lty  .  tblri  tent  an  anisotropic 

FoUrl’abtllty.  Ve  write  the  optical  field  as  %  fe  '  *  *®4  *licv  for  the 

pr-tetre  of  a  4c  field  It.  TV  d  l  net  loo  of  profnfnttoft  will  be  taken  •*  the  -r 
axis  and  the  direction  of  the  dc  field  will  define  Ur  •*  axis.  The  anisotropic 
fart  of  the  poiart'Abilltjr  my  then  be  expressed  m 


t 

a 


(7.12) 


with  a  #  t  ant  b  •  Hit  *  •  fJL  *)•  TV  field*  set*  to  calculated  from  Ur  ««**; 
•  /ly  «r*«c  envelope  approx Itm^ijn, 


if 

d~ 


'?.  ■ 


Substitution  of  Fqc.  (7*11)  and  (7.12)  Into  Rq.  (7*13)  yields  a  eel  of  two  cwpW 

equations  In  f  and  t  .  It  Is  convenient  to  write 
*  T 

id  U) 

a  . e  „  *»* 

'x*y  «t,oy  (7.1*) 

I  r-wUr^  the  wal  and  1  aginary  forte  of  Rq,  (7.13)  fives  four  elutions  ft  r 
9  ,(  pi  and  d  .  One  tntegrn  L  follows  Iwdlttely,  l.e.,  |!  |  •  |!  |  * 

JC  OX 

cum  tar.  i  as  would  be  expected  since  the  Material  le  lossless.  V*  *ny  Urn  write 

*« •  ,J~» 

l  *  Aim  (7.1*) 

oy 

Finally,  solving  for  f  *  d  •  d  we  find  that  Ur  Isotropic  fart  of  Uv  tunllivar 

y  x 

poknrl'ntrn  drops  out  and  we  obtain  two  equation  for  a  and  f . 


da 

67. 


«  d  *  vX  cos  a  sin  a  sin  at 


(7.i6) 


•  i  *  v  •  vl(cos  Sf  -  l)coc  2* 

v>  o 

vfrr e  v  •  -p  ,  v  «•  I? 

o  c 

These  two  equations  describe  completely  the  state  of  polarl ration  of  tlr  vnvr 
a?  It  i  ropafates.  The  inotropic  part  of  Ur  nor.Uimr  response  doer,  not  enter 
tenure  It  does  not  change  the  state  of  polarisation. 


rations (7. 16)  do  Ml  lead  to  a  direct  solution  for  •(£)  wt  t(£)hoMever,  It  U 
possible  to  find  a  solution  f or  *  ae  a  function  of  a.  We  find 


as 


*.) 


•to  . 


»)  | 


2 

sin  3U| 


•1 


(7.17) 


Here  and  f  nr*  oonstsau  that  nr*  deter. lord  by  the  Initial  conditions  of  polari¬ 
sation.  This  equation  yields  trajectories  of  t  v*  a,  the  sense  of  traversal  of  ther 
being  determined  by  Equation  (7-16).  These  are  plotted  schematically  in  tt*ira(l44for 

the  case  •  0  and  for  v0  l  0.  for  the  oase  vQ  •  0,  the  curve*  art  the 

tv  weenl  routing  polarisation  ellipse.  The  points  sin  t  •  H,  a  •  t  n/%  are  the 
tv.  »<tiec«  of  circular  polarisation.  They  are  neutrally  suble,  as  they  are  surrounded 
ly  the  trajectories.  Th?  case  of  linear  polarisation,  f  •  0,  a  arbitrary  is  seen  to 
be  unstable  as  alt  trajectories  diverge  froa  this  liar.  The  case  v(  /  0  corresponds 
to  the  presence  of  a  linear  birefringence.  Here  the  trajectories  split  into  two 
classes  depending  on  the  initial  conditions  ini  the  Intensity.  One  type  is  slnilar 
to  the  case  vQ  •  0  and  corresponds  to  s  rotating  (and  distorted)  elliptical  polar  It*- 
ti  *.  TV  other  c  ret  spools,  for  example,  to  an  Initial  linear  polarisation  at  an 
angle  q.  The  polarisation  then  becomes  elliptical,  rotates  and  becomes  more  circular. 
After  passing  through  the  circular  state,  th»  polarisation  continues  to  rotate  with 
Increasing  elliptic  It/  until  it  te  cones  linear  at  -a.  In  the  degenerate  case  1*0, 
these  trajectories  simply  becose  a  •  constant,  t  •  v JL  as  My  be  seen  directly  from 
Eqs.  (7.16).  This  is  Just  the  induced  birefringence  and  the  phsee  difference  between 
x  and  y  components  increases  linearly  with  distance.  TV  trajectories  on  a  plot  of 
t  vs  *  are  two  vertical  lines  at  l  I. 

It  Is  anticipated  that  further  wort  trill  be  carried  out  In  the  analysis  of  these 
polarisation  effects.  Analytical  solutions  for  a(Z)  and  #(£)  have  been  obtained  and 
will  be  discussed  in  n  later  report.  It  would  also  be  of  Inure# t  to  extend  this 
annlyels  to  nonlinear  ate  rials  with  s  norm  complicated  spatial  syne  try,  l.e.,  non¬ 
linear  crystals  such  as  oalciu.  It  Is  expected  that  these  nonlinear  polarisation 
effecu  will  be  important  In  the  propagation  of  high  pow»r  pulses  In  dielectric 
•rdla,  particularly  In  the  polarisation  effect#  in  self-trapped  fllawr.U. 
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Table  II  -  Performance  Characteristics  for  Various  ' vr,p3 
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DIAGRAM  OF  AREA  THEOREM  FOR  ULTRASHORT  PULSES 
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FIG.  1 


SCHEMATIC  REPRESENTATION  FOR  THE  BAECKLUNO  TRANSFORMATION 

GIVEN  BY  EQ.  (2-37) 
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BREAKUP  OF  6*r  PULSE  AS  DERIVED  FROM  SEQUENCE 
OF  TRANSFORMATIONS  SHOWN  IN  FIG.  4 
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FIG.  11 


PROPAGATION  OF  n  -  PULSE  IN  AMPLIFIER 
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FIC.  14 


SOURCE  MOVING  ALONG  OPTIC  AXIS 
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SOURCE  MOVING  TRANSVERSE  TO  OPTIC  AXIS 
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FIG.  II 


G.  19 


FIG.  20 


CANTED  TPF  DISPLAY  OF  A  PULSE  WITH  POSITIVE 
LINEAR  FREQUENCY  SWEEP 
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FIG.  21 


TWO-PHOTON  ABSORPTION-FLUORESCENCE  DISPLAY  OF  CHIRP 


~ H  h- 

20  p*ec 


J920479— 21 


STIMULATED  RAMAN  SCATTERING  IN  CCI4 


a)  INTERLEAVED  STOKES  AND  LASER  PULSE 
TRAINS.THE  FIRST  PULSE  IN  EACH  PAIR  IS 
THE  LASER  PULSE 


b)  TWO-PHOTON  ABSORPTION-FLUORESCENCE 
STOKES  PULSE  WIDTH  MEASUREMENT 

c)  TWO-PHOTON  ABSORPTION-FLUORESCENCE 
STOKES  PULSE  WIDTH  MEASUREMENT 
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FIG.  23 


RAMAN  SCATTERING  IN  METHANOL 


a)  SPONTANEOUS  SPECTRUM 

b)  STIMULATED  SPECTRUM  IN  THE  BENDING  MODE 
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FIG.  24 


SELF-TRAPPING  IN  METHANOL  -  STOKES  LIGHT 
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FIG.  25 


SELF-TRAPPING  IN  METHANOL 


a)  LASER 


b)  STOKES 


e)  ANTI-STOKES 


SELF-TRAPPING  IN  ACETIC  ACID-STOKES  LIGHT 


a)  COLLIMATED  BEAM 


b)  FOCUSED  BEAM 


STIMULATED  VIBRATIONAL  SCATTERING  IN  SF6 


a)  LASER  -  INCIDENT  AND  TRANSMITTED 


b)  STOKES 
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FIG.  29 


SELF-TRAPPING  IN  N20  AND  C02 


a)  N20  -  LASER  LIGHT 


b)  N20  -  STOKES  LIGHT 


e)  C02  -  STOKES  LIGHT 
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FIG.  30 


POLARIZATION  OF  ROTATIONAL  SCATTERING  IN  N2O 
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VIBRATIONAL  OECAY  MEASUREMENT  -  GENERAL  TECHNIQUE 
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FIG.  32 


PULSED  LIGHT  SOURCE 
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FIG.  33 


b)  LOW  INDUCTANCE  LAMP 
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CIRCUIT  FOR  RELIABLE  LAMP  TRIGGERING 
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FIG.  36 


o)  LAMP  PULSE 


FAST-RISE  DYE  LASER  FLASHLAMP 
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a)  OPEN  ARC 


b)  CONFINED  ARC 
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OUTPUT  OF  MODE-LOCKED  COUMARIN  DYE  LASER 


SWEEP  SPEED:  100  n»#c/div 


OUTPUT  PULSE  DURATION  VS  PEAK  TRAVEL!  NG  WAVE  GAINS 


FIG.  40 


pstc 


ABSORPTION  COEFFICIENT  K  (cm" 


J92C479-21 


FREE-FREE  ABSORPTION  VS  PRESSURE  AND  WAVELENGTH 
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COPPER  VAPOR  FLUORESCENCE 
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FIG.  4i 


BEHAVIOR  OF  POLARIZATION 
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